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1. Über Stabilität und Turbulem 
von Flüssigkeitsströmen; 
von Werner Heisenberg. 


‘ 


; 


Einleitung. 

Das Turbulenzproblem, das ganz allgemein den 
der folgenden Untersuchungen bilden soll, ist im Laufe der — 
Zeit in so vielen Arbeiten von so vielen verschiedenen Ge- — 
Miehtspunkten aus behandelt worden, daß es nicht unsere Ab- — 
sieht sein kann, einleitend über die bisherigen Resultate eine 
Übersicht zu geben. Wir verweisen zu diesem Zweck auf eine 
Arbeit von Noether!) über den heutigen Stand des Turbulenz- | 
problems, in welcher auch die meisten Literaturangaben zu 
finden sind. 

Für unseren Zweck genügt es, den gegenwärtigen Stand 
des Turbulenzproblems in ganz groben Umrissen zu skizzieren: 
die bisherigen Untersuchungen zerfallen in zwei Teile: die 
@inen von ihnen befassen sich mit der Stabilitätsuntersuchung — 
itgendwelcher laminaren Bewegung, die anderen mit der turbu- 
lenten Bewegung selbst. 

Die ersteren führten anfangs zu dem negativen Resultat, daß 
alle untersuchten Laminarbewegungen stabil seien. v. Mises?) 
und Hopf?) bewiesen auf Grund eines Ansatzes von Sommer- | 
feld‘) die Stabilität des der Couetteschen Anordnung ent- 


sprechenden linearen Geschwindigkeitsprofils, Blumenthal®) _ 
gelangte bei einem von Noether zur Diskussion gestellten _ 
Profile 8. Grades zu demselben Ergebnis. Dagegen gelang 


es später Noether‘), ein labiles Profil anzugeben — allerdings _ 


1) F. Noether, Zeitschr. f. angew. Math. u. Mech. 1. S. 125, 1921. 

2) R. v. Mises, Beitrag z. Oszillationsprobl.: Heinr. Weber - Fest- 
Schrift. 1912. S. 252. 

3) L. Hopf, Ann. d. Phys, 44, 8.1. 1914, ns 

4) A.Sommerfeld, Atti d. IV. congr. int. dei Mathem, Rom 1909. 

5) O. Blumenthal, Sitzungsber. d. ag en Akad. d. Wiss. S. 563. 
1913. 

6) F. Noether, Nachr. d. Ges. d. Wiss. Göttingen 1917. 
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ein Profil, das schon im Falle der reibungslosen Flüssigkeit 
labil ist und auch in Wirklichkeit nie als stationärer Bewegung. 
zustand realisierbar ist. In neuerer Zeit aber hat Prandtli) 
gezeigt, daß es in der Tat Profile gibt, die erst bei Berück. 
sichtigung der Reibung labilen Charakter besitzen. 


Die andere Gruppe von Arbeiten, die in jüngster Zeit 
durch die Rechnungen?) von v. Karman, Latzko u. a. große 
Erfolge erzielt hat, untersucht die turbulente Bewegung selbst 
und geht dabei auf halbempirischem Wege mit Benutzung 
der Ähnlichkeitsgesetze vor. Theoretisch fußt sie fast durch 
weg auf dem Boden der Prandtlschen Grenzschichttheorie, 
Das für uns wichtigste Resultat dieser Arbeiten ist das aus 
dem empirischen Blasiusschen Widerstandsgesetz (die Zitate 
bei Schiller®)) folgende sogenannte y''-Gesetz der turbulenten 
Geschwindigkeitsverteilung. 


Ein Hauptziel der erstgenannten Arbeiten, der Stabilitäts- 
untersuchungen, bildete stets die Ermittelung der kritischen 
Reynoldsschen Zahl. Eine befriedigende Berechnung dieser 
Zahl ist bisher nicht gelungen und es muß auch als frag 
lich betrachtet werden, ob eine solche durch Stabilitäte 
untersuchungen geleistet werden kann. Die Versuche von 
Ekman’), Ruckes*) und Schiller’) zusammen mit dem 
negativen Ergebnis Hopfs beim linearen Geschwindigkeits- 
profil legen vielmehr den Gedanken nahe, daß die kritische 
Reynoldssche Zahl nicht den Punkt angibt, an dem die 
Laminarbewegung anfängt labil zu werden, sondern den, an 
dem zum erstenmal die turbulente Bewegung als stationärer 
Zustand möglich ist. Vom Standpunkt der Theorie aus 
müssen wir also darauf gefaßt sein, unter Umständen zwei 
kritische Reynoldssche Zahlen zu finden, die eine entsprechend 


1) L. Prandtl, Physik. Zeitschr. 28. S. 19. 1922 und O. Tiet- 
jens, Dissert. Göttingen 1922. 

2) Vgl. bes. Zeitschr. f. angew. Math. u. Mech. 1. S. 233f. 191. 

3) V. Ekman, Turbulent motions of Liquids. Arch. f. Mat. och 
fysik 6. S. 12. 1919. 

4) W. Ruckes, Dissert. Würzburg 1907. Vgl. auch den Vortrag 
von W. Wien, Über turbul. Bewegungen. Phys. Zeitschr. 8. 1904 und 
-Verh. d. deutsch. phys.. Gesellsch. 9. 1907. 

5) L. und kritische Suhl. 


Physik. Zeitschr. & 
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Über Stabilität und Turbulenz von Flüssigkeitsströmen. 


dem Beginn der Turbulenz, die andere entsprechend dem 
Niederbrechen der Laminarbewegung. | 

Auch die vorliegende Untersuchung soll in die zwei ver- 
schiedenen Teile, die Behandlung des Stabilitätsproblems 
einerseits, die der turbulenten Bewegung andererseits zerfallen. * 

Das Ziel des ersten Teiles soll sein, alle bisherigen Unter- _ 
suchungen unter einem einheitlichen Gesichtspunkt zusammen- 
zufassen, d.h. möglichst allgemein die Bedingungen dafür 
aufzustellen, wann ein Profil labilen oder stabilen Charakter 
trägt und die Methoden anzugeben zur Lösung der Stabilitäts- 
gleichung für jedes beliebige Geschwindigkeitsprofil und zur — 
Berechnung der für labile Profile sich ergebenden kritischen 
Reynoldsschen Zahl. Dieses Ziel kann natürlich nur durch 
Anwendung von Näherungsverfahren unvollkommen erreicht 
werden. Dennoch hoffen wir, die qualitativ wesentlichen — 
Gesichtspunkte durch solche Rechnungen klarstellen zu 
können. Die Untersuchung irgend eines beliebigen Profils 
scheint zwar zunächst physikalisch sinnlos zu sein, da nur — 
bestimmte Profile in Wirklichkeit vorkommen. Indem wir x 
aber irgendein Profil als endliche Störung eines anderen auf- hy 
fassen können, wie dies z.B. Noether a. a. O. getan hat, — 
wir andererseits später die Untersuchungen auf das zunächst ; 


müssen, so scheint die Untersuchung eines beliebigen Profils 
doch von großer Wichtigkeit. 

Als Anwendung der Methoden soll das Parabelprofil 
durchgerechnet werden. 


stark idealisierenden Voraussetzungen Differentialgleichungen 
für die turbulenten Bewegungen abzuleiten und aus diesen 
qualitative Aufschlüsse zu erhalten über einige Eigenschaften ‘ 
turbulenten Geschwindigkeitsverteilung. 
I. Teil: Die Stabilititsgleichung. 
§ 1. Festlegung des mathematischen Problems. 


Die wesentlichste Beschränkung, die wir unseren Rech- 
nungen auferlegen, besteht in der ausschließlichen Betrachtung 
von ebenen Laminarbewegungen und nur ebenen Störungen 
derselben. Unter Zugrundelegung eines rechtwinkligen K- 
ordinatensystems X, Y, Z nehmen wir also an, daß die Ge- 
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schwindigkeit in der Z-Richtung Null und alle übrigen Größen 
von 2 unabhängig seien. Weiter aber werden wir auch nur 
die Stabilität solcher Laminarbewegungen prüfen, die zwischen 
zwei geraden, parallelen Wänden verlaufen. Dabei nehmen 
wir an, daß die Wände der X-Achse parallel laufen, daß also 
auch die zu untersuchende Laminarbewegung eine Geschwin- 
digkeitskomponente nur in der X-Richtung besitzt. Diese 
Geschwindigkeit w in der X-Richtung wird in irgendeiner 
Weise von y abhängen; über die Funktion w = w(y) behalten 
wir uns vor, später einige Annahmen über Stetigkeit, Sym- 
metrie usw. zu machen; sonst soll diese Funktion jedoch zu- 
nächst ganz willkürlich sein. 

Würden wir w=ay setzen, so würden unsere Ansätze 
gerade mit den von Hopf beim Couetteschen Fall unter. 
suchten identisch. Auf die Frage, ob sich die untersuchten 
Profile w = w (y) als stationäre Bewegungen realisieren lassen, 
gehen wir zunächst nicht ein (vgl. u. $. 581). 

Bevor wir die (bereits von Sommerfeld a. a. O. auf 
gestellte) Stabilitätsgleichung noch einmal in Kürze aus den 
Stokesschen Differentialgleichungen herleiten, führen wir in 
bekannter Weise dimensionslose Variable ein. Sei h eine 
charakteristische Länge (etwa der Abstand beider Wände), 
U eine charakteristische Geschwindigkeit des Profils, mu die 


Zähigkeit, o die Dichte, nn BR R die Reynoldssche Zahl, 
so führen wir statt z, y, u, v, t, p (unter u, v werde die Ge 
schwindigkeit in der z- bzw. y-Richtung verstanden, t sei die 
Zeit, p der Druck) neue Varibale 2, Yo, Up, Vp, to, Po ein ver 
möge der Relationen 


580 W. Heisenberg. 


Lassen wir nachträglich den Index 0 weg, so lauten die 
Stokesschen Gleichungen: 


Ou du du 1 öp 

dv dv 1 Op av 


Weil wir Inkompressibilität voraussetzen, schreiben wir 


dy 
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Durch Elimination von p erhalten wir bekanntlich: 


a ey Ö Oy 


Unter A ist hier das Differentiationssymbol = = 


verstanden. 


Problem, die einer gewissen Laminar- 
stromung, enthalten. Dementsprechend wird Gleichung (4) auch 
die Grundlage für die Rechnungen des Teils II bilden. eo 
speziell zur Stabilitätsuntersuchung überzugehen, teilen wir die 
Bewegung und damit auch das Vektorpotential y in eine 
Grundströmung und darüber überlagerte kleine Schwingungen. 


Wir machen also den Ansatz: BEN. 
= wy) = 0. 


Glieder weg, die p nicht enthalten (da wir r die Gleichung (4) 
fir p= 0 als erfüllt ansehen), lassen ferner alle in g quadra- 
tischen Glieder weg (da wir als klein voraussetzen), so lautet 
die entsprechende Differentialgleichung für 


+ ag) 


trachten, die entweder vermöge äußerer Kräfte wirklich a; 
stationär sind, oder deren zeitliche Veränderung langsam 
gegen die der kleinen Schwingungen erfolgt. a a 
Gleichung (7) ist in dieser Allgemeinheit schon de ee 
Noether a. a. O. abgeleitet. Sie ist eine gewöhnliche Diffe- 
rentialgleichung für 9 von der vierten Ordnung. Dem ent- 
spricht es, daß die Funktion g vier Grenzbedingungen er- 


beiden Wänden verschwinden. Setzen wir B/a =c, so daBc 
im wesentlichen die Wellengeschwindigkeit bedeutet, so läßt 
sich das mathematische Problem folgendermaßen formulieren: 
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Es sind die Lösungen der Gleichung $. 
(Ta) a? gy) (w c) gw" _ 242 9” + eg) 


mit der Nebenbedingung, daß an den begrenzenden Wänden 
(etwa y = 1 und y = — 1) g = 0 und ¢’= 0 sein soll, zu unter. 
suchen. Zu jedem Wert von « und R soll der zugehörige Wert 
von c und ß berechnet werden; @ sei der Einfachheit halber 
stets positiv. Je nachdem der imaginäre Teil von ß positiv, 
Null oder negativ ist, haben wir es mit einer stabilen, unge- 
dämpften oder labilen Schwingung zu tun. Es sind die Be. 
dingungen für das Profil w aufzusuchen, unter denen Glei- 
chung (7a) nur stabile, bzw. unter denen sie auch labile 
Schwingungen zuläßt. 

Bevor wir zu den Lösungsmethoden übergehen, sei noch 
auf eine besondere Eigenschaft der Gleichung (7a) hingewiesen, 
Im Limes der reibungslosen Flüssigkeit R= oo verwandelt 
sich Gleichung (7a) in eine Differentialgleichung 2. Grades 
für 9. 

(8) — pw —c)— pu" =0. 

Dem entspricht es, daß für R=w nur noch zwei Grenz- 
bedingungen zu befriedigen sind, welche besagen, daß die 
Normalkomponente der Geschwindigkeit, also » oder p an 
beiden Wänden verschwinden soll, nicht mehr aber 9’. 

Die Bedingungen für die Lösbarkeit von Gleichung (8) 
sind schon von Rayleigh!) ausführlich untersucht worden. 
Man kann, um eine einfache Bezeichnung einzuführen, 
„schwingungsfähige“ und „nichtschwingungsfähige‘“‘ Grund- 
strömungen unterscheiden, je nachdem Gleichung (8) eine 
Lösung mit reellem c besitzt, die den Grenzbedingungen ge- 
nügt, oder nicht.) Gibt es Lösungen mit komplexem c, so wird, 
wie sich später zeigen wird, auch bei Berücksichtigung der 
Reibung die Stabilitätsfrage für diese Schwingungen schon 
durch (8) entschieden, die Schwingungen sind dann immer 
labil. 

Man wird aber darüber hinausgehend zu der Vermutung 
geführt, daß das Profil w dann und nur dann unter Einfluß 


1) Lord Rayleigh, Papers I. S. 361; III. S. 575, 594; IV. S. 208. 
2) Hierbei genügt es aber keineswegs, das Profil durch Tangenten- 
polygone anzunähern, das Resultat der Scheme 
gungen würde dadurch völlig gefälscht. 
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der Reibung labile oder ungedämpfte Schwingungen zulassen 
wird, wenn es zu den schwingungsfähigen Grundströmungen 
ört. 

” Diese Vermutung liegt um so niher, als sie sich an allen 
bisher untersuchten Profilen bestätigt hat.!) Trotzdem wird 
sie keineswegs dadurch begründet, daß Gleichung (8) im 
Limes R = oo aus Gleichung (7a) hervorgeht, denn es ist z.B. 
in den Arbeiten von Oseen?) bewiesen, daß der Grenzübergang 
R=oo in den Differentialgleichungen schon mehrfach zu 
falschen Resultaten, besonders hinsichtlich der Grenzbedin- 
gungen der reibungslosen Flüssigkeit, geführt hat und daß 
man daher den Grenzübergang nur an den Integralen der 
Gleichung (7 a) vornehmen darf. Außerdem ist es von vorne- _ 
herein keineswegs zu entscheiden, ob die Reibung die un- 
gedämpften Schwingungen von (8) im Sinne einer Dämpfung 
oder einer Anfachung abändert. 

Wir werden im folgenden versuchen, unsere oben aus- 
gesprochene Vermutung zu beweisen, indem wir zeigen, daß die 
schwingungsfähigen Systeme oberhalb eines bestimmten Wertes 
der Reynoldsschen Zahl im allgemeinen labilen Charakter, alle 
michtschwingungsfähigen dagegen stabilen Charakter besitzen. 

Durch diesen Satz wird das Problem der Stabilität eines 
Profils ganz wesentlich vereinfacht; denn für sehr kleine Werte 
von « lassen sich bekanntlich die Lösungen von (8) direkt 
hinschreiben. 


$ 2. Die Lösungsmethoden und das allgemeine Verhalten 
der Integrale von (7a). 

Die wichtigste Eigenschaft der Gleichung (7a), die eine 
näherungsweise Darstellung ihrer Lösungen erlaubt, besteht 
darin, daß R als sehr groß betrachtet werden kann. Es wird 
sich nämlich zeigen, daß, falls eine Labilitätsgrenze existiert, 
diese im allgemeinen bei sehr hohen Werten von R liegt. Da 
és andererseits physikalisch ganz unwahrscheinlich ist, daß bei — 
kleinen Werten von R noch einmal Labilität des betreffenden 
Profils eintreten kann, so reicht es für unsere Zwecke zunächst + 
aus, R als sehr groß anzusehen. 


1) Vgl. auch L. Prandtl, Phys. Zeitschr. a. a. O. 
2) Vgl. z.B. C. W. Oseen, rg 2. tr Ann. d. Phys. 46. 
8.231 u. 623. 1915. 
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Diese Annahme macht es möglich, die Lösungen von 
(7a) anzunähern durch Entwicklung nach negativen Potenzen 
von R bzw., wie sich zeigen wird, von Va R. Ferner werden 
wir « als klein annehmen und die Lösungen gegebenen Falles 
nach positiven Potenzen von «* entwickeln. 

Die beiden Entwicklungsmethoden nach (« R)-"* einer. 
seits, nach «? andererseits, scheinen sich insofern zu wider. 
sprechen, als im ersten Falle «FR als groß, im zweiten a! 
als klein angenommen wird. Der Widerspruch wird jedoch 
dadurch behoben, daß R im allgemeinen als außerordentlie; 
groß betrachtet werden kann, so daß z.B. für R = 2000, 
a = 1/15; R= 200, a? = 1/19) wird, was für eine gute Kon 
vergenz beider Entwicklungen völlig hinreicht. Die Konvergenz- 
eigenschaften dieser Näherungsmethoden sind aber noch ge- 
nauer zu betrachten. Die Untersuchung lehrt, daß die Reihen 
nach (« R)-"* im allgemeinen divergieren, daß sie aber den 
bekannten Charakter der semikonvergenten Reihen tragen, 
d.h. daß die Glieder zuerst abnehmen, dann wieder wachsen 
und daß man die bestmögliche Annäherung erhält, wenn man 
die Reihe beim kleinsten Glied abbricht. Unser Näherungs- 
verfahren hat also ähnliche Konvergenzeigenschaften, wie z. B. 
die in der Astronomie verwendeten Reihen der Störung» 
theorie, deren Verhalten Poincaré, Meth. nouv. d. |. mee, 
cel. II ausführlich schildert. 

Eine erhebliche Schwierigkeit für die Verwendung der 
semikonvergenten Entwicklungen liegt aber darin, daß sie in 
der Umgebung eines bestimmten Punktes ihre Gültigkeit ver- 
lieren, so daß es nicht ohne weiteres entschieden werden kann, 
in welcher Weise die Näherungslösungen auf beiden Seiten 
des Punktes aneinandergefügt werden müssen, um ein be 
stimmtes Integral der Gleichung (7a) beiderseits anzunähern. 
Wir werden hierauf in $ 3 ausführlich eingehen. 

Die Entwicklung nach positiven Potenzen von «? scheint 
im allgemeinen wirklich zu kon: ergieren. Für spezielle Pro 
file läßt sich diese streng beweisen (z. B. für das lineare Profil); 
eine Untersuchung über die Frage, unter welchen Bedingungen 
für das Profil w diese Konvergenz wirklich stattfindet, haben 
wir jedoch nicht ausgeführt. 

Wir beginnen mit der Ableitung der a 


der Gleichung (7a): 
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(1a) cE (p"” — 20? + 
Zu diesem Zwecke setzen wir zunächst 


Die Entwicklung wollen wir auf die beiden höchsten Glieder 
in Y&x R beschränken. Es folgt: 


a (w —c) + Va R (g9' + 29% 91): (w— =ia 
+4 Va R (49° 9, + 697 go) - 


Hierbei ist vorausgesetzt, daB «? und w” von der Größen- a 


ordnung 1 oder jedenfalls < Ya R seien. Es ergibt sich jetzt ER 
durch einfache Rechnung: 


(10) =) —t(w—o), Jandy =- 


Wir erhalten also zwei partikuläre Integrale der Gleichung (7a): 


Der Punkt y, soll dadurch festgelegt werden, daß für y = yp ® 
w=c sei. y, kann also unter Umständen auch komplex sein. _ 
Das Vorzeichen der Wurzel soll so gewählt werden, daß für > ee 


Warzel wird: 
(ce 


Bemerkenswert an diesen beiden Integralen ist, daß a2 
ing in dieser Näherung nicht vorkommt (d.h. nur in der 
Kombination « R, die man, wie aus Gleichung (7a) hervor- — 
geht, in gewissem Sinne als die eigentliche Reynoldssche Zahl 
bezeichnen könnte). 

Wie wir später sehen werden, bestimmen die beiden Inte- __ 
grale (11) das Verhalten von @ in der Grenzschieht und das 
Nichtvorkommen von «? in (11) bedeutet physikalisch, daß ir 
aur solche Schwingungen betrachten, deren Wellenlänge groB 
ist gegen die Grenzschichtdicke — was bei den empirisch 
beobachteten labilen Schwingungen sicher der Fall ist. 
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586 

Zu einem vollständigen Lösungssystem von (7a) brauchen 
wir jedoch noch zwei weitere Integrale; naturgemäß werden | be 
wir dazu diejenigen nehmen, die durch Entwicklung nach | wm 
Potenzen von (« R)-! aus den Integralen der Gleichung (8) | ki: 


hervorgehen. spi 

Zu diesem Ende lösen wir zunächst Gleichung (8) durch 
Entwicklung nach Potenzen von «?. Wir setzen also: int 
(8) — a? (w — pw" =0, 
(12) 
Daraus folgt: ve 
gp” (w —c) — = 0. gr 
(0 — — = (wo), 
ga (2) an’ dei 

Es ergeben sich FIR die Methode der Variation der Kon- 
stanten die beiden Integrale: ch 


(1 


+...) 


Diese Integrale sind nun noch um Größen der Ordnungf Fli 
(« R)-1, ... usw. zu korrigieren, wenn sie der Gleichung (7a) f sie 
genügen sollen. zei 

Ohne die entsprechende Reihenentwicklung im einzelnen 
anzuschreiben, geben wir als Resultat m mit den Größen der f „2 


Ordnung (« R)-! an: vo 
_ dy Po 

(14) 
2 dy Pr 


Mit (11) und (14) ist ein vollständiges —— der f sek 
Gleichung (7a) näherungsweise gefunden. 
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Ehe wir dieses Lösungssystem zur Erfüllung der Grenz- 
bedingungen anwenden, wird es zweckmäßig sein, wenn wir 
ws den physikalischen Sinn der vier Integrale 9,, 92, 93, Pa 
klar machen und einige Resultate vorwegnehmen, die wir erst 
später begründen können. 

Die Integrale y, 9, haben die Eigenschaft, bei den uns 
interessierenden großen Werten von R sehr rasch veränderlich 
zu sein, wie aus dem Exponenten der Ordnung Y« R hervorgeht. 
Wenn also etwa p, an der einen Wand von der Größenordnung 1 
ist, so wird es in einigem Abstand von der Wand exponentiell 
verschwinden (an sich könnte es auch noch außerordentlich 
groß werden, doch dies wird naturgemäß durch die Grenz- 
bedingungen verhindert). Das hat zur Folge, daß sich, von 
der unmittelbaren Umgebung der Wände abgesehen, p aus 93 
und 9, allein zusammensetzt, d.h. dem Verhalten von @ in 
der reibungslosen Flüssigkeit sehr ähnlich ist. 

Daß «* in 9, 9, nicht explizit vorkommt (vgl. Glei- 
chung (11)), wohl aber in 9, 9, (vgl. Gleichung (14)), muß 
physikalisch offenbar so gedeutet werden, daß bei der An- 
nahme, a? sei nicht sehr groß (a? <ya R, vgl. Gleichung (7a) 
und (9)) die Wellenlänge zwar als unendlich groß gegen die 
Grenzschichtdicke, nicht aber gegen die Kanalbreite betrachtet 
werden kann. Der charakteristische Unterschied zwischen den 
,Grenzschichtintegralen“ 9, einerseits, den der reibungslosen 
Flüssigkeit entsprechenden Integralen %,, , andererseits, spricht 
sich also in dem Auftreten bzw. Nichtauftreten von «? be- 
zeichnend aus. 

Was die Konvergenz der Entwicklung nach Potenzen von 
a? betrifft, so können wir hoffen, daß sie für Werte von «? 
von der Größenordnung 1 noch sehr wohl hinreicht, um eine 
gute Näherung zu ermöglichen, denn für ein lineares Profil 
gehen die Reihen für y, 9, über in solche vom Typus der 
Potenzentwicklung von Cosa, die in der Umgebung « = 1 
noch sehr rasch konvergiert. 

Das Strömungsbild, das wir nach allen diesen Schlüssen 
erwarten müssen, entspricht den Ansätzen, die man in der 
Prandtlschen Grenzschichttheorie macht. Mit Ausnahme der 
unmittelbaren Umgebung der Wände gehorcht die Bewegung 
sehr nahe den Differentialgleichungen der reibungslosen Flüssig- 
keiten. An den Wänden selbst aber haftet eine Grenzschicht, 
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deren Dicke von der Größenordnung (« R)-" ist. In diese} _ 
nimmt die Geschwindigkeit u gegen die Wand zu rasch gegen (15 


Null ab, während » auch außerhalb der Grenzschicht schon} 4 
nahezu Null ist. Int 
8 3. Die Übergangssubstitutionen. 
Wenn wir den Verlauf der Integrale von (7a) von de die 
einen begrenzenden Wand bis zur anderen studieren wollen, 
so müssen wir damit rechnen, daß an einer Stelle y = y im (16 
Kanal w—c =O ist (oder wenigstens der reelle Teil vo 
w —c Null ist), daß also die Wellengeschwindigkeit dort mit} ent 
der Geschwindigkeit der Grundstrémung übereinstimmt. In 
diesem Punkte verlieren die Näherungsformeln (11) und (14) 
für die Integrale von (7a) ihre Gültigkeit. 
Es ist daher nötig, die Übergangssubstitutionen für 
1, Ps, Ys, 9a zu kennen, die beim Übergang von Re 
>0 zu Re (w — c) < 0 anzuwenden sind. Zu diesem Ende lin 
entwickeln wir w und » in der ,,Umgebung des kritischen nith 
Punktes y, nach Potenzen von (« R)-’* und setzen deshalb <a 
Y — Yo =: (« R)-*. Ferner nehmen wir an, daß der imag- 
naire Teil von y, von kleinerer Größenordnung als (« R)-"* sei, 
Ist er nämlich von höherer Größenordnung, so verstehen sich 
die Ubergangssubstitutionen von selbst, weil dann im ganze daf 
Bereich reeller y nirgends ein „kritischer Punkt‘ auftritt. Is Gle 
er von derselben Größenordnung, so läßt sich das Verhalten 
von p und w leicht aus den beiden eben angeführten Grenz 
fällen interpolieren. Wir können zunächst sogar Im (y,) =! 
setzen, weil sich » in unserem Falle (Im (y,) < (x R)-"*) nach 
Potenzen von Im (y,) entwickeln läßt und fürs erste nur das} für 
Verhalten von p für Im (y,) = 0 gebraucht wird. 
Wir setzen also jetzt 9 te 
= (a R)-"*; 


Dann entsteht aus (7a): on 
= egy’ +..)an te ge’ by? 
Also in erster Näherung: 
(15) i Por" n a; 
in zweiter Näherung: 
0 Für die Integrale 9,, 9, Gl. (11) entnehmen wir ul 
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(15), daß sie sich im „kritischen Gebiete“ (n Größen- 
ordnung 1) wie die von Hopf beim linearen Profile gefundenen 
Integrale, d. h. wie gewisse Zylinderfunktionen verhalten. Wir 
können also jetzt schon schließen, daß die Übergangssubstitu- 
tionen für gy, 95 avs (11) bis auf Größen der Ordnung (a By 
dieselben, wie beim linearen Profil sein müssen: EEE 


(16) A> 

entsprechen einem Ubergang von icine 

Ne (w c) <0 > Re (w c) >0. 


Für das Studium der Integrale 9, 9, in der Umgebung 
yon Yy—Y,=0 reichen aber die bisherigen einfachen Rech- 
nungen nicht aus, denn für diese würde die Näherungs- 
lösung (15) lauten 9” = 0; aus (14) wissen wir aber, daß im 
Limes R>o0 9,' an der Stelle y — y, = 0 im allgemeinen loga- 
rithmisch unendlich wird. Die Gleichungen (15) und (15a) sind 
also in dieser Form ungeeignet zur Erfassung dieser Singularität. 

Wir setzen jetzt vielmehr « =0 und w’’=0 (d.h. wir 
brechen die Entwicklung von w mit dem 2. Gliede ab), inte- 
grieren aber sonst (7a) zunächst exakt. Dabei bemerken wir, 
daB ¢ = w — c ein partikuläres Integral dieser vereinfachten 
Gleichung (7a) sein muß, und für g machen wir daher den 
aus der Theorie der linearen Differentialgleichungen geläufigen 


Ansatz: = (w — c) f y dy. Dann folgt aus 
= —iaR (p" (w — —w" g) 


= d iA Wa Re 

was nach nochmaliger Integration übergeht in: 2 

(17) (w—c)+8y' w'+ = —ia R((w—c)? yp — (). 


C ist eine Integrationskonstante. Führt man jetzt wieder ein: 
n=(y — yo) (a R)'*, e= (a R)-"*, w—c = ean + &*b7?, 


Y= Yo 
80 ergibt sich: 


Yo + a= — (a? 7? — C) 
(Ta) ! + 8y,' a+ yo" by? + bn + bd 
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Diese Differentialgleichungen enthalten natürlich imme 
noch sämtliche Lösungen von (7a). Da wir speziell 9, studieren 
wollen (p, verhält sich für « = 0 im Punkte y = y, regulär), 
so suchen wir diejenige Lösung von (17a) aus, die sich in 
einigem Abstand von 4, also für große (w —c)a R, wie 


wu verhält; denn wir wissen ja aus (14), daß 9, in einigem 
Abstand von %, durch 


d 
— 6) 


gegeben ist. Es Es wird also nach (17a) 


und 
(18) Wy" 1+ 3y, =~ Cn + any) - 


y, ist wieder dadurch völlig bestimmt, daß es sich „im 


Unendlichen“ wie — Pr verhalten soll. 

Wenn wir nun nach den Übergangssubstitutionen für g, 
(und 9,) fragen, so ist damit folgendes gemeint: In den asym- 
ptotisch (14) für 9, 9, kommt imme 
das Integral Ir wos VO: das seinen Sinn verliert, wenn & 
über die Stelle y=y, (w — c=0) hinweg zu erstrecken ist, 
In Wirklichkeit tritt aber eben an Stelle von wer in der 


Nähe des kritischen Punktes die Funktion y. Es kommt also 
alles auf das Verhalten von y (insbesondere y,) in der kritischen 
Umgebung an. Ist dieses Verhalten, also auch die Grok 
des Integrals f ydy (erstreckt über die kritische Stelle hin 
weg) bekannt, so ist diese Kenntnis gleichbedeutend mit det 
Kenntnis der Übergangssubstitutionen für 9; 9. 


Die durch (18) und die Grenzbedingung im Unendliche 
charakterisierte Lösung y, lautet: 


n 
3(— tet, Hi) dy 
+00 
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Hierin treten die Hankelschen Zylinderfunktionen vom 
Index $ und dem Argument (— ia, 7)"* auf (a, = Dabei 
ist das Vorzeichen von (— 4 &,n)":so zu nehmen, daß (— 4 a, N)": 


ai 
für = — positiv wird. Eine nähere Untersuchung von 


(19) lehrt, daß sich y, in der ganzen oberen Halbebene und 
sogar noch in der unteren Halbebene teilweise, nämlich 
fürn=re'® in den Grenzen 


im Unendlichen wie — aa verhält, wenn «, oder a positiv 


ist. Ist a negativ, so vertauschen sich obere und untere Halb- 
ebene und es gilt 


(20) lim y, (re!) = — wenn [ <E< 
r—>0o an 6 
Daraus entnehmen wer das wichtige Resultat a oe 
+00 te a>odO, 
(21) fe dn = 
_ 200 , 


Es sind also jetzt die Übergangssubstitutionen für 9, 94 bis 
auf Größen der Ordnung (x R)-" genau auch für endliche Werte 
von « gefunden; denn wir wissen jetzt — und das ist hin- 


richend — was wir nach (21) unter dem Integral f wir 
estreckt von w — c < 0 nach w —c >0, zu verstehen haben. 

Auch die Formeln (16) lassen sich aus (17) noch einmal 
ableiten. Denn den asymptotischen Lösungen (11) der Glei- 
chung (7a) wn die Integrale 


(19a) (#( —it, 


der homogenen Gliisheng (18) (C=0). Die Aufgabe, die 
Ubergangssubstitutionen der „asymptotischen‘‘ Integrale (11) 
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und (14) zu finden, ist hiermit in der für uns erforderlichen 
Genauigkeit (bis auf Größen der Ordnung (« R)-"*) gelöst, 


§ 4. Erfüllung der Grensbedingungen und die Stabilität der 
dem Lösungssystem I entsprechenden Schwingungen, 
Nachdem unsere bisherigen Betrachtungen von der Art 
des Profils ganz unabhängig waren — bis auf einige, die singe 
lären Punkte betreffenden Einschränkungen, die dem Profi 
auferlegt werden mußten — wollen wir, um uns nicht in einer 
übergroßen Anzahl verschiedener Möglichkeiten zu verlieren, 
die Natur der Grundströmung weiter spezialisieren. Die Be 
trachtungen gelten jedoch viel allgemeiner. Wir nehmen also 
an, daß die begrenzenden Wände durch die Formeln y = +1; 
y = —1 dargestellt sind, daß ferner die Wand y=1 gega 
die andere eine Relativgeschwindigkeit in der positiven X-Rich- 
tung besitze (von der Größe w (+ 1) —w(—1)) und daß die 
Laminarströmung an den Wänden hafte (was der Versuchs 
anordnung von Couette entspricht); daß schließlich im Ge 
biete — 1 < y < +1, also in der Flüssigkeit, Re (w — c) ein 
mal und nur einmal Null sei. AuBerdem werden wir im ganze 
Gebiete Stetigkeit für w und die Differentialquotienten von u 
voraussetzen und darüber hinausgehend noch annehmen, daß 
die Funktionen w, w’, w’’ usw. stets normale Größe besitzen, 
d.h. z.B. nicht an einzelnen Stellen groß werden von de 
Ordnung (« R)'*. 

Ferner betrachten wir zunächst fiir die folgenden Rech 
nungen « als so klein und « R als so groß, med = mit ge 
nügender Genauigkeit setzen können u el 


9, =w-c, 


(4a) 
a) Y, = (w 


Die Festlegung der unteren Grenze des Integrals in g, bedeute 
offenbar keine Beschränkung der Allgemeinheit unserer Lé 
sungen. Wir setzen dadurch vielmehr 9, als diejenige linear 
Kombination von 9, und 9, fest, die an der Stelle y = -I 
verschwindet. Sollte dort auch w — c verschwinden, so trill 
offenbar an Stelle von 9, die Funktion gy, = w — c, die jets 
für y= — 1 Null wird. 

Um die Grenzbedingungen zuerst an der Wand y=-| 
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zu befriedigen, bilden wir zwei Aggregate f,, f, aus 9, 9a, 
$s, 9a, für die wirklich y=9'=0 ist für 


1 


— Gs 


- 


Tw (— 1) 

i | Ps | : 
w(-1)-e 9(-1) 
Hierbei verstehen wir unter g, von jetzt ab die Wurzel 
-iaR(w— ce), nicht, wie in (10), Y—i(w—c), um uns 

das Anschreiben des Faktors Ya R zu ersparen. 

Um die Grenzbedingungen auch an der anderen Wand 
zu erfüllen, muß man zwei Konstanten A und B so zu bestimmen 
suchen, daß Ah (+1) + Bh (+1) =0, 

A fy’ (+1) + Bf (+1) =9. 
Die Bedingung fiir die Méglichkeit einer solchen Bestimmung 
lautet: 


(23) 


(22) 


fh (+1) fe (+1) 
(+1) fe (+1) 
Durch diese Bedingung wird c bzw. ß festgelegt, wenn R und « 
gegeben sind. Es handelt sich also jetzt darum, die Glei- 
chung (23) nach c aufzulösen und das Vorzeichen des imagi- 
nären Teils von $ zu bestimmen. Gleichung (23) bildet das 
vollstindige Analogon zur Sommerfeldschen Turbulenz- 

gleichung für das lineare Profil. 
Aus (16) entnehmen wir: Ber te 


| (+1) _ 
? 


= 0, 


w(-1)-e 9(-)) 
1 


1){w(— 1) — (~1) 


| Gs (+ 1) (+1) + (+ >| 
u(-D-e 9(-1)- 
Annalen der Physik. 1V. Folge. 74. 


(24) 


h(+ = 1)+ 
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fy (+ )= Ps (+ ) 9 
— |, 


1) (= 1) = + 


| Ps (+ 1) 


v(-1)-e 


%2(- 1) 


“| Diese Werte von f,, fa, fı', fe. setzen wir in (23) ein, nehmen 
aber zuerst noch eine Abschätzung der Größe der einzelnen 
Glieder vor, um die Rechnung nicht durch Anschreiben un- 
wesentlicher Glieder unnötig zu komplizieren. Zu diesem Zwecke 
bemerken wir, daß im allgemeinen entweder 9, (+ 1) >, (—1) 
oder 9 (+1)<&9,(— 1) sein wird. Dies bewirkt der Faktor 


VaR im ete von 9, a in (11), wenn nicht gerade 


sein sollte, was wir ausschlieBen. 

Welcher der beiden Fälle eintritt, läßt sich nicht vorher 
entscheiden; im allgemeinen sind beide möglich und liefern 
beide Lösungen. Im Falle eines schiefsymmetrischen Profil 
gibt der eine Fall die zu dem des anderen symmetrischen Lö- 
sungen. Jedenfalls verhalten sich die beiden Möglichkeiten in 
prinzipieller Hinsicht ganz analog und es genügt daher, eine 
von beiden zu untersuchen. Wir nehmen also an: 


%2(+1)<2(— 1), 


d. h. (vgl. S. 585) der Punkt w 
als an w = w (— 1) liegen. 


daß ,(—1) ~ 


= c soll näher an w = w (+) 


klein, also —— iD sehr groß ist. Es bleiben daher in ji 


und f,’ nur die Glieder stehen, die 9, (— 1) im Nenner haben, 
in f, und fs fallen die Glieder, die , am, fort. 


Aus (28) erhalten wir so: hob f 
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(+ 1) (+ 1)) 


(+1) 


[py (+ 1) + ig, (+ 1)] 


1 
— 


(25) 


Auch in dieser Form ist die Gleichung für c noch recht kompli- 
gert. Wir vereinfachen daher (25) zunächst weiter, indem wir 
jetzt nicht nur Größen der Ordnung e- Var , sondern auch 
solche der Ordnung (« R)-"* streichen. 

Zu diesem Zwecke stellen wir fest, daß g, (+ 1) von der 
Ordnung (« R)'* ist — wir schließen also zunächst aus, daß 
w(+1) — c sehr klein ist —, daß ferner: 


+ (+1 =- + ig, (+ 


+ 9o(+ 1) 


Von den Gliedern der Gleichung (25) behalten wir daher nur 
noch diejenigen bei, die mit dem Faktor gu (+ 1) multipli- 
ziert sind. 

So entsteht das einfache Resultat: 


od 
d 
(26) e y — 2 => - 


Diese Gleichung besitzt zwei völlig verschiedene Lösungs- 
systeme: 


+1 
L e Yo I. f w- 
-1 


Das System I stellt das vollkommene Analogon dar zu den 
Lösungen, die Hopf für das lineare Profil erhalten und a. a. O. 
§ 4 ausführlich diskutiert hat. In der Tat zeigt sich, daß die 
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dem System I entsprechenden ne stets stabilen 


27) 2 f dy = +30, 
Yo 


won eine positive (vgl. 8. 585) ganze Zahl ist. Es ist leicht zu 
sehen, daß diese Gleichung nur erfüllt sein kann, wenn «c=} 
einen positiv imaginären Teil besitzt. Also sind die durch (9) 
charakterisierten Schwingungen tatsächlich gedämpft, der Be 
trag der Dämpfung ist von der Größenordnung w (+ 1) -t 
und braucht daher keineswegs klein zu sein. 


Sr § 5. Das Lösungssystem II und die Bedingungen für die 
Labilität eines Profils. 


Die Lösungen im System II sind identisch mit den Lé 
sungen der Rayleighschen Gleichung (8) und genügen de 


(w — 


oder (vgl. die Bemerkung zu (14a) $. 592) ganz allgemein 


y 
R d 
y,=(w— 0) =0 fir y=+1 und y=-1. 


Die letztere Form unterscheidet sich von der ersteren in ge 
wissen Ausnahmefällen, die nachher besprochen werden sollen; 
außerdem stellt (28) natürlich nur eine erste Näherung (« =0) 
dar. Für die Lösungen von (28) müssen jetzt vier verschieden 
Möglichkeiten unterschieden werden: Entweder 1. besitzt (28) 
Lösungen mit komplexem c. Dann ist das Profil immer labi, 
da stets auch der konjugiert komplexe Wert zu c eine Lösung 
darstellt; oder 2. es gibt Lösungen von (28) mit reellem & 
Dann bezeichnen wir mit Prandtl a. a. O. das Profil ab 
„schwingungsfähig“‘. Dies kann nach (21) nur vorkommen, 
wenn an der Stelle w=c) w’=0 ist, wenn also das 
Profil entweder einen Wendepunkt besitzt oder aus lineare 
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Stücken zusammengesetzt ist. Oder 3. es gibt reelle Werte 
yon ¢; die wenigstens den reellen Teil von 


m Null machen. Oder endlich 4. Es tritt keiner dieser drei 
Fälle ein, (28) besitzt keine Lösung, Im Fall (8) und (4) be- 
gichnen wir das Profil als „nicht schwingungsfähig‘‘. Wir be- 
haupten nun, daß Fall 1. stets Labilitdét, Fall 3. und 4. stets 
Stabilität, Fall 2. im allgemeinen Labslität des zugrunde ge- 
legten Profils ergibt. Für 1. und 4. haben wir dies im vorher- 
gehenden schon bewiesen. Im Fall 3. setzen wir ¢ = c¢, + ic,, 
woc, denjenigen reellen Wert von c bedeutet, für den der 


reelle Teil von +1 
f 
(w — e)* 
=1 


verschwindet. Dann wissen wir aus § 3, daß für c,=0 der 


aber gleich — a ai. Also muß es aus Stetigkeitsgriinden 


(vgl. § 3) eine Stelle c,>0 geben, an welcher auch der 
imaginäre Teil des Integrals (28) verschwindet. Die vier so 
charakterisierten Lösungen von (28) liefern also eine Größe c 
mit positiv imaginärem Bestandteil, also stabile Schwingungen. 

Der Fall 2. schließlich erfordert etwas ausfihrlichere Rech- 
nungen. Bevor wir sie ausführen, bemerken wir, daß zu 2. auch 
zwei Lösungstypen von (28) gehören, die sich nicht in der Form 


darstellen lassen. Wenn nämlich w(+1)=w(-—1) ist, so 
lautet eine Lösung von (28) g=w-—w(+1); in der Tat 
ist hier gp = 0 fir y=-+1 und y=—1. Ferner kann es 
vorkommen, daß z.B. w für w=w(+1) unendlich wird. 
Dann ist auch 


y 
d 
= (w—w(+ 1) f 
—L 


tine Lösung von (8), die den Grenzbedingungen genügt. Wir 
wollen hierauf nicht weiter eingehen, da dieser Fall im Teil II 
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ausführlicher diskutiert werden wird. Wesentlich ist abe 

noch zu bemerken, daß der Unterschied zwischen Fall 2. und $, 

sehr groß ist und daß es z.B. in keiner Weise genügt, nach 

Rayleigh ein krummes Profil durch ein Polygon anzunäherm, 
Zur Veranschaulichung dieses a ist 


als a ee von c qualitativ in Fig. 1 dargestellt, wo die 
ausgezogene Kurve dem krummen, die gestrichelten Kurva 
dem aus linearen. Stücken bestehen 
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Re (J(e)) 


— fil 
----- Näherungspolygon den Profile entsprechen. Man sieht, 
iz daß jeder Knick eine neue Wurad 
Fie (J) =O verursacht, weil J im 
Punkte c = Wp beim geknickten 
" Profil sich wie „— verhält. Die 


entspricht dem bekannten Rayleigh- 
schen Satz, daß es soviel Schwin 
gungswurzeln wie Knicke gibt. Trotz 
dem besitzt das krumme Profil kein 
Schwingungswurzel. Nach diese 


| 


x = Stellen, an denen 


so erhalten wir nach leichten Umformungen: 


© = “Knick Bemerkung kehren wir zu unserer Be 
hauptung zurück, daß die schwingungs 
fähigen Profile im allgemeinen bd 
IR Berücksichtigung der Reibung lalil 


- 00 


+1 


-00 


Fig. 1. 


werden. 

Zum Beweise dieser Labilität 
kehren wir zur Gleichung (25) un 
zu den genaueren Lösungen im Sp 
stem II zurück. Da wir wissen, dab 
c bis auf Größen der Ordnung (« = , 
reell ist, können wir 


\ 


annehmen. Vernachlässigen wir ferne 
die Glieder der Ordnung (« R)-’* in 1 
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Setzen wir weiter c = cy + 6, wo Cy die Nullstelle von J, 6 eine 
kleine Größe der Ordnung (x R)-"* bedeutet. « nehmen wir 
der Einfachheit halber positiv an, dann können wir auf der 
rechten Seite c durch c, ersetzen und die linke in eine Tay- 
lorsche Reihe nach 6 entwickeln. So ergibt sich, wenn wir 
die Taylorsche Reihe mit dem 2. Gliede abbrechen, was wir 
als hinreichende Näherung voraussetzen: 

J(c) = Ic) + 


und aus (25) wegen 


0; 9 =V —iaR(w—e) 
(über das Vorzeichen vgl. 8. 585). ERBEN 
dJ zu 1-i 1+i 
-V2akeu= 
de — w(— + [w (+ 1) — 


Hieraus folgt wegen cg — w(— 1) > w(+ 1) — «, (cy soll näher 
an w (-+ 1) liegen!), daß der imaginäre Teil von 4, also auch 


der von c und von ß dasselbe Vorzeichen hat, wie 7, REN 


- aJ 
und daß Schwingungen, welche einem negativen Werte von 7 
entsprechen, labilen Charakter tragen. Wenn also unser teil- 


weise lineares Profil noch die Eigenschaft hat, dap 47 <0 


ist an der Stelle w = c, so ist es labil. Diese Bedingung sf <0 


ist aber sehr häufig, z. B. immer dann erfüllt, wenn der Punkt 
w=c in der Nähe der einen Wand (z.B. y = + 1) liegt und 
das Profil von der Stelle w=c bis zum Rand linear verläuft. 

Zusammenfassend folgern wir: 

Die Labilität bzw. Stabilität eines Profils kann bei allen 
von uns bisher betrachteten Profilen schon durch das Verhalten 
derselben bei reibungsloser Flüssigkeit entschieden werden. Profile, 
die im letzteren Falle ungedämpfter Schwingungen fähig sind, 
werden bei Berücksichtigung der Reibung unter bestimmten 
Voraussetzungen labil. Die letzteren Profile müssen, wie oben 
gezeigt, ganz spezielle Eigenschaften haben, z. B. entweder 
teilweise aus linearen Stücken zusammengesetzt sein, oder sie 
müssen einen Wendepunkt w’’ = 0 besitzen (vgl. oben). 

Diese Profile der Gattung 2. sind aber zugleich die einzigen, 
die noch ein physikalisches Interesse beanspruchen können. Denn 
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nor sie besitzen ein Verhalten hinsichtlich ihrer Stabilität, 
das etwa den Vermutungen von Reynolds entspricht. Wir 
werden im folgenden zeigen, daß diese Profile wirklich im 
allgemeinen eine „kritische“ Reynoldssche Zahl besitzen (mit 
Ausnahme der geknickten Profile). 


86. Die Reynoldssche Zahl der Stabilitätsgrenze; 
numerische Rechnung am Parabelprofil. 


Wenn also ein Profil gegeben ist, das. bei reibungsloser 
Flüssigkeit ungedämpfte Schwingungen zuläßt und mit Reibung 
labil ist, so taucht die Frage auf: Welches ist der kleinste Wert 
der Reynoldsschen Zahl, für den Labilität eintritt? Um diese 
Frage zu beantworten, reichen die vereinfachten Gleichungen 
(25), (26) usw. nicht aus. Wir müssen zur Gleichung (23) und 
den Formen (11) und (14) für die Integrale 9, g. 9, 
zurückkehren. Es ist aber natürlich ganz unmöglich, allgemein 
für ein beliebiges Profil w die kritische Reynoldssche Zahl 
als Funktion von w und Integralen über w darzustellen; unsere 
Aufgabe soll es nur sein, den Weg anzugeben, auf dem man zur 
kritischen Geschwindigkeit gelangt und dann an einem spe 
ziellen Beispiele die Rechnung durchzuführen. 


a 

A 


Da in unseren letzten Rechnungen « und R nur in der 
Kombination « R vorgekommen sind (weil wir «? als klein an- 
nahmen), können diese Rechnungen auch bestenfalls nur einen 
kritischen Wert für « R, nicht aber für R allein liefern. Wir 
müssen daher zunächst das Verhalten der Wurzeln von (%) 
bei wachsendem «? untersuchen. Von den Wurzeln von (2) 
interessieren uns hierbei nur die im Lösungssystem II, welche 
der Gleichung 9, (+ 1) = 0 genügen. 


An Stelle von (28) müssen wir also die Gleichung 


( 
g =0 fir y=—1, y=+1 
diskutieren. Besteht das Profil. wie bei Rayleigh, aus 


linearen Stücken, so gibt es (vgl. S. 598) stets eine Wurzel von 
(28a) für jeden Knick und diese Wurzeln bleiben für jeden 
Wert von «? bestehen. Das geknickte Profil liefert also heinen 
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Maximalwert von «? und kann daher auch niemals auf eine 
kritische Reynoldssche Zahl führen.!) 

Anders ist dies, wenn (vgl. S. 597) eine Lösung von (28) 
bzw. (28a) mit reellem c deswegen möglich ist, weil entweder 
irgendwo im Profil w’’ =0 ist oder weil w (+ 1) = w(— 1), 
g=w-— w(+ 1) eine Lösung von (28) darstellt. Diese letz- 
teren Lösungstypen liefern nämlich immer nur für einen ganz 
bestimmten Wert von «* eine Lösung von (28a). Für w’= 0 
ist ¢ schon festgelegt dadurch, daß eben w =c sein soll für 
w’'=0; also definiert die Gleichung (28a) einen ganz be- 
stimmten Wert «®. Für den Fail w (+ 1) = w (— 1) aber gibt 
es offenbar nur für «® = 0 eine Lösung von (28a). 

Für diese Art von Lösungen der Gleichung (28) bzw. (28a) 
die im Limes R = oo durch einen ganz bestimmten Wert von 
«charakterisiert sind. werden wir erwarten, daß sich bei Berück- 
sichtigung der Reibung « auch nur um geringe Beträge von seinem 
„bestimmten‘‘ Wert entfernen kann, das Auftreten eines Maxi- 
malwertes (und im Falle w” = 0 auch eines Minimalwertes) für « 
ist bei diesen Profilen sehr verständlich. Alle Schwingungen also, 
deren Wellenlänge kleiner ist, als eine bestimmte kritische Wellen- 
länge, sind in solchen Fällen für alle Werte von « R gedämpft. 

Nachdem eine obere Grenze für «? gefunden ist, wird man 
die ungefähre Größe für die untere Grenze von « R zu be- 
stimmen suchen. Eine einfache Untersuchung der Glei- 
chung (25) lehrt, daß wesentliche Veränderungen im imagi- 
nären Teil von c erst eintreten, wenn der Exponent von e in 
der Naherungsdarstellung (11) in 9, (+ 1), 9, (+ 1) auf Werte der 
Größenordnung 1 gesunken ist. Ist dies aber der Fall, so 
haben wir den kritischen Wert, bei dem der imaginäre Teil 
von ¢ von negativen zu positiven Werten übergeht, sehr bald 
erreicht, wie auch das Zahlenbeispiel zeigen wird. Nehmen 
wir an, daß w zwischen w = c, (cy = reeller Teil von c) und 
w(+1) wesentlich linear verläuft, so lautet daher die Be- 
dingung für die ungefähre Größe von « R: 

(29) (a [0 (+ 1) — 
w’ (+ 1) 

1) Hierbei ist immer noch vorausgesetzt, daB R und « R groB und 
t€R sei. Möglicherweise treten also kritische Reynoldssche Zahlen 
dann auf, wenn diese Voraussetzungen nicht mehr gelten. Die betreffen- 
den Reynoldsschen Zahlen R dürften aber dann bei so kleinen Werten 
liegen, daß sie physikalisch sicher keine Bedeutung haben. 


u 
tät, 
Wi 
it 
in 
mit 
oser 
ung 
ert 
iese 
gen 
und 
rein 
zahl 
Bere 
zur 
pe 
er 
an- 
a 
(2 ) er 
iche 
| q 
1 
+ i 
von 
den 
| ES 
nen 


602 

N 2 

(29a) 


Da w (+ 1) — cy in den uns interessierenden Fällen klein. sein 
wird, so kann man aus (29) vermutungsweise auf große kritische 
Reynolds sche Zahlen schließen. Zugleich bemerken wir hier, 
daß es für einen bestimmten Wert von R stets nicht nur einen 
Maximalwert, sondern auch einen Minimalwert für « de 
labilen Schwingungen geben wird. Dies folgt daraus, daß wir 
ja einen Minimalwert von « R (nicht R) gefunden haben. 

Als numerisches Beispiel für unsere bisherigen allgemeinen 
Rechnungen wählen wir das Parabelprofil, weil es physikalisch 
am interessantesten ist. Es gehört zur Klasse der „schwingungs 
fähigen Profile“ vom Typus w (+ 1) = w(— 1). 

Wir wollen auch hier nur die ebene Bewegung betrachten, 
d.h. nicht die Poiseuillesche Strömung in Rohren, sondem 
die Strömung, die zwischen zwei parallelen ruhenden Wände 


(y=+1, y=-—1) unter Einfluß eines konstanten Druck 
gradienten herrscht. Wir setzen also 
(30) v=1-— 


Aus der Symmetrie von w und w — c läßt sich ableiten, 
daß 9 eine gerade Funktion von y sein muß.!) Von den Lö 
sungen der Gleichung (7a) greifen wir also zwei symmetrische 
partikuläre Integrale heraus und suchen die Grenzbedingungen 
an der einen Wand, etwa y = — 1, zu befriedigen. Die an der 
anderen Wand sind dann von selbst erfüllt. Als das eine dieser 
symmetrischen Integrale kann offenbar einfach 9, genommen 
werden. Für das andere wählen wir: 

9 (y) 
(0) (0) 

1) Zerlegt man 9 in einen in y geraden und einen ungeraden Teil, 
so muß, wegen der Symmetrie von w — c und w, jeder Teil von @ für sich 
der Differentialgleichung (7a) und den Grenzbedingungen genügen. Fir 
die allgemeine Stabilitätsuntersuchung des Profils 1 — y? ist es dahe 
ausreichend, die beiden Fälle „‚p gerade‘ und ,, ungerade‘ getrennt und 
nur diese beiden Fälle zu behandeln. Es ist aber leicht zu sehen, daß die 
Annahme symmetrischer Schwingungen, d.h. „p ungerade‘ zu keine 
Lösung von (28) und daher nicht zu labilen Schwingungen führt. Di 
Annahme ,,g gerade“ genügt also zur Stabilitätsuntersuchung. Die 
ist insofern bemerkenswert, als hiernach alle symmetrischen Schwit 
gungen stabil und nur unsymmetrische Störungen labil sind. 
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Wie aus (29a) hervorgeht, wird c bei unserem Profil in 
Nähe der kritischen Geschwindigkeit klein sein von der Ord- 
nung (x R)="*, wir werden also in den folgenden Rechnungen 
Glieder von höherer als erster Ordnung in c streichen. Ferner 
stellen wir fest, daß 9, (0) > ¢, (0) sein wird, so daß sich in 
der Umgebung von w = 0 und w = c die zweite symmetrische 
Funktion p einfach auf 9, (y) reduziert. Aus (16) geht dann 
hervor, daß wir zur Befriedigung der Grenzbedingungen für 
y=-— 1 die beiden Integrale 9, und gy, — 4g, zur Verfügung 
haben. Gleichung (23) geht damit über in 


9 (— 1) 9 (— 1) — 4+, (—1) 
(— 1) #1 (— 1) (— 1) ; 
In 9, ist hierbei als untere Grenze für die auftretenden Inte- 
grale immer y = 0 zu nehmen, um die Symmetrie von 9, zu 
gewährleisten. In 9, werden wir ferner nur bis zu Größen der 
Ordnung «* entwickeln und bei der Entwicklung nach (« R)-! 
mit den Gliedern der Ordnung (« R)-! abbrechen. Wir schreiben 
nun (31) in der Form 


(31) 


gi’ (— 1) — 1) _ 
Setzt man (11) und (14) ein, so ergibt ich 00° 


-1 


2fV-iaRw— ejay 
Yo 
9 
+ *ViaRe=— 
2fV-siaRw- gay 
uh z dy(w + 
aR 
dy a 
0 0 


Da c sehr klein wird, nehmen wir in erster Näherung w von o 
bis ¢ als linear an: w~2(y+1). Dann wird ~~ 
| 


c-wdy=m—-c. | 
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-(2@R)", so entsteht aus (83) 


-1 


P 
Setzen wir z= 3° 


Ä 


Diese Gleichung ist zu (26) völlig analog. Uns interessiert vor 
allem der Grenzwert R bzw. z, bei dem die labilen Schwin- 
gungen sich in stabile verwandeln, also der imaginäre Teil 
von c gerade Null ist. Dieser Grenzwert z wird natürlich noch 
von « abhängen. Wir nehmen also jetzt c als reell an und er. 
halten dadurch den Grenzwert von z bzw. R als Funktion 
von «: Den Minimalwert von R auf dieser R («)-Kurve wird 
man als charakteristische Reynoldssche Zahl für das Parabel- 


EEE profil ansprechen. Die Durch- 
ae es rechnung zeigt nun, daß man 
Stabilität durch die Form (35) der Stabili- 
tätsgleichung nur den oberen 

200-: . 

‘Teil der Kurve R= 'R (a) (in 

or; Labilitét der Fig. 2 ausgezogen) erhält — 
120- \ was ja nach den Betrachtungen 

\ 

Mr 3 von $ 6 zu erwarten war, daß 

f sr ines oe man aber den unteren Teil der 
| et, _KurveR=R (x) nur berechnen 
ao UY @ @ % % kann, wenn man für 9,, % 
Fig. 2. andere Näherungen [vgl. (19a)] 


als die asymptotischen For- 
meln (11) verwendet. Die kritische Reynoldssche Zahl be 
zeichnet eben gerade das Gebiet, in welchem die asym- 
ptotischen Formeln ihre Gültigkeit verlieren. Da dieser Um- 
stand zu sehr komplizierten numerischen Rechnungen führen 
würde und wir (vgl. $ 7) allgemein der hier charakteri- 
sierten Art von Labilität keine wesentliche physikalische Be 
deutung beimessen können, haben wir zur Berechnung de 
unteren Teiles (in Fig. 1 gestrichelt) der Kurve R= R(@) 
rohe Schätzungen benutzt, die keinen Anspruch darauf machen 
können, quantitative Ergebnisse zu liefern. Das qualitative 
Verhalten der Kurve wird aber wohl richtig wiedergegeben. 
Wir entnehmen also aus der Fig. 2: BE 
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1. Es gibt sowohl einen Maximalwert von «, wie einen 
Minimalwert von R, nach deren Über- bzw. Unterschreitung 
die Labilität aufhört. 

9. Für einen bestimmten Wert von R existiert sowohl ein 
Mazimal- wie auch ein Minimalwert von «; innerhalb dieser 
Werte herrscht Labilität, außerhalb Stabilität. 

8. Der Maximalwert von « liegt etwa bei « = 0,7 (a? = 4). 
Der Minimalwert von R bei Größen der Ordnung 10%. Eine 
einigermaßen genaue Berechnung dieses Minimalwertes aus der 
Figur ist nicht möglich. 


§ 7. Physikalische Diskussion der Resultate des Teiles I. 


Fassen wir noeh einmal ausführlich alle Ergebnisse zu- 
sammen, die wir für das Stabilitätsproblem gefunden haben. 
In erster Linie hat sich im Verlaufe der Rechnung gezeigt, 
daß man allgemein die Frage nach der Stabilität eines Profils 
bei reibender Flüssigkeit entscheiden kann, wenn man, wie 
Lord Rayleigh, den Grenzfall reibungsloser Flüssigkeit 
(Gleichung 8) behandelt. Profile, die dann, d.h. für R = co 
labil sind, bleiben dies auch für hinreichend große endliche 
Werte von R ($ 4) — wie dies wohl von vornherein zu er- 
warten war. Ebenso erweisen sich Profile, die im reibungs- 
losen Falle keiner Schwingungen fähig sind, als stabil ($ 4) und 
Profile, die bei der Untersuchung durch Gleichung (8) ungedämpfte 
Schwingungen zulassen, im allgemeinen als labil ($ 5). Dieser 
letztere Fall ist offenbar der einzige, der physikalisch gegenüber 
der reibungslosen Hydrodynamik etwas Neues bedeutet. Wir 
müssen aber hervorheben, daß er, entgegen dem, was man aus 
Rayleighs Arbeiten zunächst schließen möchte, einen Aus- 
nahmefall darstellt. Wenn man von den Möglichkeiten w (+ 1) 
=w(—1), w’pana= 0 ($5) absieht, so ist es für das Ein- 
treten dieses Ausnahmefalls notwendige Bedingung, daß irgendwo 
in der Flüssigkeit w’’ = 0 ist. Die von Rayleigh eingeführten 
geknickten, aus linearen Stücken bestehenden Profile gehören 
zu diesen Ausnahmeprofilen; hieraus ist aber nur zu schließen, 
daß eben gekrümmte Profile zum Zwecke der Stabilitätsunter- 
suchung nicht nach Rayleigh durch Polygone angenähert 
werden dürfen (S. 598). Auch für überall gekrümmte Profile 
(w” + 0) findet man allerdings (§ 5) bei Benützung der Diffe- 
rentialgleichungen mit Reibungsgliedern eventuell Schwin- 
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gungen, die für jeden Wert von « R gedämpft sind, bei denen 
aber im Limes R =o der Betrag der Dämpfung wie (« R)=" 
gegen Null geht; man hat also auch hier fir R=w m. 
gedämpfte Schwingungen. Diese Schwingungswurzeln gehen 
aber verloren, wenn man die vereinfachte Differentialgleichung($) 
ohne Reibungsglieder zugrunde legt. Insofern liegt also auch 
hier keine Ausnahme der Regel vor, nach welcher dann, 
wenn die reibungslose Gleichung (8) ungedämpfte Schwingungen 
zuläßt, die Berücksichtigung der Reibung eine Anfachung er. 
gibt. Die Möglichkeit ungedämpfter Schwingungen für die 
Gleichung (8) muß aber, wie oben gesagt, als Ausnahmefall 
betrachtet werden. Das Parabelprofil gehört zu diesen Aus 
nahmeprofilen ($ 6). Fragen wir bei den „labilen Profilen“ 
weiter nach dem Wertebereich für R, innerhalb dessen Labilität 
stattfindet, so zeigt sich, daß auf „kritische Reynoldssch 
Zahlen“ auch im allgemeinen nur Profile der letzten Klasse 
führen können, d. h. nur solche, die ohne Reibung ungedämpfk 
Schwingungen zulassen. Von diesen lieferten aber die Rayleigh- 
schen ,,geknickten Profile‘ dennoch keine (vgl. auch $. 601, 
Anm. 1) kritische Reynoldssche Zahl. Bei diesen gibt es näm- 
lich wohl einen Minimalwert von « R, nicht aber eine 
Maximalwert von «, also auch keinen Minimalwert von R für 
die Labilität (§ 6). Nur diejenigen Profile der letzen Klasse, 
bei denen im reibungslosen Falle nur ein bestimmter Wert von a 
auf ungedämpfte Schwingungen führt (z. B. die Typen w” =0 
für w= ec), w (+1) =w(-— 1), ergeben einen Mazimalwert von 
a und einen Minimalwert von « R, also auch einen Minimalwer 
für R. Für einen bestimmten Wert von R gibt es daher bé 
diesen Profilen sowohl einen Mazimal- wie einen Minimalwer 
von « für die labilen Schwingungen. Alle diese Resultate sind 
in Übereinstimmung mit den bisherigen Stabilitätsunter- 
suchungen hydrodynamischer Profile.!) 

Es fragt sich jetzt, wie si’ h experimentell diese mathe 
matischen Resultate äußern werden. Es erscheint verwunder 
lich, daß sich empirisch die stabilen (z. B. das Couettesche‘) 


1) Vgl. die in der Einleitung zitierten Arbeiten. 

2) Vgl. jedoch hierzu die interessanten Stabilitätsuntersuchungen der 
Couetteschen Bewegung gegen dreidimensionale Störungen, G. J 
Taylor, Stability of, a viscous liquid contained between two rotating 
cylinders. Phil. Transact. of the Royal Society London 223, S. 289-3 
1922. 
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und die labilen (z. B. das Poiseuillesche) Profile ganz 
gleich verhalten. Oberhalb einer bestimmten Reynolds- 
schen Zahl tritt bei hinreichenden Störungen Turbulenz 
dn; macht man die Störungen möglichst klein, so läßt 
sch das Laminarprofil bei beliebig hohen Reynoldsschen 
Zahlen erhalten. Besonders die letztere Tatsache, die von 
Ekman (a. a. O.) am Parabelprofil geprüft worden ist, scheint 
bei den labilen Profilen im Widerspruch zur Theorie zu stehen. 
Es-ist aber leicht zu erkennen, daß dieser Widerspruch nur 
scheinbar ist.!) Je kleiner nämlich die äußeren Störungen sind, 
um so länger dauert es (besonders bei hohen Reynoldsschen 
Zahlen, da die Anfachung von der Ordnung (« R) =" ist (vgl. 
$5), bis sie die Bewegung merklich beeinflussen; bei der 
$tröomung in Rohren wird man also diesen Zeitpunkt immer 
soweit hinausschieben können, daß die betreffende Flüssigkeits- 
menge, um deren Stabilität es sich handelt, das Rohr schon 
wieder verlassen hat, wenn ihre Labilität in Erscheinung tritt. 
Die von uns berechnete Reynoldssche Zahl könnte also nur 
an in sich geschlossenen Rohren geprüft werden, bei denen 
immer dieselbe Flüssigkeitsmenge strömt.?) Die Versuche von 
Schiller (a. a. O.) andererseits, durch die gezeigt wird, daß 
# unterhalb einer bestimmten Reynoldsschen Zahl nur 
laminarbewegung gibt, können durch Stabilitätsuntersuchungen 
überhaupt nicht erfaßt werden, da man hier gar nicht als ur- 
sprüngliche Bewegung die laminare hat, sondern da es sich 
nur um die Existenz oder Nichtexistenz einer turbulenten 
Bewegungsform dabei handelt. Jedenfalls ist aus allen diesen 
Betrachtungen nur zu schließen, daß sich durch Stabilitäts- 
betrachtungen allein das Turbulenzproblem eben durchaus 
nicht lösen läßt. 

Immerhin können uns die vorhergehenden Untersuchungen 
auch für das eigentliche Ziel, die Berechnung der turbulenten 


1) Auf diese Möglichkeit, die Ekmanschen Versuche als eine Art 
Anlaufeffekt aufzufassen, hat mich Hr. Professor Prandtl hingewiesen, 
dem ich hierfür und für viele andere wertvolle Hinweise an dieser Stelle 
den herzlichsten Dank aussprechen möchte. 

2) Vielleicht sind aber die Störungen der Stabilität, die Ruckes 
& a. 0. bei ziemlich kleinen Reynoldsschen Zahlen beobachtet, durch 
die Labilität nach $ 7 hervorgerufen. Dies wäre dann sehr wohl denkbar, 
wenn die kritische Reynoldssche Zahl nach $7 unter derjenigen liegt, 
bei welcher zum ersten Male Turbulenz (vgl. Teil II) möglich ist. 
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Bewegung, wesentliche qualitative Ergebnisse liefern. Wenn 
wir nämlich die turbulente Bewegung auffassen als eine be 
stimmte Grundstrémung mit überlagerten ungedämpften Schwin- 
gungen, so können wir aus unseren Rechnungen schließen, 
daß der Minimalwert R, für den dieser Bewegungstypus mig. 
lich ist, wahrscheinlich ebenfalls bei Werten der Größen 
ordnung 103% liegt; daß die Wellenlänge der ungedämpften 
Schwingungen bei 22 h/2, nämlich « bei Werten der Ordnung] 
liegt und daß « für ein gegebenes R zwischen bestimmte 
eventuell sehr enge Grenzen eingeschlossen ist; daß ferne 
diese Schwingungen, wie aus der Kleinheit von w (+ 1) -t 
zu schließen ist, den Charakter einer Wandstérung haben, 
Diese qualitativen Ergebnisse sind von der speziellen Form 
der Grundströmung ganz unabhängig. Uber solche qualitativen 
Anhaltspunkte hinaus aber liefern die bisherigen Rechnunge 
nichts, was uns der eigentlichen Lösung des Turbulenzproblems 
näherbrächte. 

a II. Teil: Die turbulente Bewegung. 
” § 1. Formulierung des mathematischen Problems. 
led Diejenige Reynoldssche Zahl, die man für gewöhnlich ak 
die „kritische‘‘ bezeichnet, die z.B. in den Versuchen vo 
Schiller gemessen wird und die den Eintritt der Turbulem 
bei hinreichend großen Störungen angibt, hat mit Stabilitäte 
fragen und mit der Laminarströmung nichts zu tun; sie ist 
durchaus eine charakteristische Konstante der turbulenten Be 
wegung. Ebenso lassen das Blasiussche Widerstandsgesets 
und die aus ihm abgeleitete bekannte Folgerung, daß de 
turbulente Geschwindigkeit in der Nähe einer Wand mit de 
1/, ten Potenz der Wandentfernung anwächst, deutlich erkennen, 
daß die sogenannte „turbulente‘‘ Bewegung ihre eigenen woll 
definierten Gesetzmäßigkeiten hat und daß sie eine zweile 
mögliche Bewegungsform der reibenden Flüssigkeiten darstellt 
Zur Lösung des Turbulenzproblems kann man also nur ve 
suchen die turbulente Bewegung ihrer Undefiniertheit zu ent 
kleiden und sie soweit zu idealisieren, bis sie der mathematischal 
Analyse durch die Stokesschen Gleichungen zugänglich wird. 

Das „Turbulenzproblem‘“ der Hydrodynamik ist @ 
Problem der energetischen, nicht der dynamischen Stabilität 
Es gibt zwei verschiedene Bewegungsformen der reibendél 
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Flüssigkeit, von denen jede einen bestimmten Wertebereich der 
Reynoldsschen Zahl besitzt, innerhalb dessen sie möglich 
ist. Die Laminarströmung ist möglich von R = 0 bis R=m, 
wird aber unter Umständen oberhalb eines bestimmten Wertes 
von R dynamisch labil. Die turbulente Bewegung dagegen 
existiert erst oberhalb eines bestimmten kritischen Wertes von 
R, ist dann aber immer!) energetisch stabiler als die Laminar- 
bewegung. In demjenigen Bereich von R, in welchem beide 
Bewegungsformen möglich sind, kann man daher stets die 
Flüssigkeit aus dem laminaren Zustand durch hinreichend 
große Störungen, in den turbulenten Zustand ,,herunterfallen“ 
lassen. 

Um die turbulente Bewegung näherungsweise mathe- 
matisch behandeln zu können, fassen wir wieder die Strömung 
gwischen zwei parallelen Wänden ins Auge und machen zu- 
nächst folgende Annahmen?): Die Strömung soll a) um die 
X-Achse symmetrisch, die begrenzenden Wände in Ruhe sein; 
sie soll b) periodisch in der X-Richtung mit der Periode 
9x/a*) und c) periodisch in der Zeit mit der Periode 2/ß sein. 
d) Alle Störungen sollen sich mit der Geschwindigkeit ß/« 
relativ zur X-Achse fortpflanzen, d. h. wenn wir die Bewegung 
ineine Fourierreihe entwickeln, so sollen in den Exponenten 
von e nur Vielfache von i (ft — ax) vorkommen.*) 

Die Fourierentwicklung der Stromfunktion soll also lauten: 


(36)) = Poly) + [y) + gq, 
Die mathematische Aufgabe besteht dann in der Er- 
mittelung der (nach a) ungeraden Funktionen 9, 9, 9--- 


1) Vgl. F. Noether, a. a. O. 

2) Auch dieser Ansatz, der eine einfache Verallgemeinerung des 
Sommerfeldschen Stabilitätsansatzes darstellt, wurde zur Untersuchung 
der turbulenten Bewegung selbst von F. Noether: Zur Theorie der 
Turbulenz, Jahresber. d. deutsch. Math. Vereins 28. S. 138. 1914 ohne 
weitergehende Folgerungen angegeben. 

3) Zur Annahme eines bestimmten « berechtigt uns das Ergebnis 
von Teil I, daß « zwischen bestimmte Grenzen eingeschlossen ist, die 
insbesondere in der Nähe des Minimalwertes von R sehr eng sind, 

4) Wir brauchen nicht hervorzuheben, daß die wirklichen Be- 
wegungsvorgänge zweifellos viel komplizierter sind; trotzdem kann man 


wohl erwarten, daß diese Ansätze qualitative Aussagen über die Tur- 
bulenz gestatten. 
Annale 
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(Gi, konjugiert zu ,, 9). Der Grad der Konvergenz 
der Reihe (36) ist zunächst völlig unbekannt, die Konvergenz. 
frage kann erst entschieden werden nach der Berechn » 
von 9%, 9, --- Wollen wir die Genauigkeit bis zur n™ 
Näherung treiben, d. h. wollen wir 9,...9, berechnen, so be 
kommen wir offenbar (n + 1) simultane Differentialgleichungen 
für die (n + 1) unbekannten Funktionen g,...@,- Im folgen- 
den werden wir teils der ersten, teils der zweiten Näherung be 
dürfen. Wir gehen also mit dem Ansatz (36) in Gleichung (4% 
ein, vergleichen beiderseits die Koeffizienten der periodischen 
Funktionen und brechen mit dem Gliede e?'#!-«) und daher Win 
mit 9, ¢2 ab. Für 9, schreiben wir (wie in (6) für ©) w Aor 
So entstehen drei simultane Differentialgleichungen (die gleich ih 
zeitig entstehenden konjugierten Gleichungen brauchen wir 
nicht hinzuschreiben): 


aye | Pi Pi — Fr Pr +22 = hal 


+ 2H; — 


we 


” ” 1 


= + 160g), 


Von diesen Gleichungen läßt sich die erste zweimal integrieren 
und liefert: 


(37a) — + 22 Fs — Fa G2) = —C-Gj; 
C und C, bedeuten willkürliche Integrationskonstanten. ' 

e 
Da die linke Seite von (37a) und w nach Forderung (a) ain 
ungerade in y sind, so muß in unserem Falle C verschwinden} 


Gehen wir von der zweiten auf die erste Näherung zurich} ‚$t 
so reduziert sich unser Gleichungssystem auf zwei simultane itr’ 
Differentialgleichungen für w und 9,: aad ee lag 
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rgenz- 
— p,)(w — 2) — = — Brg,” 

— 
lah linschaltungsweise wollen wir uns jetzt überlegen, was an 


ng be. Helle von (38) tritt, wenn wir nicht eine um die X-Achse 
ng (4) etrische Strémung (Forderung a), d.h. die Strémung 
ischen it Flüssigkeit unter Druckgefälle zwischen zwei ruhenden 
dahe|Wamden, sondern eine um die X-Achse antisymmetrische . 
frömung (d.h. Strömung zwischen zwei relativ zueinander 
Dewegten Wänden ohne Druckgefälle wie beim Couetteschen 
fall) betrachten. Forderung (b) und (c) sollen bestehen bleiben. 
Der Ansatz (36) wird dann nicht mehr genügen, da 9, 
, usw. bei beliebigem ß/«& keine geraden Funktionen mehr 
ind; wir müssen, um das ganze Strömungsbild ungerade zu 


” + Ighalten, auch noch die symmetrischen Schwingungen der 
form e‘(-#t-.a2) in den Ansatz für y aufnehmen d. h. w muß 
beginnen mit den Gliedern 

alg) ht Ft + (— y eit tax) 4 p,(y)e 


y,), [Das hat zur Folge, daß in y schließlich alle Glieder der Form 
ge#t-naz) auftreten (Wegfallen der Forderung (d)). 


7 An Stelle von (38) ergibt sich: ar 

a 

C, - 

+ a*g,). 


Die beiden Gleichungen des Systems (38) bzw. (39) haben eine 
ng (länfache anschauliche Bedeutung. 


Die zweite Gleichung ist nichts anderes, als unsere frühere 
arück| Stabilitatsgleichung‘ (7), welche die Amplitude der einer Grund- 
ultane] #trömung w überlagerten Schwingung bestimmt und die Grund- 
lage für unsere Untersnchungen im Teil I bildete. Die erste 
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Gleichung aber stellt den Impulssatz dar. Die linke Seite 
dieser Gleichung gibt nämlich im wesentlichen den durch dig 
turbulente Wirbelung im Mittel übertragenen Impuls an!), da 
Glied mit w’ rechts stellt die laminare Impulsübertragung dar, 
die Konstante C bzw. C, y ist die ,,Konstante des Impuk- 
satzes. 

An den Wänden ist wegen der Grenzbedingungen 9, =,’ =0, 
Also ist dort w’=C bzw. = (, ywana; an den Wänden über 
wiegt also der laminare Impulstransport den turbulenten, 
W’wana wird im allgemeinen sehr groß sein (vgl. den nächsten 
Paragraphen). In der Kanalmitte aber, d. h. im ganzen Kanal 
außerhalb der unmittelbaren Wandnähe, ist w’ von der Größen- 
ordnung 1, also sehr klein gegen C bzw. C, y. Hier wird daher 
der turbulente Impulstransport den laminaren völlig über 
wiegen. 

Der Bauart der Systeme (38) und (39) entspricht es, daß 
wir eine triviale Lösung derselben sofort angeben können, 

7 = 9, 


nämlich 

d.h. wir kommen so zur Laminarbewegung zurück. 4 
Unser Ziel muß es jetzt sein, über die nichttrivialen 

_ Lösungen von (38) und (39) bestimmte Resultate zu erhalten, 


w=( bzw wW=(,y, 


8% Die turbulente Bewegung in Wandnähe und das 
Widerstandsgesetz. 

Das wichtigste Ergebnis über das Verhalten von w in un 
mittelbarer Nähe der Wände ist das von v. Kärmän (a. a. 0) 
durch Dimensionsbetrachtungen aus dem Blasiusschen Wider- 
standsgesetz abgeleitete Gesetz, daß w in der Nähe eine 
Wand wie 7'* anwächst (wenn 7 die Entfernung von der Wand 
darstellt). Wir wiederholen kurz den Kärmänschen Ge 
dankengang, da wir uns dadurch allgemein darüber orien- 


1) Wir meinen hier den mittleren Impuls in der X-Richtung, der 
bei unserem Problem in der Y-Richtung übertragen wird. Der Impuls 
in der X-Richtung ist im wesentlichen gleich u, die Geschwindigkeit des 
den Impuls tragenden Teilchens in der Y-Richtung ist v, also der in der 
Zeiteinheit übertragene Impuls u-v, im Mittel u-v, was im Falle (36) 
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tieren können, was wir — auch ohne das Blasiussche Gesetz 
zu kennen — für das Verhalten von w in Wandnähe zu er- 
warten haben. 

Die an einer Wand angreifende Schubspannung r (d.h. also 
der Widerstand) muß sich, wie aus Dimensionsbetrachtungen __ 
leicht zu ersehen ist, darstellen lassen in der Form: 11: 


U 


yo x eine bestimmte dimensionslose Konstante bedeutet. = = 


Wenn wir speziell ein Potenzgesetz annehmen, so ist on 


Aus (40) folgt umgekehrt: rsh 
E (ay 406 
Die Geschwindigkeitsverteilung in Wandnähe muß sich dann 
durch eine Gleichung der Form darstellen: 
1 


(w ist hier wieder [vgl. (1) und (6)] dimensionslos gewählt und 
enthält daher U im Nenner, 7 bedeutet den Abstand von der 
Wand). 

Nehmen wir wieder in 1. Näherung ein Potenzgesetz an 
(6 sei eine dimensionslose Konstante): 


1-8]1+é e 
Fordert man jetzt, daß die Geschwindigkeitsverteilung in un- 
mittelbarer Nähe der Wand nur von 7, u, oe, nicht aber von h 
abhängen soll, was physikalisch sehr plausibel ist, so folgt: 
(41) = é. 
Fir & = 3, wie es dem Blasiusschen Gesetze entspricht, 
ergibt sich also « = +. 
Um uns über die physikalische Bedeutung dieses Resul- 
tats klar zu werden, bemerken wir: w ~~" bedeutet, daß 
W“=dw/dn am Rande unendlich ist, daß also w sich un- 
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endlich an die Wand anschmiegt. Nun ist es aber klar, daß 
in Wirklichkeit w an der Wand nicht unendlich sein k 
denn w’ bedeutet ja gerade im wesentlichen die Schub. 
spannung an der Wand und ist daher nach (40) bis auf einen 
von R unabhängigen Zahlenfaktor gleich R* zu setzen, 
(42) w' pana ~ RS. 
w’ an der Wand ist also bei den uns interessierenden großen 
Werten von R sehr groß. Die Geschwindigkeitsverteilung 
w-n" wird ihrer Ableitung entsprechend nur im Grenz 
fall unendlich großen Wandabstandes oder reibungsloser 
Flüssigkeit (R= oo) in Strenge gültig sein. Noch einfacher 
läßt sich dieser Sachverhalt: überblicken, wenn wir das ge 
fundene Gesetz w n"' in der Form schreiben n = w’, Au 
der Tatsache endlicher Schubspannung wissen wir, daß das 
erste Glied der Potenzentwicklung n (w) von der Form y,u 
sein muß. Dieses Glied ist aber sehr klein, nämlich im wesent- 
lichen gleich dem reziproken Werte von w’, also von de 
Ordnung R~* (vgl. 42). 
Die Meinung des abgeleiteten Gesetzes w ~7'* ist also 
offenbar die, daß die Reihenentwicklung von 7 (w) anfangen 
soll mit den Gliedern 
wo y, außerordentlich klein ist, und daß en das PR ied 
y, w für einigermaßen große Werte von w gegen das zweile 
gestrichen werden kann. 
Nach den vorhergehenden Darlegungen erwarten wir um 
abhängig von der Gültigkeit des +-Gesetzes für das Aus 
sehen der Grundströmung der turbulenten Bewegung in der 
Mitte kleine Krümmungen, in Wandnähe Anschmiegen der 
Grundströmung an die Wände. 
Für ein solches Profil gelten die Untersuchungen de 
Teiles I nicht unmittelbar, denn wir haben dort w’, w” usw. 
als endlich vorausgesetzt. Doch lassen sich diese Unter 
suchungen leicht auf solche Profile wie das hier vorkommende 
verallgemeinern (vgl. $ 5, 8. 597). Insbesondere wird die Lö 
sung der reduzierten Gleichung (8) (also im R= oo) mit 
Befriedigung der Grenzbedingungen hier besonders einfach; 
das soeben charakterisierte Profil gehört nach § 5, 8. 597 zu 
den „schwingungsfähigen‘“, es existiert eine Lösung von (29) 
mit reellem ¢; dies ist außerordentlich wichtig, zeigt es doch; 
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daß die „turbulenten Profile“ immer „labil‘“ sind nach § 5, 
oder anders ausgedrückt, daß gerade die Abweichung (42) vom 
jaminaren Widerstandsgesetz die Möglichkeit eines „labilen“ 
Profils, also der Turbulenz gibt. Pell 2 

Die Lösung von (28) heißt für « = 0 rn = 


Daß g für y= — 1 Null wird, haben wir durch Wahl der 
unteren Grenze des Integrals erreicht; daß es aber auch für 
y=+ 1 Null wird, wegen der Wahl c=w(+1) ist leicht 
zu sehen aus folgender Umformung: 


w(-1) 
Das Integral der rechten Seite wird nämlich im Punkte 


w=w(+1) von niedrigerer Ordnung unendlich als Feet 
da w’ dort (im Limes R = oo) unendlich wird. Also ist y, = 0 
fir y= +1. Durch (44) haben wir im Limes reibungsloser 
Flüssigkeit die Amplitude der turbulenten Schwingungen dar- 
gestellt und aus den Grenzbedingungen den Wert für ß/a, 
nämlich $/& = w (+ 1) abgeleitet. Es ist aber selbstverständ- 
lieh, daß auch die zu (44) symmetrische Lösung 


4a) 


den Grenzbedingungen genügt; also gilt hier c=w(-—1]). 
Bei der Couetteschen Anordnung schließen wir daher aus 
(44) und (44a), daß zwei zueinander symmetrische Schwingungs- 
syteme existieren, deren Wellengeschwindigkeit bzw. mit den 
Geschwindigkeiten der beiden Wände (w (+1) und w (—1)) 
übereinstimmen. 

Für die symmetrische Strömung zwischen zwei ruhenden 
Wänden dagegen ist w(+1)=w(—1)=0. Aus (44) und 
an schließen wir r dann, daß jedes Integral « der Form 


den ie genügt. Aus der Forderung (a), daß 
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gm ungerade sein soll, ergibt sich aber, daß wir als unten 
Grenze des Integrals f ey 


w 


= 0 wählen müssen; also 


y 
dy 
w 


(44b) 


y=w 

0 
Im Falle symmetrischer Strömung gilt also insbesonder 
für B/a: 


(44c) 


w(+ 1} = w(—1)=£ 


Von der turbulenten Grundströmung interessiert uns zu- 
nächst am meisten die Art der Singularität derselben an den 
Wänden, also wenn wir annehmen w ~ 7 , der Exponent e. Wir 
wollen zu zeigen versuchen, daß aus den Differentialgleichungen 
(38) bzw. (39) im Limes R = oo wenigstens in unmittelbarer 
Wandnähe tatsächlich ein solches Potenzgesetz mit dem 
Exponenten e=}+ folgt. Dabei ist allerdings der Konvergenz- 
bereich der benutzten Potenzreihen nicht sicher gestellt, so 
daß die Schlüsse, soweit sie auf die Gestalt des Profils in 
einigem Abstand von der Wand angewendet werden, unsicher 
sind. Wir entwickeln 9, und w in der Umgebung von 7 =0 
nach ganzen und positiven Potenzen von 7 — dies ist mög. 
lich für jedes endliche — und dann umgekehrt 7 nach ganzen 
Potenzen von w. So werden wir direkt zu der Formel (48) 
für n (w) geführt. 


Wir behaupten — und dies ist das wichtigste Resultat, 


das wir später brauchen —g, läßt sich in 1. Näherung darstellen 
durch eine Reihe der Form 


Pi = + a5 + +..., 


WO %g... reelle, &;, a,,... rein imaginäre Konstanten sind; 
ferner ist w von der Form: 


(45) 8, 7n'+... 
Diese Behauptung läßt sich für die Differentialgleichungen (88) 
direkt durch Ansatz von » und w mit unbestimmten Koeffi- 
zienten beweisen, wenn die Glieder a, &, &, &%, und Ay Bı 
gegeben sind. Suchen wir also vorerst diese Glieder zu be 
stimmen. Zunächst muß 9, und 9,’ für n=0 wegen der 
Grenzbedingungen Null sein, also fängt die Reihe für 9, mit 
%n7 an (% =a, = 0). Daß ferner das folgende Glied a, 7? 


Grd: 


die 
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wegfällt, d. h. sehr klein wird gegen die übrigen Glieder, können 
wir nachträglich verifizieren. Einschaltungsweise wollen wir 
hier zu diesem Zwecke durch einmalige näherungsweise Inte- 
gration der zweiten Gleichung (38) beweisen, daß «, von der 
Größenordnung « R wird. Für «®—= 0 lautet Gleichung ey 


woraus folgt: 


Die Konstante A ist hierin von derselben Größenordnung, 
wie die linke Seite von (46) in der Mitte des Kanals, also 
höchstens von der Größenordnung 1 (vgl. I, § 2). 9,” am 
Rande ist daher wegen der Grenzbedingungen von der Ord- 


mung a R. Dasselbe gilt für «,. 


Nehmen wir also einstweilen «, als klein an und suchen 
dies später zu rechtfertigen. Von den Konstanten ß,, ß, ist 
die erste 8, wegen Forderung (a) $1 gleich Null. 

Die Konstanten «, und f, sind zunächst willkürlich!) und 
wir haben auch keine Möglichkeit, sie aus der Lösung der 
Differentialgleichungen (38) und (39) in Wandnähe abzuleiten. 
Diese Möglichkeit bestände erst, wenn es uns gelänge, die 
Lösung (45) analytisch bis zur anderen Wand fortzusetzen. 
Diese Aufgabe jedoch ist mathematisch außerordentlich kompli- 
ziert, schon deswegen, weil, wie sich zeigen wird, die ver- 
änfachten Gleichungen (38), (39) nicht ausreichen, um 9, und 
win der Mitte des Kanals zu bestimmen. Obwohl wir dem- 
nach auf die Lösung dieser Aufgabe verzichten müssen, können 
wir doch erwarten, daß wir durch die bloße Entwicklung von 
9, und w in der Nähe einer Wand mit unbestimmten Koeffi- 
menten a, 8, doch diejenigen qualitativen Eigenschaften von 
® und 9, in Wandnähe erhalten können, die erfahrungsgemäß 


ganz unabhängig sind von dem Verhalten der Flüssigkeit 


der Kanalmitte wie z. B. das Gesetz w = n'". 


1) Wir wollen «, als reell annehmen; dies bedeutet keine Be- 
schränkung der Allgemeinheit, denn ¢, ist nur bis auf einen Faktor der 
Form e'X bestimmt, da man in (36) den Anfangspunkt der Zeit beliebig 
wählen kann. 
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Gehen wir mit dem Ansatz: dit 

w= + Bau? + Ban? +: 

in Gleichung (88) ein, wobei wir die zweite Gleichung durch 
(46) ersetzen. Wir nehmen also wieder « als sehr klein 
was hier nur bedeutet (vgl. I, § 2), daß die Wellenlänge de 
Schwingungen groß sein soll gegenüber der u 


Außerdem setzen wir nach c) 


Für die erste Gleichung (38) schreiben wir ferner 
(47) — ta RG,’ — Fy — Pai Par) 
=w—Cy. 
Hierin bedeutet g,, den reellen, y,; den imaginären Tel 
von 
Aus (46) und (47) folgen jetzt die Rekursionsformeln 


(48) n(n — 1)(n — 2)«, o(n — 28) 
(49) nd, = 2« — _, 8, ; 


Hierin bedeutet «,, «,' den reellen bzw. imaginären Tel 
von &,. 


Aus (48) folgt zunächst 

Aus (49) ergibt sich dann: 


Bs = By = B, = By = 9. ob 
Auch f, können wir ‘veils Null setzen. 
Aus (49) folgt nämlich zuhören 


Pr 
das Glied f,7* ist also für sehr kleine n gegen das erste 
Glied ß,n zu vernachlässigen; für größere aber überwieges 
die höheren Glieder usw. völlig. 
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Nehmen wir also auch 6,=0 an und rechnen so die 


höheren Glieder der Reihen für $, und w aus, sh 
Es folgt 
ia @,=a,=0; a= — (« Ry? 


= — (aR) (5 ay” By a, 6,° ); 


2-7-9-10-11 8+5-6-8-9-10+11 


= = 0. 
B, = Bs = Po = = = By = 9; fox 


4 a,* a,* B, BIE. 
Die von uns behauptete Darstellung (45) für w ist also be- 
wiesen und es läßt sich leicht zeigen, daß auch von den weiteren 


Gliedern 6 nur jedes sechste Glied einen endlichen Wert hat. 
Für die Darstellung von n als Potenzreihe eat w dpc 


hieraus: 
(50) 4/1 = 47 ( ) qv ? Sch ifs 
= (@ R) 


Die Glieder y, bis Y,, y, bis 734 usw. sind alle gleich 
Null.?) 

Die Entwicklung (50) stimmt nun in der Tat völlig mit 
(48) überein und auch ohne Kenntnis der Konstanten «, und ß, 
scheinen wir so zu dem von Kärmän halbempirisch abge- 
leiteten Gesetz 7 ~ w* zu gelangen. Die Koeffizienten y, und y, 
lassen sich jedoch nicht berechnen. Umgekehrt können wir 
vielleicht aus dem empirischen Befund für die Koeffizienten ß, 
und x, schließen, daß y, von der Größenordnung 1, y, von der 
Ordnung (x R)-"*, also &, ~ (x R)”* ist. Nachträglich bestätigen 
wir so auch unsere frühere Behauptung «,n<«,. Fragt man 
sich nach der Größenordnung der Werte von w, für die das 


1) Diese Potenzreihe 7 (w) läßt sich natürlich auch direkt aus (46) 
und (47) ableiten, ohne den Umweg über die Reihe von w(n), wenn w 
als unabhängige Variable eingeführt wird. Doch sind die dazu nötigen 
Rechnungen etwas komplizierter. 
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8. Glied in (50) klein ist gegen das zweite, fiir die alg 
das w’-Profil wirklich gilt, so ergibt sich w ~ £,-‘*, alg 
= R-‘*, Aus den Differentialgleichungen (88) folgt nach 
alledem das Profil w ~7'* zunächst nur qualitativ; darüber, 
daß das }-Profil fast bis zur Kanalmitte hin beobachtet ist, 
geben unsere Rechnungen keinen Aufschluß. Das war jedoch 
auch nicht zu erwarten, denn die übrigen Konstanten, die 
in das Gesetz eingehen, hängen auch von dem Verhalten der 
Flüssigkeit an der gegenüberliegenden Wand ab. 

Fassen wir einschaltungsweise noch einmal kurz zusammen, 
welche Vernachlässigungen wir bei der Ableitung von (50) 
aus (48) und (49) begangen haben und suchen dadurch m 
bestimmen, mit welcher Genauigkeit die aus (50) gezogenen 
Schlüsse richtig sind. Zunächst haben wir das System (8) 
statt (37) benutzt, also Größen der Ordnung 9,/w gestrichen, 
Ferner haben wir a, = 0, B/@«=0, ß, = 0 gesetzt und damit 
Größen der Ordnung 


vernachlissigt. Die Genauigkeit unserer Rechnungen wird 


durch das größte der hier vernachlässigten Glieder bestimmt 
sein. Einfache Größenordnungsbetrachtungen, die wir hier 
nicht durchführen, machen es wahrscheinlich, daß von diesen 
Gliedern 9,/w das größte ist, daß aber auch dieses Glied mit 
R->oo gegen Null geht. Durch Wahl eines hinreichend großen 
Wertes für R wird man also die Genauigkeit der aus (48),( 49) 
und (50) abgeleiteten Resultate beliebig weit treiben können. 

Was das Blasiussche Widerstandsgesetz betrifft, so läßt 
sich dieses natürlich umgekehrt nach dem oben geschilderten 
Verfahren mittels Dimensionsbetrachtungen aus dem Gesetz 
n-w" herleiten, wenn man, wie geschehen, annimmt, daß 
das Verhalten von w in Wandnähe von der Kanalbreite un 
abhängig ist. Die Frage aber, ob auch diese letztere physikalisch 
sehr plausible Annahme aus den Differentialgleichungen (38) 
bzw. (39) folgt, müssen wir aus den oben angegebenen Gründen 
(Unmöglichkeit der analytischen Fortsetzung) unbeantwortet 
lassen. 

Eine bemerkenswerte direkte Folgerung über das Wider- 
standsgesetz können wir jedoch aus (38) und (39) durch Größen 
ordnungsbetrachtungen ziehen. Im Kanal, die unmittelbare 
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Wandnähe und den Punkt y=0 [vgl. u. (66)] ausgeschlossen, 
kann man nämlich statt der ersten Gleichung (38) wegen der 
Größe von C, schreiben (vgl. S. 612) 


(51) ia R Fi 91) 
Da nun die Amplitude 9, nicht mit R> oo gegen unendlich 
gehen kann — dies würde allen unseren Rechnungen den 


ay 


physikalischen Sinn nehmen —, so folgt, daß C, höchstens 
groB von der Ordnung « R ist, daß also der Exponent £ von 
Gleichung (40) < 1 sein muß (was ja auch in gewisser Weise 
aus (41) hervorgeht). Daraus folgt [vgl. Gleichung (42) und 
(40)|, daß das in der Hydraulik gewöhnlich angenommene 
Widerstandsgesetz r = const- U? eine obere Grenze für alle 
denkbaren Widerstandsgesetze der Turbulenz darstellt, die 
wnabhängig ist von der Wandbeschaffenheit. Man kann ver- 
mutungsweise schließen, daß das Gesetz t — U nur gültig ist 
für glatte Wände — nur für diese haben wir ja auch n — w’ 
erhalten —, daß sich aber das Widerstandsgesetz für rauhe 
Wände immer mehr dem quadratischen annähert.!) Denn für 
ruhe Wände wird die Amplitude gy, von R unabhängig und 
von der Größe der Wandstörungen sein — außerdem werden 
fir rauhe Wände die Grenzbedingungen nicht mehr ver- 
wsachen, daß 9, in erster Näherung reell ist, wie es der 
Gleichung (44) entspricht. 

An den Schlüssen dieses Paragraphen ändert sich nichts, 
wenn man an Stelle der Differentialgleichungen (38) die Glei- 
chungen (39) zugrunde legt. 

$3. Die turbulente Bewegung außerhalb der unmittelbaren 
Wandnihe. 

Fir die Bewegung in der Kanalmitte ist es wesentlich, 
daß sich 9, hier aus denjenigen beiden Integralen von (7a) 
zusammensetzt, die schon bei der reibungslosen Flüssigkeit, 
also für Gleichung (8) auftreten (vgl. I, $ 2). Die wichtigste 
hieraus folgende Eigenschafti:von 9, ist die, daß es, bis auf 
Größen der Ordnung g, und (« R)-! der Bedingung may 
(52) $1 — Fi = const. 
genügt. 

1) Vgl. die genaueren Untersuchungen von v. Karman a. a. O. 


und die experimentellen Untersuchungen von Schiller, dieselbe Zeit- 
schrift 8, S. 2, 1923. 
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Dies ergibt sich vermöge des Abelschen Satzes aus der 
Tatsache, daß, bis auf Größen der Ordnung 9, und (« R)-! g, 
und 9,; (der reelle bzw. imaginäre Teil von 9,) Lösungen de 
Differentialgleichung (8) sind. 

Daraus ist zu schließen, daß die Gleichungen (38) bzw, (89) 
nicht ausreichen, um die Bewegung in der ganzen Kanalbreite 
festzulegen, daß wir vielmehr zu Gleichung (37) und zu dem 
ihr im Couetteschen Fall entsprechenden Gleichungssystem 
zurückgreifen müssen. 

Dies bringt im allgemeinen eine Komplizierung der mathe 
matischen Aufgabe mit sich. Nur im Couetteschen Fall 
läßt sich die Aufgabe verhältnismäßig leicht lösen, weil die 
erste Gleichung (39) bis auf Größen der Ordnung w’/C algo 
(a R)-": und 9,2, vgl. (37), mit Gleichung (52) übereinstimmt, 
Während also Gleichung (52) ihrer Ableitung aus dem Abel. 
schen Satz gemäß nur bis auf Größen der Ordnung 9, richtig 
ist, soll hier im Couetteschen Fall noch (39) bis auf Größen 
der Ordnung 9,? gelten. 
leistet, wenn wir setzen: 


(58) 


Sister 

Diese Gleichung ist also als Lésung unserer zugrunde ge 

legten Differentialgleichung für den Couetteschen Fall an 


zusehen. 
Aus (53) würde für 9, aus (37) folgen 


(54) 


Nun gilt für den Couetteschen Fall nicht das System (37), 
sondern ein komplizierteres, das wir hier nicht anschreiben 
wollen. Wir geben aber von ihm an, daß es ebenso wie (37) 
für 9, = 0 zur Lösung (54) führt und daher zum Ergebnis: 


(55) 9, be-rY 


Hier sind a, b, y irgendwelche komplexe Konstanten. Fit 
w folgt dann aus der 2. Gleichung (39) bzw. aus ihrer redı. 
zierten Form (8): 

(56) w—c=a,er¥ +b, eM, 


Da nun w —c =O sein soll an der einen Wand, da anderer 
seits w ungerade sein soll um y = 0, so folgt, daß w einfach 
durch Nullwerden von y, und geeignetes Wachsen von a, und 
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Wir erhalten also das wichtige Resultat, daß beim Couette- 
schen Fall das Grundprofil w der turbulenten Bewegung in der 
anzen Kanalbreite wesentlich linear verläuft (— jedoch stark 
abweichend vom laminaren Profil; es wird viel flacher sein, 
als das laminare —) [daß es sich aber (vgl. II. § 2) am Rand 
wieder wie n'' an die Wände anschmiegt). 

Gehen wir nun zum komplizierteren Fall einer Strömung 
gischen zwei ruhenden Wänden, also genau zum System (37) 
üibr. Zur Lösung müssen wir uns naturgemäß mit rohen 
Näherungen begnügen. Zunächst können wir in (37) die 
rechten Seiten aller drei Gleichungen, nämlich die Reibungs- 
glieder, streichen — dies wird durch die Überlegungen von I, 
$2 vollauf gerechtfertigt. Dann setzen wir [vgl. Gleichung (44c)] 
pla = 0 


Für gy, erhalten wir so an Stelle der zweiten Gleichung (37) 
67) — — Py — F,)— 29, (F," 

— a? Gy) + (9 — 40? py) + — 40? = 0. 
Entwickeln wir p, als Lösung der Gleichung (57) noch Potenzen 
von «? einerseits, von g, andererseits, bemerken ferner, daß 
9, ungerade sein soll (vgl. 44b) und schreiben 9, = 90 + Pu» 
so ergibt sich bei en nur der linearen Glieder: 


d , ” s 

69) = + Ps Pio + 292 Pio 
0 0 


— 2910 Ps — Pio Pr ) 


Hierein ist 9, natürlich nur bis auf einen konstanten Faktor a 
bestimmt, den wir hier, was keine Beschränkung bedeutet, 
reell angenommen haben. 

Gehen wir mit diesem Werte von 9, in die vereinfachte 
1, Gleichung (37), nämlich: 


(60) Yı —Fi Gi = const. 
ein, so ergibt sich, wenn g,; den imaginären Teil von 9, be- 
zeichnet : 


y 
[2 [2 mn d 
(61) (Vai Pro +292: P10 —2Pi0 Poi — Pro Pai =const. 
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Jetzt genügt aber 9,,, wie aus der dritten Gleichung (87) und 


daraus, daß 9, in erster Näherung reell ist, folgt, der Gleichung 
(62) Po; W Po, 

0 


Gehen wir mit diesem Werte von 9,,; in (63) ein, bedenken wir 
ferner, daß für y= 0 die linke Seite von (61) und damit die 
Konstante der rechten Seite Null ist (dies bedeutet für die 
Konstante der rechten Seite von (60) nur, daß sie in erster 
Näherung Null, d. h. klein von der Ordnung 9, 93? bzw. a? 9,9, 
oder «* p, ist), so erhalten wir 


(64) fo — Pio Pio = 9; q 
was vollständig mit (54) übereinstimmt. — 

Diese Gleichung wird nun allerdings, eins wie (54), in 
der Umgebung des Punktes y = 0 trivial, sie ist, weil 9 eine 
ungerade Punktion von y ist, dort identisch erfüllt. Sie kam 
also dort auch keine Bestimmung von w liefern. Dies führt beim 
symmetrischen Profil (64) zu einer merkwürdigen Unstetigkeit 
im Punkte y=0. (Beim ungeraden Profil ist eine solche Up 
stetigkeit aus den Differentialgleichungen nicht zu erkennen), 
Intergriert man nämlich (87a), so erhält man, wie oben gezeigt, 
nach einmaliger Integration die Gleichung: 


(65) 2& RUG", Pr” Gr + 2 (G2 Pa G2) + 
[wo € (vgl. 8. 612) bei Gültigkeit des Blasiusschen Gesetz& 
von der Größenordnung (« R)**, also jedenfalls sehr groß ist} 
Die linke Seite von (65) verschwindet aber mit gy, und 9, (die 
ja ungerade Funktionen von y sind) an der Stelle y = 0. Dott 
muß also gelten 

(66) wo. = C. 


Dies bedeutet, daß w,_.” sehr groß ist (~ (« R)"«), daß alsow 
an der Stelle y= 0 einen scharfen Knick !) (Krümmungsradiw 


1) Auf diesen Knick hat mich freudlicherweise Hr. Prof. Prandil 
auf Grund des empirischen Materials hingewiesen. Der Knick scheint 
dabei empirisch nicht so scharf zu sein, wie ihn die Rechnung ergibt, wa 
sich wohl zwangslos dadurch erklärt, daß eben die Annahme (a) 8, 6@ 
über die Symmetrie auch der Wirbel und Störungen. nicht genau de 
Wirklichkeit entspricht. 
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~(« R)-**) aufweist. In geringer Entfernung von diesem 
Punkt muß w nach (64) wieder wesentlich linear verlaufen. 

Wir erhalten das Resultat: 

Auch bei der Strömung zwischen zwei ruhenden Wänden — 
was sich wohl auch auf die Strömung im Rohr übertragen läßt — 
verläuft das Profil näherungsweise über die ganze Kanalbreite 
linear, in der Mitte jedoch besitzt es einen 
scharfen Knick (an die Wände schmiegt Kajch 
6 sich mit dem y''-Gesetz an). (Vgl. 
Fig. 3.) 

Die physikalische Ursache des scharfen 
Knicks ist die, daß der Gradient der tur- 
bulenten Impulsübertragung für y=0 aus 
$ymmetriegründen verschwindet, daß also, 
wel der Gradient der gesamten Impuls- Fig. 3. 
ibertragung über die Kanalbreite kon- 
sant ist, der Gradient der laminaren Impulsübertragung, 
d.h. w’’ dort sehr groß sein muß. 


§ 4. Schlußbemerkungen und Zusammenfassung 5 | 
der physikalischen Ergebnisse. 


Unsere Untersuchungen haben noch zwei wesentliche 
Lücken: Erstens liefern sie nicht den Übergang von dem 
n*-Profil in das im mittleren Teil gültige lineare Profil 
umd zweitens beschränken sie sich auf große Werte von 
R, liefern also auch nicht den Minimalwert von R, falls ein 
solcher existiert, für den die turbulente Bewegung noch möglich 
ist. Die erste dieser beiden Lücken ist am schwersten aus- 
mfiillen (vgl. S. 617) und wir können auch kein Verfahren 
angeben, das in befriedigender Weise die hier gestellte Auf- 
gabe lösen würde. Man kann versuchen, die beiden von der 
Wand und der Kanalmitte herkommenden Näherungen an- 
tinander zu stückeln, was durch die Bedingung geschehen 
müßte, daß an der betreffenden Ansatzstelle 9,, 91’, 9, 9» 
wund w’ stetig sein sollen; die Konvergenz der Entwicklungen 
(45) und (50) reicht aber kaum aus, um so eine einigermaßen 
definierte Näherung zu garantieren. Jedenfalls hängt das End- 
tesultat, das Profil w, noch stark von der Art der Zusammen- 
setzung der beiden Näherungen ab. Es muß schließlich auch 
ils fraglich betrachtet werden, ob eine solche exakte Durch- 

Annalen der Physik. IV, Folge, 74. 42 
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führung der Ansätze (S. 609) wesentlich neue physikalische 
Resultate in Übereinstimmung mit der Erfahrung liefern würde, 
da diese Ansätze zweifellos eine sehr starke Idealisierung de 
Wirklichkeit bedeuten. 


Die Ausfüllung der zweiten Lücke dagegen bietet pringi. 
piell gar keine Schwierigkeiten, alle Hilfsmittel dazu sind 
schon im Teil I enthalten, und wenn einmal das Profil w voll 
ständig bekannt ist, so reichen die im Teil I geschilderten 
Methoden prinzipiell aus, um nach Art von I, §6 de 
Minimalwert von R zu berechnen, für den die turbulente 
Bewegung möglich ist. Man könnte z. B. die kritische 
Reynoldssche Zahl für ein Profil berechnen, das nach dem 
oben erwähnten Verfahren durch Anstückelung der beide 
Näherungen gewonnen ist oder man könnte für diese Unter 
suchung das empirisch beobachtete Profil zugrunde legen und 
so die Reynoldssche Zahl auf halbempirischem Wege aus 
rechnen. In jedem Falle wird man wohl — das haben die 
Untersuchungen im Teil I wahrscheinlich gemacht und direkte, 
hier nicht angegebene Rechnungen bestätigt — zur selben 
Größenordnung der kritischen Reynoldsschen Zahl, nämlic 
R ~10* kommen; der genaue Wert von R wird allerding 
noch zu sehr von der Art abhängen, wie das Profil gewonnen 
wurde, als daß wir Vergleiche mit der Erfahrung ziehen 
könnten. Deshalb haben wir hier noch keine derartige 
Berechnung von R durchgeführt. 


Fassen wir zum Schluß zusammen, was als physike 
lisches Ergebnis aus unseren Untersuchungen über das 
Turbulenzproblem zu schließen ist: Im Teil I haben wir er 
kannt, daß die laminare Bewegung und ihre - Stabilitits 
verhältnisse keine wesentliche Bedeutung für das Turbulenz 
problem und die kritische Reynoldssche Zahl besitzen. Im 
Teil II aber haben wir die turbulente Bewegung selbst unter 
sucht und können daraus einige Angaben über den turbulenten 
Bewegungszustand machen: Allgemein besitzt die Ge 
‚schwindigkeitsverteilung über den ganzen Kanal einfachsten 
Charakter, sie ist — je nach den Versuchsbedingungen — 
linear oder konstant (§ 8), in der Mitte besitzt sie bei der 
(symmetrischen) Strömung zwischen zwei ruhenden Wänden 
einen scharfen Knick; an den Wänden schmiegt sie sich im 
n'“Profil den Wänden an ($2). Darüber, daß das "-Profil 
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bis weit in das Kanalinnere hinein Gültigkeit besitzt, geben 
die Rechnungen keinen Aufschluß. Die turbulenten Schwin- 
gungen verlaufen beim Couetteschen Fall im Innern des 
Kanals nahezu harmonisch ($ 8, Gl. (58)), in der Nähe der 
Wände treten alle Oberschwingungen auf. Die Geschwindigkeit 
der Wellen stimmt mit der Wandgeschwindigkeit überein ($ 2 
Gl. 44—44c), beim Couetteschen Fall gibt es zwei Gruppen 
yon turbulenten Schwingungen, deren eine hinsichtlich ihrer 
Fortpflanzungsgeschwindigkeit mit der einen Wand überein- 
stimmt, während die andere die Geschwindigkeit der anderen 
Wand besitzt. Die turbulenten Störungen tragen also äußer- 
lich den Charakter einer Wandstörung. Dabei ist aber hervor- 
zuheben, daß diese Störungen unabhängig von Wandrauhigkeit 
und ähnlichen Einflüssen als freie Schwingungen existenz- 
fahig sind. Die Amplitude der turbulenten Wellen nimmt 
gegen die Wände (folgt aus (44) $ 2) stark zu, um erst in un- 
mittelbarster Wandnähe zu Null zu gehen. 

Die Wellenlänge der auftretenden Schwingungen (I, $ 8) 
ist der Größenordnung nach gleich, eher etwas größer als die 
Kanalbreite. Der Minimalwert der Reynoldsschen Zahl (I, § 8), 
für welchen noch die Turbulenz möglich ist, liegt größenord- 
nungsmäßig bei 10%. Als Widerstandsgesetz ($ 2, Gl. (50)) 
scheint sich für glatte Wände aus dem Profil 7 unter ge- 
wissen Voraussetzungen das Blasiussche r = u” zu er- 
geben. Für rauhe Wände nähert es sich wahrscheinlich ($ 2) 


liegenden Arbeit bildete weniger die Aufstellung dieser zum 
großen Teil schon bekannten Gesetzmäßigkeiten, als vielmehr 
der Nachweis, daß sich alle bisherigen, teilweise scheinbar sich 
widersprechenden Ergebnisse durch einfache Grundannahmen 
einheitlich mathematisch beschreiben lassen. 


Meinem hochverehrten Lehrer, Herrn Professor Sommer- 
feld, möchte ich für die Anregung zu dieser Arbeit und für 
vielfache Förderung den herzlichsten Dank aussprechen. 


München, Institut für theoret. Physik. id pet 
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2. Über die Ungültigkeit des Gleichverteilungssatsa 
bei Oszillatoren und Rotatoren; 
ah; von George Jaffe. 4 


Die Gültigkeit des Gleichverteilungssatzes ist im wesent 
lichen an die Vorzugsstellung geknüpft, die dem Energieintegnl 
gegenüber anderen Integralen der betrachteten Systeme ep 
geräumt wird. Eine solche Bevorzugung konnte berechtig 
erscheinen, solange man sich ausgesprochenermaßen auf di 
Ergodenhypothese berief, wie Maxwell und Boltzmann da 
in ihren grundlegenden Arbeiten stets getan haben; seitdem 
die Ergodenhypothese aber als ungeeignete Grundlage für 
statistisch-mechanische Untersuchungen erkannt ist, kann aud 
der Gleichverteilungssatz nicht in derjenigen Allgemeinhät 
als gültig angesehen werden, die häufig für ihn in Ansprud 
genommen wird. Diese Verhältnisse kommen in den kritischen 
Darstellungen von Ehrenfest!) und Hertz?) vollständig kla 
zum Ausdruck, es fehlt aber noch an konkreten Beispielen, 
bei denen sich auf Grund der Prinzipien und Methoden de 
klassischen statistischen Mechanik der Gleichverteilungssats 
als ungültig erweist. Im folgenden sollen zwei solche Beispiek 
gegeben werden; bemerkenswerterweise handelt es sich dabä 
um oft behandelte und einfache dynamische Modelle, bei denen 
die allgemeine Gültigkeit des Gleichverteilungssatzes stets aus 
_ driicklich angenommen worden ist. 


§ 1. Oszillatoren im Grenzfall eines großen » und tiefer 
Temperatur. 


Wir betrachten zuerst den Fall von Oszillatoren und 
denken uns dafür eine sehr ausgedehnte starre Hülk 
von zwei Sorten von Atomen erfüllt. Die Atome erste 


1) P. u. T. Ehrenfest, Math. Enc. IV. 32. 

2) P. Hertz in R. H. Weber und R. Gans, Repertorium de 
Physik. I. Bd. 2. Teil. 8. Buch. S. 436ff. Leipzig und Berlin 1916 
Man vgl. auch S. H. Burbury, Phil. Mag. (5) 50. 8S. 584. 1900 und 
G. H. Jeans, Phil. Mag. (6) 4. S. 585. 1902. _ 
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Art sollen die Masse m, haben und durch elastische Kräfte 
von der Intensität a an feste Ruhelagen gebunden sein; 
2, Yi» 2 seien die Komponenten der Entfernung eines solchen 
Atoms aus seiner Ruhelage. Die Atome zweiter Art sollen die 
Masse m, haben und frei beweglich sein. Die Anzahl der 
Teilchen beider Art soll sehr groß sein. Ob die Ruhelagen 
der Atome m, regelmäßig in Form eines Gitters angeordnet 
oder nach den Gesetzen des Zufalls über den ganzen zur Ver- 
figung stehenden Raum verteilt gedacht werden, ist für das 
Folgende ohne Bedeutung, jedoch soll der mittlere Abstand 
zweier Atome — gleicher oder verschiedener Art — groß 
gegen die in Frage kommenden Wirkungssphären sein. 

Unter diesen Umständen gelten nach Boltzmannschen 
Sätzen die folgenden Verteilungsfunktionen’) : 


m, . a 


fy (4a) 2,) = Age , (@,A,,A, konstant), 
aus denen die Gleichheit der mittleren kinetischen und poten- 
tiellen Energie pro Freiheitsgrad hervorgeht. Sieht man das 
gekennzeichnete System als angenähertes Modell eines starren 
Körpers an, so ergibt sich aus der Verteilungsfunktion (1) zu- 
gleich die Gültigkeit des Gesetzes von Dulong und Petit. 

Die angeführten Betrachtungen gelten scheinbar ohne Ein- 
schränkungen hinsichtlich der Größe der elastischen Kraft a, 
und so gilt es allgemein als ,,Folgerung aus der klassischen 
statistischen Mechanik“, daß ein oszillatorischer Freiheitsgrad 
beim Übergang zu immer größeren Werten von a den gleichen 
mittleren Energieinhalt aufweisen würde. Erst im Grenzfalle 
absoluter Starrheit müßte der Energieinhalt unstetig auf Null 
heruntergehen. 

Abgesehen davon, daß die Erfahrung einer solchen Auf- 
fssung widerspricht, muß sie auch logisch als wenig be- 
fedigend angesehen werden. Helmholtz hat mit allem Nach- 
druck darauf hingewiesen, daß der starre Körper nur als 
Grenzfall des elastischen Körpers aufgefaßt werden darf und 
hat in Erfüllung dieser Forderung eine Ableitung des Prinzips 


1) L. Boltzmann, Wien. Sitzungsber. 58. S. 517. 1868 (Ges. Abh. 
1, 8.49); 68. S. 731. 1871 (Ges. Abh. 1. 8.307). Vorlesungen über 
Gastheorie II. Anm. auf S. 126, 127. 
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der virtuellen Verschiebungen gegeben, bei der die sonst ühf gese 
lichen starren Verbindungen durch elastische Kräfte ersetg§ der 
werden.!) Die gleiche Forderung, die Helmholtz für dyn 
mische Betrachtungen erhoben hat, muß nun auch für staf. 
stische Untersuchungen geltend gemacht werden, denn die 
beruhen auf jenen. Wenn ein Körper sich dynamisch für a 
in Betracht kommenden Fälle mit hinreichender Genauigkei 
wie ein ideal starrer Körper verhält, so muß dasselbe ang 
in statistischer Hinsicht gelten, oder mit anderen Worte 
ideale Starrheit muß sich auch statistisch als stetiger Grenatall 
elastischer Deformierbarkeit ergeben. 


Wie kommt es aber dann, daß diese Forderung in dem oba 
angegebenen klassischen Beispiel verletzt wird? Die Antwor 
auf diese Frage kann nach unserer Auffassung nicht zweild 
haft sein: weil bei der Ableitung der Verteilungsfunktionen (j 
und (2) als einzige Bewegungsinvariante die Energie Berück 
sichtigung gefunden hat, und weil mit dieser einen Inv 
noch nicht genügend von dem dynamischen Charakter dé 
Systems zum Ausdruck gebracht ist. Die Energiegleichug 
läßt unter Umständen — hier z. B. für große Werte von a- 
Fälle statistisch noch zu, die dynamisch — durch andere Be 
wegungsintegrale — ausgeschlossen werden. 

In dieser Beziehung geben schon die Verhältnisse bei 
Ableitung des Maxwellschen Verteilungsgesetzes zu denken 
Setzt man die Gasmoleküle als ideal elastische Kugeln von ver: 
schwindendem Durchmesser voraus, nimmt man außerdem di 
Gültigkeit des StoBzahlansatzes und molekulare Unordnung 
an?), so ist das Ergebnis völlig bestimmt. Die alleinige Berück 
sichtigung des Energieintegrals liefert aber die richtige Ver 
teilung nur für den Fall, daß die Gasmasse als Ganzes ralt 
Die allgemeinste Verteilungsfunktion, auch für den Fall, dab 
eine translative oder rotatorische Massenbewegung vorhanda 
ist, erhält man nur unter Mitberücksichtigung der Schwef — 
punkts- und Flächenintegrale. Nun wird allerdings in diesem ® 
Falle der Charakter der Wärmebewegung durch die Masser 
bewegung nicht beeinflußt, wohl aber die Form des Verteilung 8) 


1) H. v. Helmholtz, Vorlesungen über die Dynamik diskreie 
Massenpunkte. § 61. 8.281ff. Leipzig 1911. 
2) P-'u. T. Ehrenfest, a. a. O. 
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gesetzes, und so zeigt sich auch schon hier die Wichtigkeit 
der Berücksichtigung aller Bewegungsinvarianten. 

_ Solange bei unserem Modell die von der elastischen Bindung 
herrührenden Kräfte klein sind gegen die beim Stoße wirkenden 
Kräfte, wird der Stoß zwischen einem Atom m, und einem 
Atom m, genau so erfolgen, als ob auch m, frei wäre. Dann 
ist das Ergebnis des Stoßes vollständig durch Energie- 
und Schwerpunktssatz bestimmt, und Formel (1) liefert die 
richtige Verteilungsfunktion. Sobald aber die Bindungskräfte 
mit den Stoßkräften vergleichbar werden, oder — was das- 
selbe besagt — sobald die Stoßdauer mit der Eigenperiode 
des Oszillators vergleichbar wird!), hängt das Ergebnis des 
$toßes auch noch von der elastischen Kraft a ab, und diese 
Abhängigkeit muß berücksichtigt werden. Genau wie bei der 
Aufstellung des Maxwell schen Verteilungsgesetzes muß man 
also bis auf den Elementarvorgang, d.h. in diesem Falle den 
Stoß zwischen einem gebundenen und einem freien Atom, 
mrickgehen, und die vollständige Behandlung würde die Auf- 
stellung sämtlicher Bewegungsinvarianten, also die vollständige 
Integration des dynamischen Problems erfordern. Eine voll- 
ständige Integration ist nicht durchführbar, handelt es sich 
doch gewissermaßen um ein Dreikörperproblem, bei dem die 
Störungskräfte mit der Zentralkraft vergleichbar werden. Wir 
werden uns darum in diesem Paragraphen auf den Grenzfall 
beschränken, daß die Bindungskraft groß gegen die Stoßkraft 
ist; denn gerade dieser Grenzfall ermöglicht hinterher den 
Übergang zur starren Bindung. 

Wir denken uns die Atome als Massenpunkte und setzen 
voraus, daß die Kräfte, die bei einem Stoß zwischen einem 
Atom m, und einem Atom m, wirksam werden, in die Richtung 
der Verbindungslinie fallen und Funktion des Abstandes 


Ak V(22 — 2) + (Ye — + (% — %)? 


snd. Die Bewegungsgleichungen lauten dann en A 
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632 G. Jaffe. 


nebst zwei weiteren Paaren analoger Gleichungen für die 
y- und z-Richtung. 
Um zu einer Lösung der Gleichung (3) zu gelangen, die 
näherungsweise für große Werte von a gilt, verfahren wir 
folgendermaßen. Beachten wir, daß das erste Glied auf de 
rechten Seite implizit eine der Zeit ist: 
(6) — =X, 3 
so läßt sich Gleichung (8) als eh Differentialgleichung mit 
rechtem Glied auffassen: 

az az, X(t) 
(6) 


die formal integriert folgenden Ausdruck für x, liefert: 


z + 2im, v (d 
Hier bedeuten C, und C, Konstanten, ki 4 


8) 

die Eigenfrequenz!) des Oszillators. Durch wiederholte par- 
tielle Integration gelangt man von (7) zu der gesuchten Nahe 
rungslösung 
xy XO, 
mv mw m, v® 
Es ist zu beachten, daß X (t) und seine zeitlichen Diffe 
rentialquotienten sich durch die sechs Koordinaten z, bis 4 
und die sechs Geschwindigkeiten #, bis z, ausdrücken lassen. 
Für X (t) folgt das aus seiner Definitionsgleichung (5); durch 
einmalige Differentiation erhält man 


“Ys — +  — — — — 


(9) xz, = Gert 4 


(10) 


1) Es ist zu beachten, daB in der vorliegenden Arbeit zur Verein- 
fachung des Satzes mit » durchweg das 2r-fache der Größe bezeichnet 
wird, die sonst häufig, z.B. in der Strahlungsformel, darunter ver- 
standen wird. Das Argument hv/kT des Planckschen Gesetzes lautet 
also in unserer Schreibweise hv/2n kT. 
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und indem man mit Differenzieren fortfährt und immer die 
auftretenden zweiten Differentialquotienten der Koordinaten 
durch ihre Werte aus den Differentialgleichungen (3), (4) usw. 
ersetzt, gelangt man zu der erwähnten Darstellung. 

Setzen wir von f(r) voraus, daß es mit allen seinen 
Differentialquotienten für große Werte von r verschwindet, 
so gilt dasselbe von X (t) für negativ und positiv große Werte 
von t. Ferner werden X (t) und seine Differentialquotienten 
während des Stoßvorganges um so größere Werte annehmen, 
je näher die Masse m, an die Masse m, herankommt, d.h. je 
größer die lebendige Kraft L, von m, vor Beginn des Stoßes 
war. Die Gleichung (9) wird also bei gegebenem » erst für 
hinreichend kleine Werte von L, eine brauchbare Approxima- 
tion liefern und bei gegebenem L, erst für hinreichend große 
Werte von v. Im allgemeinen gilt also die Reihenentwicklung (9) 
erst für hinreichend große Werte von » und hinreichend kleine 
Werte von L,. Sind aber diese Voraussetzungen erfüllt, so er- 
kennt man aus (9) und der oben angegebenen Eigenschaft 
von X (t), daß die Bewegung der Masse m, für Zeiten, die hin- 
reichend lange nach dem Stoße liegen, die ungestörte Fort- 
setzung der Bewegung vor dem Stoße ist. Es treten zwar 
während des Stoßvorganges Störungen auf, aber insgesamt 
wird weder die Energie noch die Phase der Bewegung von m, 
geändert. Ebenso erlangt dann m, schließlich wieder seine 
ursprüngliche lebendige Kraft. Es liegen also — in Überein- 
stimmung mit der Forderung von Helmholtz — Verhältnisse 
vor, bei denen der Rückprall von m, sehr ähnlich erfolgt, wie 
an einer festen Wand. Dabei kann m, selbst noch Schwingungen 
ausführen, deren Energie erst im Grenzfall » =o zu ver- 
schwinden braucht. 


Gehen wir nun von (9) zu reellen Werten über, so finden 


wir für x, und #, 
pr 
(11) +9 +20 (t) 


m, v? m, »* 


(12) #, =v A cos (vt + + 


wo A die Amplitude und @ eine Phasenkonstante ist. Denkt 
man sich nun nach dem oben Gesagten X (t) und seine Diffe- 
rentialquotienten durch die Werte der Koordinaten und Ge- 


die 
die 
wr 
Ä 
der 
> 
er, 
= 
- 
Zu 
3 
places 
- 
= 
iffe- 
sen, EIERN, 
rch 
yı) 
z,); 
hnet 
4) 
ver 


schwindigkeiten ausgedrückt, so läßt sieh t eliminieren. De 
so gewonnene Ausdruck 


| 2 4, m, »* 
| 


und zwei in y und z analog gebaute, stellen also Bewegung. 
invarianten für hinreichend große Werte von » und hinreichend 
kleine Werte von L, dar. 

Wir wollen nun im folgenden annehmen, daß für alle vor. 
kommenden Fälle die Glieder mit »~* eine ausreichende An 
näherung darstellen und wollen außerdem die drei gefundene 
Invarianten zu einer einzigen zusammenfassen.!) Schreiben 
wir noch i 


1) 


so erhalten wir demnach in 


einen Ausdruck, der neben der Gesamtenergie 


Sifedintet 

| +5 (#7 + 3,” + 2,7) — 


konstant sein muß. In diesem Grenzfall läßt übrigens der 
Ausdruck (15) eine besonders anschauliche Deutung zu. Die 
Bewegung des Atoms m, findet so statt, als ob es Pendel 
schwingungen um einen Nullpunkt ausführt, der um die Kom- 
ponenten 

X by Z 
m, m” 


verschoben ist, d.h. also in jedem Augenblick in Richtung 
der gerade wirkenden Kraft f(r) um den Betrag §=— is er 


1) Man vgl. hierzu die Ausführungen auf S. 639 und 640. 
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so daß sich die elastische Kraft a& und die Stoßkraft f(r) 
das Gleichgewicht halten. Die Beeinflussung der Bewegung 
von m, durch den Stoß mit m, ist also in dem betrachteten 
Grenzfall gewissermaßen als ‚statisch‘ zu bezeichnen. 

Wir wollen noch die Bedingung für die Gültigkeit der 
Näherungsformel (15) mit der ,,StoBdauer“ in Beziehung = | 
bringen. Die Dauer eines Stoßes ist bei der Annahme stetig = 
wirkender Kräfte ein Begriff, dessen scharfe Definition — rg = 
Schwierigkeit bereitet, genau wie etwa der Begriff der = 
Wirkungssphäre. Es ist aber einleuchtend, daß die Zeitdauer 
merklicher Wechselwirkung zwischen einem Atom m, und 
einem Atom m, von entscheidender Bedeutung für die hier __ 
untersuchten Verhältnisse ist. Man könnte etwa so vorgehen, 
daß man unter der Stoßdauer r diejenige Zeit versteht, 
während welcher die Kraft f(r) größer ist als ein bestimmter —_ 
Bruchteil des maximalen Betrages, den sie während des Stoß- RT 
vorganges erreicht. 

Für die folgende Betrachtung, die lediglich der Ab- 
schätzung einer Größenordnung dient, soll die einfache Vor- _ 
stellung zugrunde gelegt werden, daß die Komponenten der , a 
Kraft während einer ganz bestimmten Zeit r niherungsweise — 
konstante Werte X, F, Z haben. Bezeichnet man außerdem 
Mit Uy, Vg, We bzw. wg die Geschwindigkeitskomppo- 
nenten von m, vor und nach dem Stoße, so ist also ae 
(17) Xr=m, (u,’—u,), Yr = me, (% —v2), Zt=m, (w —w,). 


Der Grenzfall idealer Starrheit der Atome m, ist um so DEN ER 
näher realisiert, je geringer die während des StoBes uf mn 
übertragene Energie im Vergleich mit der lebendigen Kraft BR 
von m, ist. Diese übertragene Energie läßt sich auf Grund 
von (11), (12) und der analogen Gleichungen für die y- und 
#-Richtung berechnen; sie ist bei der Beschränkung uf Glieder 
von der Ordnung »~? nur potentieller Art, und ergibt sich Eee. ” 
Dabei ist der Mittelwert über alle Werte der Phasen- __ 


(19) E, + v,? + w,?) 
des stoßenden Meleküls — genommen.!) __ 


1) Man vgl. hierzu auch unten Formel (37). 
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Der von uns behandelte Grenzfall ist um so genauer ver 
wirklicht, je kleiner das Verhältnis 


z 


ist. Dann ist aber auch näherungsweise 


(21) =—wu,, = —w,, 
so daß man aus (17) bis (21) erhält dial; deat 
(22) 


Sofern also nicht der triviale Fall vorliegt, daß m, onl 
gegen m, ist, kann e nur klein sein, wenn vr groß, d.h. die 
Stoßdauer t groß gegen die Schwingungsdauer 2r/» ist. 
Wir gehen nunmehr dazu über, eine Zustandsverteilung 
für unser dynamisches Modell aufzustellen, die sich stationär 
aufrecht erhält, wenn bei allen Stößen zwischen Atomen ver- 
schiedener Art die Voraussetzungen für die Konstanz von I 
erfüllt sind.!) Der mittlere Abstand zweier Atome soll so groß 
gegen die Wirkungssphäre sein, daß nur binäre Stöße berück- 
sichtigt zu werden brauchen. 

Da die Invariante J sich nicht additiv aus zwei Bestand- 
teilen zusammensetzt, die einzeln von den Koordinaten und 
Geschwindigkeiten je eines Atoms abhängen, lassen sich auch 
keine Verteilungsfunktionen nach der Art von (1) und (2) an- 
geben, bei denen die Häufigkeit einer bestimmten Phase des 
einen Atoms unabhängig von der Phase der Atome anderer 
Art ist. Man kann aber genau so verfahren,.wie Boltzmann 
zu einer stationären Zustandsverteilung für zusammengesetzte 


1) Nach dem Maxwellschen Verteilungsgesetz kommen Geschwin- 
digkeiten beliebiger Größe — wenn auch entsprechend selten — vor, so 
daß bei strenger Gültigkeit des Maxwellschen Verteilungsgesetzes für 
die Atome m, die Voraussetzung der Konstanz von / nicht bei allen 
Stößen zutreffen könnte. Es ist aber zu bemerken, daß das Max wellsche 
Gesetz selbst nur ein Grenzgesetz für eine unendlich große Zahl von 
Atomen ist; wenn die Gesamtzahl der letzteren und die Gesamtenergie 
endlich sind, kann auch die Energie des einzelnen Teilchens nie über alle 
Grenzen steigen. Im übrigen genügt es, durch Vergrößerung von » und 
Verkleinerung der mittleren kinetischen Energie die Ausnahmefälle auf 
einen zu vernachlässigenden Bruchteil herunterzudrücken, da nur ein 
Grenzgesetz aufgestellt werden soll. 
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Gasmoleküle gelangt, wenn Wechselwirkung zwischen je zwei 
Molekülen besteht.!) 

Wir setzen nämlich die Zahl der Paare von Atomen ver- 
schiedener Art, für welche das Atom erster Art in der Phasen- 
zelle 


(28) do, = dz, dy, dz,diz,dy,dz,, 

und das Atom zweiter Art in der Phasenzelle ke or 
liegt, in der folgenden Weise an: | 
(25) adN,, = 


wo A,., « und # konstant sind. Genau wie an der zitierten 
Stelle bei Boltzmann, läßt sich nach der Kirchhoffschen 
Methode beweisen, daß die Verteilung (25) durch die Stöße 
zwischen Atomen verschiedener Art nicht gestört wird. Es 
braucht also nur noch der Nachweis erbracht zu werden, daß 
dasselbe auch für die Stöße zwischen Atomen m, gilt. 

Um zu dem Verteilungsgesetz für die Atome m, zu ge- 
langen, hat man nur den Ausdruck (25) nach den Differen- 
tialen (23) zu integrieren. Nun ist nach unseren Voraussetzungen 
der Bruchteil der Atome zweiter Art, die jeweils in merklicher 
Wechselwirkung mit Atomen erster Art stehen, verschwindend 
klein; es genügt also zur Ausführung der angegebenen Inte- 
gration solche Paare zu betrachten, die nicht gerade in Wechsel- 
wirkung befindlich sind. Für diesen Fall reduziert sich aber 
(25) bei passender Verfügung über die Konstante in 9 (r) auf 


gesetzt wird. Führt man nun in (26) die Integration nach do, 
aus, so ergibt sich für die freien Atome die Maxwellsche 
Verteilung, die bekanntermaßen durch Stöße der Atome m, 
unter sich nicht gestört wird. 


1) L. Boltzmann, Vorl. üb. Gastheorie. II. IV. Abschn., §§ 36—40. 
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Es ist damit gezeigt, daß die durch (25) dargestellte Zu. 
standsverteilung stationär ist. Es muß aber mit besonderem 
Nachdruck betont werden, daß die benutzte Methode hier, 
wie bei der Anwendung auf zusammengesetzte Moleküle, nur 
eine hinreichende, keineswegs aber eine notwendige Bedingung 
für stationären Zustand liefert, und auch das nur, insoweit 
das System der berücksichtigten Invarianten als erschöpfend 
angesehen werden kann. Die vorliegende Untersuchung ver- 
folgt aber auch nur den Zweck zu zeigen, daß die Erfüllung 
der Helmholtzschen Forderung mit den Methoden der 
klassischen statistischen Mechanik erreichbar ist. 

Schon die Gleichung (26) läßt erkennen, daß der Gleich. 
verteilungssatz mit unserem Ansatz in Widerspruch steht!), 
und daß die statistische Unabhängigkeit der Energie der 
Oszillatoren von der mittleren lebendigen Kraft der freien 
Atome als Folge der vorausgesetzten dynamischen Unab- 
hängigkeit in die Erscheinung tritt. Wir wollen aber die Be 
rechnung der Mittelwerte an die genauere Gleichung (25) an- 
knüpfen und konsequent Größen von der Ordnung »~? bei- 
behalten. Es ist dabei zu berücksichtigen, daß sich wegen des 
Faktors a = m,v? die Ordnung einiger Glieder in (15) er 
niedrigt; andererseits sind 2,, y,, 2, als Größen anzusehen, 
die mindestens von der Ordnung »-! klein sind, so daß auch 
keine weiteren Glieder aus (3) berücksichtigt zu werden 
brauchen. 


Führt man noch die Substitution a 
es) | ©, Yo 
AG — 2, =rsin dsiny 


aus, so lassen sich sämtliche Integrationen bis auf diejenige 
nach r ausführen. Mit der angegebenen Genauigkeit findet man 


28 
2a +8) 


Ser" (dr 
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(go) 


1) Weil m, &,?/2 und m,&,?/2 im Exponenten von (26) mit ver- 


schiedenen Faktoren multipliziert auftreten. Kr 
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f(r) 
(mera) 

Demnach ist der Mittelwert der zur Koordinate z, ge- — 
hörigen lebendigen Kraft für solche Paare, die den Abstand r 
haben, 


(32) | 2 )= 2(a +B) + + (r+ 
und die gleiche Formel gilt als allgemeiner Mittelwert, wenn 


man unter r auf der rechten Seite einen passend gewählten 
Mittelwert versteht. 


(31) 


In gleicher Weise ergibt sich 
2], 2%@+) (a + 8) m, »? 


auch diese Formel kann als Gesamtmittelwert aufgefaßt 
werden, dabei bedeutet dann aber im allgemeinen r einen 
anderen Mittelwert als bei der entsprechenden Auffassung der 
Formel (32). 

Gleiche Formeln wie (32) und (33) gelten für die y- und 
zKomponente; hierzu ist jedoch folgendes zu bemerken. Wir 
haben auf S. 634 drei Invarianten zu einer einzigen zusammen- 
gezogen und bei Aufstellung des statistischen Gesetzes nur 
diese eine berücksichtigt, was scheinbar im Widerspruch mit 
unserer Forderung steht, daß alle Invarianten berücksichtigt 
werden müssen. Hätten wir nun die drei unabhängigen In- 
varianten beibehalten, so würden wir jetzt gefunden haben, daß 
die mittleren kinetischen und potentiellen Energien für die drei 
Koordinatenrichtungen voneinander unabhängig sind. 

Dieses Ergebnis steht vollständig in Übereinstimmung mit 
unseren dynamischen Voraussetzungen. Denken wir uns etwa 
einen Anfangszustand, bei dem die Atome m, räumlich von den 
m, getrennt sind und einen gewissen Energieinhalt haben; 
dann mögen zu irgendeiner Zeit die m, mit den m, in Wechsel- 
wirkung treten, und es werde nach dem schließlich sich ein- 
stellenden stationären Zustand gefragt. Ist die Gesamtenergie 
der m, hinreichend niedrig und » hinreichend groß, so behält 
jedes Atom m, für alle Zeiten in jeder der drei Koordinaten- 
fiehtungen die Energie bei, die es zu Beginn hatte, und ein 
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Ausgleich erfolgt nur zwischen den m, unter sich. War algo 


für alle Zeiten bleiben. Der innere Grund dafür ist, daß bei 
einem harmonischen Oszillator die Energien der Bewegung in 
den drei Koordinatenrichtungen Einzelinvarianten sind. 

Indem wir also die drei auf 8. 634 gefundenen asym. 
ptotischen Invarianten zusammenzogen, haben wir lediglich 
Fälle von Anisotropie in den Anfangsbedingungen ausge. 
schlossen. 

Für die mittleren kinetischen Energien der freien Atome 
erhalten wir mit der angegebenen Genauigkeit 


(84) 3 2” 
Es sind nun noch die Konstanten A,, « und ß zu be 
stimmen.!) A,, hängt in bekannter Weise mit der Gesamt- 
zahl der Atome zusammen; wegen der Beziehungen (84) 
werden wir 


setzen müssen, wo k die Boltzmannsche Konstante ist, und 
6 bestimmen wir so, daß der Helmholtzschen Forderung 
genügt wird. 

Die Formeln (32) und (33) zeigen, daß unter unseren 
Voraussetzungen die mittlere kinetische Energie der Atome m, 
von ihrer mittleren potentiellen Energie verschieden ist und 
daß beide vom Gleichverteilungswerte 1/2«& abweichen. Sie 
lassen sich durch Vergrößerung von ß beliebig verkleinern, 
und der Helmholtzschen Forderung kann jedenfalls genügt 
werden, wenn man annimmt, daß lim 6 = oo ist. Setzen wir 


also in den Gleichungen (32) und (88) ß = oo, so wird 


az, f* (r) 
87) 


Es ist also dann die Energie, die von den Oszillatoren 
aufgenonmmen wird, nur potentieller Art; die kinetische Energie, 


1) Man vgl. hierzu auch die Ausführungen auf S. 654, 655. 
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daß sie sich stationär erhält, aber anisotrop, so würde sie das | 


be 
die 
| Ordı 
Ene 
blei 
auß 
Seh 
and 
Kre 
ia per 
Ate 
und 
|, > 
Gre 
dal 
& 
} 
rat 
lin 
ve 
| 
be 
de 
di 
di 
la 
: 
ten 
g 
N 
p 
| 


oren 
rgie, 


| die notwendigerweise mit den zeitlichen Änderungen von 
| 2% Ya» 2, verknüpft ist, erweist sich als klein von höherer 
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| wegung der auf sie stoßenden freien Atome. Auch dieser Rest 
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Ordnung. In der Tat ist mit der Annahme 8 = oo über die 
Energie der Oszillatoren in einem wie oben (8. 639) charakteri- 
sierten Anfangszustande so verfügt, daß sie verschwindet; sie 
bleibt dann auch im stationären Endzustande dauernd gleich Null, 
außer während der Perioden der Wechselwirkung, wenn f(r) + 0 
ist. Unter diesen Umständen führen die Oszillatoren gar keine 


von Beeinflussung verschwindet einerseits mit zunehmendem », 
andererseits aber auch mit abnehmender mittlerer lebendiger 
Kraft der stoßenden Atome, d. h. also mit abnehmender Tem- 
peratur. Denn je kleiner die lebendige Kraft eines stoßenden 
Atoms ist, desto größer ist der Wert von r, bei dem es umkehrt, 
und um so kleiner ist dementsprechend im Mittel der Wert 
von f (r). 

Unsere bisherigen Betrachtungen, die auf den extremen 
Grenzfall »r>1 beschränkt sind (vgl. S. 636), beweisen nur, 
daß man der Helmholtzschen Forderung durch die An- 
nahme lim $ = oo genügen kann, ohne mit den Grundlagen 


der klassischen statistischen Mechanik in Widerspruch zu ge- 
raten. Der Beweis, daß auf Grund der gleichen Prinzipien 
im f = oo sein muß, gelingt erst durch Aufstellung einer um- 
‚m=00o 
fassenderen Invariante und wird im dritten Paragraphen er- 
bracht werden (man vgl. S. 655). 
$ 2. Rotatoren im Grenzfall kleinen Trägheitsmoments Kar 
und tiefer Temperatur. 


Ähnliche Verhältnisse wie bei Oszillatoren, liegen auch 
bei Rotatoren vor, wenn man die Umlaufszahl bei diesen mit 
der Schwingungsfrequenz bei jenen in Parallele setzt. Solange 
die Dauer eines Stoßvorganges klein gegen die Umlaufszeit 
des Rotators ist, kann ein vollständiger Austausch der Energie 
mit den freien Molekülen stattfinden; wenn aber die Um- 
laufszeit klein gegen die Stoßdauer ist, läßt sich voraussehen, 
daß die Rotationsenergie näherungsweise unbeeinflußbar durch 
Stöße wird, weil dann während des gesamten Stoßvorganges die 
positiven und negativen Drehimpulse sich annähernd aufheben. 


der Physik. IV. Folge. 74 
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Es besteht der wesentliche Unterschied gegen den Oszil. 
lator, daß die Umdrehungszahl nicht wie die Eigenfrequenz 
ein Systemparameter ist; als solcher kommt vielmehr nur das 
Trägheitsmoment in Frage. Nach der üblichen Auffassung 
müßte die mittlere Energie, die zu einem rotatorischen Freiheits- 
grade 9 gehört, den normalen Betrag k T/2 haben, unabhängig 
vom Werte des Trägheitsmoments J; nach dem soeben Aus. 
geführten sind jedoch Abweichungen von der Gleichverteilung 
zu erwarten für solche Werte der Temperatur und des 
Trägheitsmoments, für welche die aus der Gleichung 
J# kT 
2 


(38) = 
berechnete Umlaufsdauer klein gegen die Stoßdauer ist. Dabei 
hängt natürlich auch die mittlere Stoßdauer von der Tempe 
ratur ab (vgl. 8. 645). 

Wir beschränken die folgenden Rechnungen auf den ein- 
fachsten, zweidimensionalen Fall, der das prinzipiell Wichtige 
bereits hervortreten läßt. Wir denken uns also in einer Ebene 
Massenpunkte m,, die nur im festen Abstand @ um gegebene feste 
Punkte rotieren können; außerdem gebe es Massenpunkte m, 
die frei beweglich in der Ebene sind. Wir untersuchen einen 
Stoß zwischen einem Atom m, und einem Atom m, und be 
schränken uns auf den Fall, daß der Rotator vor Beginn des 


Stoßes bereits eine sehr große Umlaufsgeschwindigkeit d—=» hat, 


Das Rotationszentrum O von m, machen wir zum Ur- 
sprung eines rechtwinkligen Koordinatensystems und ver- 
wenden außer den kartesischen Koordinaten 2,, Yı, Tg Ih 
Polarkoordinaten o, 9 für m, mit O als Pol, und Polarkoordi- 
naten r, y für m, mit m, als Pol. Es sollen also die folgenden 
Beziehungen gelten 


=008d, y=osind, = const; 4 


(39) ı=r08Y %—Yızrsiny, 
Ve; — 2)? + (% — yı)* 


_ Mit diesen Bezeichnungen lauten die Bewegungsgleichungen 
(man vgl. die Figur) 


(40) 


d’ x, F 
= f(r) cos w, m Ta = sin w, 
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wo » die oben definierte Konstante ist. 
stitution erhalten wir aus (41) für 9, die Differentialgleichung 


(44) 


Indem wir nun hier, ähnlich wie oben (vgl. 8. 632), f(r), 9, und 
p auf der rechten Seite als Funktionen der Zeit auffassen, 
erhalten wir durch schrittweise partielle Integration unter 
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sin 5 = sin (9 — y) - 


Um außer der Energie 
» My 


J 


+) - 
gr) = [fear 


einen Ausdruck zu gewinnen, der unter unseren Voraussetzungen 
näherungsweise eine Stoßinvariante darstellt, setzen wir 


(8) 


0 


Fig. 


jon 


dt? 


= of(r) 


Berücksichtigung von (48) 


r) cos (# 


+ COs (f(r) sin (#, — 


648 


Durch diese Sub- 


r)sinvtcos(#, — w)+ of {r) cos vétsin (9, — w). 
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Uber die Integrationskonstante ist bereits so verfügt, daß yor 
Beginn des Stoßes dd#/dt=» ist, und die Formel (45) läßt 
erkennen, daß innerhalb ihres Gültigkeitsbereiches nach dem 
Stoß die Umlaufsgeschwindigkeit schließlich wieder ihren ur 
sprünglichen Wert » annimmt; nur während der Periode der 
Wechselwirkung findet eine vorübergehende Störung statt, 
Die Approximation gilt — analog wie oben (vgl. S. 688) — für 
hinreichend große Werte von » und hinreichend niedrige Werte 
der lebendigen Kraft des stoßenden freien Atoms. 

In sämtlichen Gliedern der rechten Seite läßt sich 

wiederum das explizite Vorkommen von ¢ eliminieren. Wir 
beschränken uns im folgenden auf Glieder, die von der 
Ordnung »~? gegen das Hauptglied sind, und gewinnen somit 
die gesuchte Grenzinvariante 
(46) T= 2 5 + 9 f(r) cos (9 — y) = Ir? = const. 
Im Anschluß an die Gleichung (45) läßt sich in genau 
derselben Weise, wie bei Oszillatoren (vgl. S. 635), der Nach- 
weis erbringen, daß die hier gewonnene Näherung um so ge 
nauer gilt, je kleiner die Umdrehungszeit 27/v gegen die Stoß- 
dauer T ist. 

Indem wir jetzt von der dynamischen zur statistischen 
Betrachtung übergehen, setzen wir für die Zahl der Paare von 
Atomen verschiedener Art, die eine bestimmte Phase haben, 


Unter E und I sind naturgemäß hier die Ausdrücke (42) 
und (46) zu verstehen; alle darin auftretenden Variablen sind 
durch #, x, und y, ausgedrückt zu denken. Es läßt sich genau 
wie früher zeigen, daß die Verteilung (47) stationär ist; für 
die freien Atome gilt das Maxwellsche Gesetz. 

Scheidet man wieder diejenigen Paare aus, die gerade 
in Wechselwirkung stehen, so erhält man 


(48) dN., = Aye dddtd zx, di, dy, dj,. 
Von diesem Ausdruck unterscheidet sich (47) durch Glieder 
im Exponenten, die nur Lagekoordinaten enthalten. Dem- 
nach erhält man für die Energiemittelwerte ohne weitere 
Rechnung 
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2 2 2a’ 


und zwar gelten diese Ergebnisse bis auf Zusatzglieder von der 
Ordnung »~*. 

Im Falle des Oszillators haben wir die Konstante 8 aus 
der Forderung bestimmt, daß im Grenzfalle » = oo die mittlere 
Energie des Oszillators verschwinden soll, Diese Forderung 
kann beim Rotator natürlich nicht erhoben werden, denn bei 
ihm steigt die Energie mit » ins Unendliche; unsere Betrach- 
tungen gelten jetzt auch gar nicht für den Grenzfall » = oo, 
sondern lediglich für den Grenzfall»r = w. Es muß darum 
auch zunächst festgestellt werden, unter welchen physikalischen 
Bedingungen dieser Fall vorliegt, und dafür muß wenigstens 
eine Schätzung der StoBdauer auf dynamischer Grundlage ver- 
sucht werden. 

Untersucht man nun den Stoßvorgang zwischen einem 
freien und einem starr gebundenen Atom auf Grund eines 
Kraftgesetzes f(r) = yr-" und definiert die Stoßdauer so, wie 
wir es auf 8.635 getan haben, so erweist sich für zentrale 


Btößer —L, , wenn L, wieder die kinetische Energie 


des freien Atoms bedeutet. Demnach ist imMitteler—T (3 ; 
da andererseits in demjenigen Temperaturbereich, wo Gleich- 
verteilung herrscht, »~ 7"* ist, liegt der klassische Fall 
yt <1 bei hohen Temperaturen, und der von uns untersuchte 
Grenzfall bei tiefen Temperaturen vor. Dabei werden Ab- 
weichungen von der Gleichverteilung und somit auch der von 
uns behandelte Grenzfall bei um so höheren Temperaturen 
eintreten, je kleiner das Trägheitsmoment des Rotators ist. 
Unsere Formeln sind mit der Annahme’vereinbar, daß 
eine Nullpunktsenergie E, bei Rotatoren existiert; man braucht — 
dafür nur für ihren Giiltigkeitsbereich 
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zu setzen. Dabei soll E, nicht von der Temperatur abhängen, 
kann aber möglicherweise noch Funktion des Trägheits 
moments J sein!) Dynamisch ist die Nullpunktsenergie x 
zu erklären, daß bei der Abkühlung eines Körpers, der Ro. 
tatoren enthält, schließlich eine Temperatur erreicht wird, bei 
der die Stöße zu langsam erfolgen, um den Rotatoren den 
Rest ihrer Energie zu entziehen. 


§ 8. Lineare Oszillatoren für beliebige Werte von » und 7. 


Die Ausführungen der beiden ersten Paragraphen zeigen, 
daß bei Oszillatoren und Rotatoren auf Grund der dynamischen 
Verhältnisse Abweichungen vom Gleichverteilungssatze be 
tiefen Temperaturen zu erwarten sind, wenn die endliche 
Dauer der Stöße berücksichtigt wird. Die Abweichungen 
müssen bei um so höheren Temperaturen einsetzen, je größer 
die Eigenfrequenz der Oszillatoren und je kleiner das Trägheits- 
moment der Rotatoren ist. Unsere bisherigen Rechnungen er. 
geben aber nur eine asymptotische Lösung für den Grenz 
falvr>1. 

Es ist einleuchtend, daß ein stetiger Übergang zwischen 
unseren Lösungen (vr > 1) und der normalen Gleichverteilung 
(vt <1) existieren muß; aber dieser Übergang würde sich 
nur durch vollständige Integration der dynamischen Pro 
bleme (3), (4) und (40), (41) gewinnen lassen. Es müßten die 
jenigen allgemein gültigen Bewegungsinvarianten aufgestellt 
werden, die für große Werte von vr in unsere Gleichungen (18) 
bzw. (46) übergehen. Da diese Aufgabe an mathematischen 
Schwierigkeiten scheitert, wollen wir, um wenigstens einen 
qualitativen Überblick für beliebige Werte von »t zu ge 
winnen, den Fall des Oszillators im folgenden unter verein 
fachenden rechnerischen Annahmen behandeln. 

Ein Blick auf die Gleichungen (7) und (4) lehrt, daß das 
dynamische Problem integrabel wird, wenn die beim Stoß wirk- 
samen Kraftkomponenten als Funktionen der Zeit bekannt sind. 
Wäre das Problem gelöst, so ließen sich diese Zeitfunktionen be 
stimmen. Wir wollen umgekehrt annehmen, daß die Kräfte explizit 


1) Wollte man es z.B. als erwiesen ansehen, daß Rotatoren eine 
vom Betrage hv/4n haben (val. Anm. S. 
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als Funktionen der Zeit gegeben sind und wollen über diese ee 
Funktionen geeignete Annahmen machen. Zur weiteren Ver- __ 
einfachung der Rechnung wollen wir ein noch elementareres _ 
dynamisches Modell als im ersten Paragraphen betrachten, 
indem wir die Bewegungsfreiheit auf eine Koordinate be- 
schränken. An dem so entstehenden ‚linearen Körper‘‘ treten 
die prinzipiell wichtigen Punkte bereits zutage, es sei aber 
bemerkt, daß diese Vereinfachung für die Durchführbarkeit 
der Rechnung nicht wesentlich ist. 

Wir nehmen also nunmehr an, die elastisch gebundenen 
Atome m, und die freien Atome m, seien nur auf einer Geraden 
beweglich, die wir zur z-Achse wählen. Der mittlere Abstand 
sei, wie bisher, groß gegen die Wirkungssphäre, und zwischen 
je zwei gebundenen Atomen mögen sich viele freie befinden. 
Unsere früheren Gleichungen (3) bis (8) bleiben bestehen, 
wenn man sich r durch |x, — 2,| ersetzt denkt. 

Die beim Stoß wirksame Kraft denken wir uns in der 
Form 


(53) X() = X, X,9(+ ) 

gegeben. X, soll eine Konstante sein und g eine gerade Funk- 
tion!) der dimensionslosen Zahl t/r, die für Werte +|>ı 
mit wachsendem |= | schnell gegen 0 abfällt und integrabel — 


ist. Damit ist t als ein Maß der Stoßdauer charakterisiert. \ 
Ferner sei noch über X, und g so verfügt — was immer mög- 
lieh ist —, daß 


1) Die Beschränkung auf gerade Funktionen ist wiederum nicht 
wesentlich; die Formeln werden aber etwas schwerfälliger, wenn man diese 
Annahme fallen läßt. Wir haben ferner mit dem Ansatz (53) so über den 
Nullpunkt der Zeit verfügt, daß er mitten im Stoß angenommen wird; 
auch diese Einschränkung ist unwesentlich und kann durch eine Sub- 
stitution t — t, behoben werden. Wher 
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G. Jaffe. ’ 
wird. Dann ist X, ein Maß für die Größe der Kraft während 
des StoBes. 

Die Stoßdauer ist natürlich nicht als Konstante, sondern als 
Funktion der Anfangsbedingungen zu denken; als solche können 
angenommen werden die Einzelenergien E, und E, der beiden 
Massenpunkte, die Lage von m, und die Phase von m, zu irgend- 
einem Zeitpunkte vor dem Stoße. Die Stoßzeit darf aber nicht 
als beliebige Funktion der Anfangsdaten angesehen werden, 
Der umgekehrte Stoß, der durch Umkehr aller Geschwindig. 
keiten in irgendeinem Zeitpunkte nach dem Stoße entsteht, 
muß dynamisch ebensogut möglich sein, wie der direkte, und 
muß die gleiche Stoßdauer ergeben. Demnach darf r nur als 
Funktion solcher Kombinationen der Anfangsdaten an- 
gesehen werden, die beim Stoß ungeändert bleiben; oder mit 
anderen Worten: die Stoßzeit ist selbst eine Stoßinvariante, 

Wir wollen mit der Gleichung (7) zu reellen Werten über- 
gehen; dann wird unter Einführung der Amplitude A und 
einer Phasenkonstante ®: 


sin”! (cos vt X(t)dt 
vm, 


cos vt 
>= fain vtX (edt. 


Substituiert man hier (53) und unterscheidet die Variabeln 
nach dem Stoß von denen vor dem Stoß durch einen Akzent, 
so erhält man unter Berücksichtigung des Umstandes, daß g 


gerade vorausgesetzt wurde, x 
(57) g, = Asin (yt + As 
x7 (vr) sin vt 


(56) z = Asin (vt + + 


vor dem StoB, und 
(58) 2,’= A’sin(vt + 3’) = Asin(vt + 9)+ 


nach dem StoB. eg 
Dabei ist nach Ersetzung von t/t durch s 
+» 
(59) zer) = [eos (vr+s)g(s)ds 
gesetzt. Die Integration in (59) ist streng genommen nut 
von einem negativ großen Werte von t/t („vor dem Stoß“) 


bis zu einem positiv großen Werte („nach dem Stoß‘) zu er- 
strecken; dafür können aber ohne Bedenken die Grenzen — @ 
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und +o gesetzt werden. Von der Funktion x (yr) läßt sich 
aussagen, daß sie im lim yt = 0 gegen 1 konvergiert, und im 
jim yt = 00 selbst unendlich wird. Erstere Eigenschaft folgt 
aus der Festsetzung (54), letztere aus der Tatsache, daß 
imy'(er)=0 ist, wie entweder aus der früher durch- 


geführten Betrachtung dieses Grenzfalles hervorgeht!), oder 
aus der Anwendung eines bekannten Satzes von Dirichlet 
auf die Gleichung (59). 

Durch Differentiation von (57) und (58) erhält man Are 


(60) #, =v A cos (vt + 8), 
(61) cos (vt + = Acos (vt + + = vr) cos vé. 


Nun gelten die Gleichungen (58) und (61) für alle Zeiten 
nach dem Stoß; setzt man also für ¢ ein ganzzahliges Viel- 
faches von 27/y, also etwa t=z'2r/v, wo z eine ganze Zahl 
und so groß ist, daß z-2r/» der Zeit „nach dem Stoß“ an- 
gehört, so findet man aus (58) und (61) BT SEEN, 

(62) A’ sin A sin 8, 
(63) cos =v A cos + ges), 


Aus der zweiten dieser läßt sich 
X, eliminieren. Durch Integration der Gleichung (4) nach 
der Zeit findet man nämlich unter Berücksichtigung von (55) 
(64) Meg te Mg = Xo T ’ 
so daß man aus (63) erhält 
(65) m, (vt)v A’ cos = m,x%(vt)v A cos mi. 

Sehen wir also von der StoBperiode ganz ab und ver- 
gleichen nur die Werte der Variablen vor und nach dem StoBe, 
80 gestatten die Gleichungen (62) und (65) zusammen mit dem 
Energieintegral 


die Werte A’, ©, 2,’ aus A, ®, z, zu berechnen, sobald r als 
bekannt angenommen wird. 


1) Man vgl. etwa die Gleichung (58) mit (11). 7 2 a. om ae 
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Zugleich haben wir mit den linken oder rechten Seite 
von (62) und (65), wie die Form dieser Gleichungen zeigt, die 
gesuchten ‚Stoßinvarianten‘‘ gefunden, sofern wiederum yo 
der Stoßperiode selbst abgesehen wird. Man erkennt ohm 
weiteres, daß (65) im Grenzfall y= 0 in das Schwerpunkts 
integral übergeht, denn es ist lim xy (vr) =1, und, wie ay 
vt=0 

(60) hervorgeht, » A cos # die konstante Geschwindigkeit, in 
die #, im Grenzfall » = 0 übergeht. Da x (v r) mit wachsendem» 
über alle Grenzen wächst, kann man die Invariante (65) w 
auffassen, als ob die für die Berechnung der Schwerpunkts 
bewegung maßgebende Masse m,x(»r) mit zunehmende 
elastischer Bindung unbegrenzt zunimmt. 

Die Form der Gleichungen (62) und (65) ist für da 
Folgende nicht brauchbar; wir quadrieren darum noch beide 
Gleichungen und erhalten dann durch Addition folgende 
ebenfalls invarianten Ausdruck 


(67) I=m,? v? 4? (vr) A?+2m, (v t) cos + m,7 


Diese Invariante geht für den Grenzfall »=0 in da 
Quadrat des Schwerpunktsintegrals und im Grenzfall » =a, 
in Übereinstimmung mit dem früher Abgeleiteten, in die mit 
2m, x” (vr) multiplizierte asymptotische Invariante 


3) 


Wir wollen noch für die Funktion X (t) (Gleichung (58) 
zwei spezielle, im übrigen völlig willkürliche, Beispiele angeben. 
Nimmt man für das Kraftgesetz den unstetigen Verlauf 


t T 3? 
so ergibt sich Mir 
(70) (v r) = ie 
sin 
2 
Setzt man dagegen ary 
(71) 
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so wird 
(72) 4 (vt) 


Wir gehen jetzt zur statistischen Betrachtung über und 
gewinnen eine stationäre Zustandsverteilung in derselben 
Weise, wie in dem früher betrachteten Grenzfall. Im Interesse 
der später auszuführenden Integrationen empfiehlt es sich, 
von vornherein statt der Variablen x, und z, die Amplitude A 
und die Phasenkonstante 9 einzuführen. Es ist dann das Pro- 
dukt der Differentiale dz, dz, durch vAdAd® zu ersetzen. 


Wir setzen also für die Zahl der Paare von Atomen erster 
und zweiter Art, deren Phasenpunkt in der Phasenzelle 


— -ak— 

(74) AN, e do, 


wo E und J durch (66) bzw. (67) gegeben sind. ——— 


Es ist zu bemerken, daß die ,, Schwerpunktsinvariante (65)* 
für die Aufstellung des statistischen Gesetzes nicht brauchbar 
ist, weil sie eine von Null verschiedene mittlere Geschwindig- 
keit für die Atome m, liefern würde. Bei dem Ansatz (74) 
ergibt sich die mittlere Geschwindigkeit beider Atomsorten 
gleich Null, aber er liefert für die freien Atome nicht genau 
die Maxwellsche Verteilung. Darüber sind einige Worte zu 
sagen. 

Wir gehen dafür von dem Grenzfall » = 0 aus; in diesem 
Falle reduziert sich unser „linearer Körper‘ auf ein ,,lineares 
Gas“, das aus zwei Atomsorten besteht, mit dem Verteilungs- 
gesetz (vgl. S. 650) 


Dieses Verteilungsgesetz liefert zwar formal für die Ge- — 
schwindigkeiten x, und #, einzeln genommen das Maxwellsche 
Gesetz, wie man durch Integration nach di, bzw. dz, fest- 
stellen kann; es besteht aber der wesentliche Unterschied, 
daß d N,, nicht einfach proportional dem Produkt der Einzel- 
wahrscheinlichkeiten für die Geschwindigkeiten z, und # ist. 
Der Stoßzahlansatz ist also aufgegeben und eine Koppelung 
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der Geschwindigkeiten angenommen, wie sie von Burbury 
nachdrücklich gefordert worden ist.!) 

Bei Maxwell wird als unabhängiges Elementarereignis das 
Auftreten einer Geschwindigkeit von gegebenen Komponenten 
&, n, € angesehen; daher wird die Existenz einer Verteilung: 
funktion F (£,n,£) vorausgesetzt, derart, daß 


(76) dAN=F(£, n, 


die Anzahl der Atome von der gegebenen Geschwindigkeit ist, 
Nun ist aber das Auftreten einer bestimmten Geschwindigkeit 
kein unabhängiges Ereignis, weil das einzelne Atom kein 
dynamisch unabhängiges System ist. Es muß also F so be 
stimmt werden, daß es durch die Stöße nicht gestört wird, 
was nur mit Hilfe des Stoßzahlansatzes gelingt. 

Dem Ansatz (75) liegt dagegen die Auffassung zugrunde, 
daß als unabhängiges Elementarereignis, und damit als nieht 
mehr auflésbares Wahrscheinlichkeitselement, das Auftreten 
eines Paares von Atomen mit gegebenen Geschwindigkeiten 
zu gelten hat. Diese Auffassung ist allgemeiner und zulässiger 
als ein Ansatz von der Form (76), weil ein Paar von Atomen 
mit höherem Rechte als dynamisch unabhängiges System an- 
gesehen werden kann, als ein einzelnes Atom. In der Tat gehen 
ja bei Stößen in verdünnten Gasen immer Paare von Ge 
schwindigkeiten (vektoriell verstanden) in andere Paare über. 
Es muß aber ausdrücklich betont werden, daß auch ein Ansatz 
von der Form 


(77) AN. = (é, N1> ay) dé dn dt dé,d dd, 


nur eine Approximation darstellt, die nur dann zulässig ist, 
wenn das Gas so verdünnt ist, daß ausschließlich binäre Stöße 
berücksichtigt zu werden brauchen. Im allgemeinsten Falk 
muß als unabhängiges Wahrscheinlichkeitselement das Aut 
treten eines Zustandes gelten, bei dem sämtliche Koordinaten 
und Geschwindigkeiten des betrachteten Systems vorge 
schriebene — natürlich mit den Bedingungen des Systems ver- 
trägliche — Werte haben. Es mag im übrigen betont werden, 
daß der „Paaransatz‘‘ keineswegs wesentlich für unsere Be 
trachtungen ist. Der Ansatz (74) liefert auch eine stationär 
Verteilung, wenn man eine Gesamtheit von Systemen be 


1) 8. H. Burbury, Phil. Mag. . (5) 50. S. 584. 1900. 
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trachtet und unter dN,, die Zahl derjenigen Systeme ver- 
steht, deren sämtliche Koordinaten und Geschwindigkeiten in 
einer gegebenen Zelle des Gesamtphasenraumes liegen. Man 
muß dann nur unter dw diese Zelle, unter E die Systemenergie 
und unter J die Summe sämtlicher binärer Invarianten von 
der Form (67) verstehen. Wesentlich für unsere Betrachtungen 
ist dagegen die Beschränkung auf binäre Invarianten, und 
dafür ist der Paaransatz ein bequemer Ausdruck. 

Der Ansatz (75) liefert für die beiden Atomsorten ver- 
schiedene mittlere lebendige Kräfte. Man findet nämlich 


2 

2 ” ’ Pr 

(18) 

1 
1+ (m, + mg) 


wd der Gleichverteilungssatz wäre somit auch für translative 
Bewegungen ungültig. Da kein Grund dafür vorliegt, eine so 
weit gehende Annahme zu machen, wenigstens für Tempe- 
raturen, die nicht sehr tief sind, wollen wir im folgenden an- 
nehmen, daß f’/« eine endliche Konstante und 


(9) (m <1 


ist. Dadurch reduziert sich (75) auf das Maxwellsche Gesetz 

und die Konstantenbestimmung (35) bleibt bestehen. 
Unter der Annahme (79) vereinfacht sich das Verteilungs- 

gesetz (74) zu!) ; 


dM ace 


Durch die Anwesenheit der Atome m, wird also dann das ; 
Maxwellsche Verteilungsgesetz der freien Atome nicht gestört. — 
Es bleiben nun noch die Energiemittelwerte zu berechnen. 

Dazu ist zu bemerken, daß die Berechnung nicht durchführbar 
1) In einem Häufigkeitsgesetz von der Form ze" ?*""?4?420<xz<®) 

kann der Exponent 29x neben px? vernachlässigt werden, wenn 


2¢z<p x? ist, und das führt zur Bedingung q <p”. 
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ist, solange nicht r als Funktion der Zustandsvariablen bekannt 
ist. Entsprechend dem orientierenden Charakter der vorliegen. 
den Rechnungen wollen wir für t seinen Mittelwert gesetzt 
denken, d.h. also annehnen, daß r nur Funktion der abso. 
luten Temperatur sei. 

Indem man zuerst nach #,, dann nach A und schließlich 
nach # integriert, findet man aus (74) 7 


(1438 m; + 2. 


| ) | + 
28 

20 1 + 26° m + 

— 

Aus den beiden letzten Formeln erkennt man, daß a 


für alle Werte von yr endlich bleibt und mit wachsendem yr 
gegen 1/2« konvergiert, während Z, mit wachsendem » r gegen 0 
strebt; und zwar gilt das, wie immer die Konstanten « und f 
bestimmt werden, wenn sie nur endlich bleiben. 

Insbesondere unter der Annahme (79) erhält man aus 
(82), (83) oder unmittelbar aus (80) 


1 1 
(85) 23 


Damit ist also gezeigt, daß auf nd des Verteilungs- 
gesetzes (74) die mittlere Energie der linearen Oszillatoren mit 
unbegrenzt zunehmendem »r gegen Null konvergiert, auch 
dann, wenn über die Konstante 8’ so verfügt wird, daß für 
freie Atome das Maxwellsche Verteilungsgesetz praktisch be 
stehen bleibt. Im übrigen erweist sich der Energieinhalt det 
Oszillatoren als Funktion von vr; er geht für kleine Werte 
von »r in den Gleichverteilungswert k T über, und zwar mit 
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derselben Annäherung, mit der die Gleichverteilung für freie 
Atome gilt. 

Die drei Konstanten unserer Verteilungsfunktion C, « 
und f’ sind aus den Bedingungen des Problems zu bestimmen. 
Wesentlich für unsere Betrachtungen ist die Annahme eines 
yon Null verschiedenen ,,Koppelungsfaktors ß’‘; über seine 
Größe kann die Theorie nichts aussagen, sondern sein nume- 
rischer Betrag muß aus empirischen Daten ermittelt werden. 
Nehmen wir also an, unser „linearer Körper‘‘ bestehe aus 
N, Atomen m, und N, Atomen m,, so muß für irgendeine 
Temperatur, bei der bereits Abweichung von der Gleich- 
verteilung vorliegt, bekannt sein, wie sich die Gesamt- 
energie E auf den „festen‘‘ und den „gasförmigen‘‘ Bestand- 
teil des Systems verteilt. Bedeuten E, und E, die betreffenden 
Anteile (E = E, + E,), so bestimmen sich die drei Kon- 
stanten aus den Gleichungen 


In die rechten Seiten dieser Gleichungen sind streng ge- 
nommen die Ausdrücke (81), (82) und (83) einzusetzen. Aus 
den oben (S. 653) angegebenen Gründen muß aber ß’ der ein- 
schränkenden Bedingung (79) genügen, und wir werden dem- 
gemäß nur solche Kraftgesetze zuzulassen haben, die unter 
Wahrung der Bedingung (79) Anschluß an die Erfahrung geben. 
Dann dürfen auch die Formeln (82) und (83) durch die ein- 
facheren Ausdrücke (84) und (85) ersetzt werden, wobei « 
durch (35) bestimmt ist. 


Bei der Betrachtung des Grenzfalles »t > 1 mußten wir, 
um der Helmholtzschen Forderung genügen zu können, die 
Annahme machen, daß ß mit wachsendem » über alle Grenzen 
wächst (S. 640, 641). Wir haben jetzt gezeigt, daß auf Grund des 
Ansatzes (74) der Helmholtzschen Forderung stets und ohne 
weitere Annahme genügt wird; unsere gegenwärtige Unter- 
suchung beantwortet auch die Frage, in welcher Weise das 
früher benutzte 6 mit » ins Unendliche wächst. Denn ß£ steht 
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mit unserem jetzigen ß’, wie ein Vergleich der Ansätze (74) 


und (25) zeigt, in dem Zusammenhang ies. 
(87) 7 


Die Folgerungen über den Energieinhalt linearer Oszilla 
toren, die wir aus der Theorie gezogen haben, stehen in Uber. 
einstimmung mit der Erfahrung, wenn man die Hypothes 
einführt, daß unser yt mit dem hv/2r kT !) der Planck. 
schen Formel identifiziert werden darf.2) Damit wäre m. 
gleich eine dynamische Deutung der Planckschen Konstante} 
gegeben; es würde sich nämlich folgern lassen, daß di 
Stoßdauer umgekehrt proportional der absoluten Temperatur, 
und das Produkt der Stoßdauer r in die Energie 2xkT 
gleich dem elementaren Wirkungsquantum ist. 

Wenn auch die ausgesprochene Hypothese naheliegend ist, 
so möchten wir sie doch nur mit allem Vorbehalt aussprechen; 
ihre Zulässigkeit hängt im wesentlichen davon ab, ob es ge 
lingen wird, die hier entwickelten Gedanken auch auf anderen 
Anwendungsgebieten der Quantentheorie zur Deutung de 
Erscheinungen heranzuziehen. In diesem Zusammenhang 
sei nur daran erinnert, daß schon Sommerfeld?) die Brems 
zeit der Elektronen bei der Erzeugung der Röntgenstrahla 
mit Erfolg zu dem Planckschen Wirkungsquantum in Be 
ziehung gesetzt hat. 

Wir sehen davon ab, die Formeln, die sich für de 
Energieinhalt eines linearen Oszillators ergeben, wenn man it 
(84) etwa die speziellen Ansätze (70) oder (72) einführt, mit 
dem Planckschen Gesetz zu vergleichen; denn die Kraft 
gesetze (69) und (71), die zu (70) und (72) führen, sind völlig 
willkürlich. Man kann umgekehrt die Frage aufwerfen, wie 
die Funktion x (vr) bestimmt werden müßte, und welche 
Form des Schwerpunktsintegrals (65) danach anzunehme 
wäre, damit für den Energieinhalt des Oszillators die Planck 
sche Formel resultiert. 


1) Man vgl. die Anmerkung auf S. 632. 

2) Etwas allgemeiner wäre die Annahme, daß die beiden dimension 
losen Zahlen vr und hv/2rkT sich nur um einen numerischen Fakta 
von der Größenordnung 1 unterscheiden. Wegen der Unbestimmtheit, 
die schon der Erklärung von + anhaftet, sehen wir von der Einführung 
eines solchen Faktors ab. 
3) A. Sommerfeld, Sitzungsber. d. kgl. Bayr. Akad. 1911. &1 
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Da unsere Formel (84) für »t = 0 nicht gegen 1/a =k T 
konvergiert, sondern gegen den nur wenig davon verschiedenen 
Wert 


geiorde 
von 
89 vr= hr — <= 2 m, = (048 
gesetzt wird. Daraus ergibt sich 
(2) = 
(90) 8% 


Das so bestimmte x (z) konvergiert für z = 0 gegen 1 und 
für z = co gegen oo, wie es sein muß. Die Funktion x (2) läßt 
sich demnach in der gewünschten Weise ermitteln, wie immer & 
vorgeschrieben wird, also auch wenn e eine kleine Zahl ist, 
wie es für die Gültigkeit von (84) erforderlich war. 

Wenn x(z) bekannt ist, läßt sich die Bestimmung des 
zeitlichen Verlaufs der Stoßkraft auf eine Quadratur zurück- 
führen. Denn die Gleichung (59) stellt bei gegebenem x (2) 
tine Integralgleichung für g (s) dar, deren Lösung nach dem 
Fourierschen Theorem durch 


(91) (8) = "1 (z) cos (sz) dz 
J 


gegeben ist.) Wir möchten aber auf die so gekennzeichnete 
Bestimmung von x (z) und g (s) keinen sonderlichen Wert legen, 


weil man auch aus der Kenntnis von g (s)=g (=) noch keinen 


Rückschluß auf die beim Stoß wirkenden Kräfte als Funktion  —__ 
des Abstandes der stoßenden Teilchen ziehen kann.?) Es müßte — 


1) Man vgl. H. Weber, Die partiellen Differentialgleichungen der 
mathematischen Physik. Bd. I. $ 19. Braunschweig 1910. 

2) Es ist auch darauf hinzuweisen, daß das angegebene Verfahren 
nur dann zu einer eindeutigen Bestimmung von x (x) führt, wenn ander 
Forderung festgehalten wird, daß #’/« eine Konstante sein soll (vgl. S. 653). 
Andererseits braucht es nicht zu befremden, daß dieser Quotient in die 
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vielmehr gezeigt werden, daß man aus plausiblen Annahmen 
über die Wirkungsweise der Kräfte zu einer Funktion yfg 
gelangen kann, die sich nicht merklich von (90) unterscheidet, 
Wesentlich scheint uns dagegen der Nachweis, daß der 
Energieinhalt eines Oszillators Funktion von vr ist. Macht 
man die schon oben ausgesprochene Annahme, daß »r mit 
hv/2rk T identifiziert werden darf, so wird die Stoßdauer, 
wie auch bereits erwähnt, umgekehrt proportional mit der ab- 
soluten Temperatur. Hält man diese Folgerung mit der früher 
(S. 645) angegebenen Schätzung der Stoßdauer zusammen, 99 
sind die beiden Werte nur vereinbar, wenn man n = 8 setzt, 
d.h. also, wenn man annimmt, daß die Stoßkräfte um- 
gekehrt proportional mit dem Kubus der Entfernung ab- 
nehmen. Eine solche Folgerung ist einigermaßen überraschend, 
da man Grund zu der Annahme hat, daß die Abnahme der 
beim Stoß wirksamen Kräfte mit der Entfernung sehr viel 
schneller erfolgt. Demgegenüber möchten wir betonen, daß 
einerseits die Schätzung der Stoßdauer auf S. 645 auf sehr 
vereinfachten Annahmen und einer einigermaßen willkürlichen 
Definition des Begriffs der Stoßdauer beruht; andererseits ist 
für das r vielleicht nur die mittlere Phase des StoBvorganges 
maßgebend, bei der die Kräfte am stärksten sind. Während 
der größten Annäherung zweier stoßender Moleküle ist aber 
ein Kraftgesetz mit n = 3 durchaus nicht unmöglich, zumal 
wenn man an einen elektrischen Ursprung dieser Kräfte denkt.!) 
Für den Zahlenwert der Stoßdauer erhält man aus der 
Gleichsetzung von vr und hv/2r kT die Aussage 


(92) t = 7,59- 10-13 T-1 sec. 


Von den mit dieser Formel berechneten Stoßzeiten läßt 
sich nur so viel sagen, daß sie mit den Forderungen der 
kinetischen Gastheorie nicht im Widerspruch stehen. Fit 
T = 273 ergibt sich r = 2,78-10-™; dagegen ist z. B. die 
Zeit, die zwischen zwei Stößen vergeht, in Sauerstoff bei 
Normalumständen t’= 1,65- 10-10, Erstere Zahl ist klein gegen 
letztere, wie es die Gastheorie für ihre Gültigkeit voraussetzt; 


Bestimmung von 7 (x), und damit von g(s) eingeht [vgl. Formel (M)]; 
denn diese Funktionen stellen nach der Art, wie sie eingeführt sind, nur 
einen gewissen mitileren Funktionsverlauf dar. 

1) P. Debye, Physik. Zeitschr. 22. S. 302. 1921. 
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erst bei außerordentlich niedrigen Temperaturen würde das 
nicht mehr zutreffen. 

Zum Schluß möchten wir nochmals betonen, daß alle von _ 
uns aufgestellten Häufigkeitsansätze nur hinreichende Be- 
dingungen für Stationarität des Zustandes darstellen. Wir 
glauben aber den Beweis erbracht zu haben, daß es mit den 
allgemein anerkannten und gebrauchten Prinzipien der klassi- 
schen statistischen Mechanik möglich ist, Abweichungen von 
der Gleichverteilung und das Absinken des Energieinhalts der 
Oszillatoren auf Null zu erklären. 

Mit dem Gleichverteilungssatz fällt natürlich auch die ft 
ausgesprochene Folgerung, daß die statistische Mechanik mit hee 
Notwendigkeit zum Rayleighschen Strahlungsgesetz führt. ae 
Macht man, wie wir es getan haben, die Hypothese, daß die 
Stoßdauer umgekehrt proportional der absoluten Temperatur 
ist, und kombiniert unseren Ausdruck (84) mit der Planck- 
schen Formel, die den Energieinhalt eines monochromatischen 
Oszillators mit der Strahlungsdichte im umgebenden statio- 
nären Strahlungsfelde in Beziedung setzt, so gelangt man in 
bekannter Weise zu einer Strahlunsgformel. Die so gewonnene ~~ 
Strahlungsformel steht unter allen Umständen mit dm Wien- | 
schen Verschiebungsgesetz im Einklang und unterscheidet sich 
insbesondere beliebig wenig vom Planckschen Strahlungs- 
gesetz, wenn man die Funktion x (vr) aus der an (90) 
bestimmt. 


Zusammenfassung. 


$ 1. Es wird das statistische Gleichgewicht zwischen frei 
beweglichen und elastisch an Ruhelagen gebundenen Atomen 
behandelt, und zwar für den Grenzfall, daß die Eigenschwin- __ 
gungsdauer klein ist gegen die Stoßdauer. Die Untersuchung 
geht auf die Mechanik des Stoßes zwischen einem freien und 
einem gebundenen Atom zurück, und es wird für den be- 
trachteten Grenzfall eine Bewegungsinvariante gewonnen. 
Unter Berücksichtigung dieser Bewegungsinvariante neben der 
Energie wird eine stationäre Zustandsverteilung aufgestellt. 
Die so gewonnene Verteilung widerspricht dem Gleichver- 
teilungssatze und läßt sich mit der Forderung in Einklang 
bringen, daß der Energieinhalt der Oszillatoren mit zunehmen- __ 
der Frequenz gegen Null konvergieren muß. 
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$ 2. In derselben Weise wird das statistische Gleich. 
gewicht zwischen freien Atomen und Rotatoren in einer Ebene 
behandelt, wiederum unter Beschränkung auf den Fall, daß 
die Umdrehungsdauer klein ist gegen die Stoßdauer, Eine 
Abschätzung der letzteren zeigt, daß diese Voraussetzung für 
hinreichend tiefe Temperaturen zutrifft. In diesem Gebiete 
gilt der Gleichverteilungssatz nicht, und die Existenz einer 
Nullpunktsenergie ist mit der Theorie vereinbar. 

$ 8. Die Bedingungen des ersten Paragraphen werden 
durch die Annahme vereinfacht, daß die Stoßkräfte explizit 
als Funktion der Zeit gegeben, und daß alle Bewegungen auf 
eine Koordinatenrichtung beschränkt sein sollen (‚linearer 
Körper“). Unter diesen Vereinfachungen läßt sich die Auf 
gabe für beliebige Werte der Schwingungsfrequenz 2» und 
der Stoßzeit r behandeln. Es wird ein stationäres Verteilung: 
gesetz gewonnen, demzufolge der Energieinhalt der linearen 
Oszillatoren mit unbegrenzt wachsendem »r gegen 0 kon 
vergiert. Unter der Annahme, daß yt mit dem hv/2nkT 
der Planckschen Formel identifiziert werden darf, läßt sich 
eine unbestimmt gebliebene, durch den zeitlichen Verlauf der 
Stoßkräfte bedingte Funktion so bestimmen, daß sich die 
Plancksche Formel für den Energieinhalt der Oszillatoren 
und damit für das Strahlungsgesetz mit beliebiger An- 
näherung ergibt. 
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Leipzig, im Februar 1924. 


(Eingegangen 27. Februar 1924.) E 


es 
“ 
ER, 
= 
: 
of 7 
ur 
af 
| 


von B. Walter.) 
(Dem Andenken C. H. Winds gewidmet.) 


Se xd ‚(Hierzu Tafel I u. IL.) 


Wenn auch die der Röntgenstrahlen heute 
als sicher bewiesen angesehen werden kann, so ist es = 
aus mehreren Gründen erwünscht, die früher im Interesse des 
Nachweises jener Natur — zumal von C. H. Wind — soviel — 
gesuchte Tatsache, daß die Strahlen ebenso wie diejenigen des 
gewöhnlichen Lichtes auch beim Durchgang durch einen ein- 


zustellen; vor allem, um auch auf diesem direktesten Wege 
die absolute Größe der Wellenlänge eines bestimmten Röntgen- _ 
strahls und damit zugleich auch die absolute Größe der damit a 


stellen. 

Die unmittelbarste und sicherste Art nun, die Beugung __ 
einer Strahlung mit einem einfachen Spalt nachzuweisen, ist = er 
offenbar die, ein Bild dieses Spaltes mit Beugungsfransen, d.h. 
mit Intensitätsmaximis und -minimis neben dem eigentlichen __ 


Spaltbild, herzustellen; und dieses Ziel glaube ich nun mit = 


auch noch in sehr unvollkommener, so doch in m. E. durchaus 2 
überzeugender Weise erreicht zu haben. Ay et, 

Zu den Versuchen diente, soweit die Spalte in Frage 
kommen, die gleiche Apparatur, wie sie bereits vor 15 Jahren sy 


era 


1) Den größten Teil des Inhaltes der vorliegenden Abhandlung habe 
ich bereits am 10. Februar d. J. auf der Tagung des ee re 
Gauvereins der Deutschen Physikalischen Gesellschaft in ee 
vorgetragen; inzwischen sind aber noch einige Tatsachen hinzugekommen, 
welche meine damals noch etwas unsichere Auffassung der in meinen 
Beugungsaufnahmen in Frage kommenden Erscheinungen m. E. en. 
allen Zweifel erhoben haben. 
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zu den analogen Versuchen von R. Pohl und mir!) benutzt 
wurde; d. h. von den Strahlen der Röntgenröhre wird zunächst 
durch einen vertikalen Spalt von etwa 6 u Weite ein sehr 
schmales Bündel ausgeblendet, das zur Abbildung eines in etwa 
75 cm Abstand von dem Spalt aufgestellten, damit parallelen 
zweiten Spaltes dient, der also den eigentlichen Beugungsspalt 
darstellt. Sein Querschnitt hatte bei diesen Versuchen stets 
die Form eines Keiles, der sich bei meinen hier zu beschreiben- 
den neueren Aufnahmen von einer oberen Weite von etwa 40 u 
auf eine Länge von etwa 18 mm nach unten hin bis zu voll. 
ständigem Schlusse verjüngte. 

Hinsichtlich des Aussehens sowie auch der Ausrichtung 
der beiden Spalte muß ich auf die genannte Abhandlung ver- 
weisen. 

Die photographische Platte, auf welcher der Beugungsspalt 
abgebildet wurde, befand sich wieder wie früher in ebenfalls 
75cm Abstand hinter diesem Spalt. Es wurden stets Hauffsche 
Ultrarapidplatten benutzt, die in ziemlich viel gebrauchter Hydro- 
chinonmethollösung solange entwickelt wurden, bis in ihren 
nicht belichteten Teilen im roten Licht der Dunkelkammer 
deutlich die ersten Spuren des Schleiers zu erkennen waren. 
Ich führe dies an, weil nach meinen Beobachtungen der Grad 
der Entwicklung der Platte für die möglichst deutliche Hervor- 
rufung des gesuchten Phänomens von ziemlicher Bedeutung zu 
sein scheint. 

Ein grundlegender Unterschied zwischen meinen neueren 
Beugungsversuchen und den älteren bestand nun aber in der 
Art der angewandten Röntgenröhre. Es ist nämlich klar, daß 
man zur Erzielung möglichst deutlicher Beugungsfransen mit 
einer möglichst monochromatischen Strahlung arbeiten muß. 
Die früher benutzten Röhren, in denen die Antikathode aus 
Platin bestand, genügten nun dieser Bedingung in der denk- 
bar schlechtesten Weise, da das K-Linienspektrum dieses Me 
talls bei den zum Betriebe dieser Röhren benutzten Spannungen 
noch nicht angeregt wurde, das hierbei angeregte L-Spektrum 
aber wieder eine so weiche Strahlung darstellt, daß sie zum 


1) B. Walter u. R. Pohl, Ann. d. Phys. 29. S. 331. 1909 (vgl. be 
sonders die Figg.1, 2 und 8 daselbst)) 
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größten Teil von der Glaswand der Röhre absorbiert wurde — SER 
abgesehen davon, daß dieses Spektrum selbst aus einer größeren = 
Zahl von Linien mit so verschiedener Wellenlänge besteht, daB 
hier schon deswegen von einer monochromatischen Strahlung 
keine Rede mehr sein kann. Die Strahlung jener Röhren be- 
stand daher in der Hauptsache aus dem von ihnen aus- 
gesandten kontinuierlichen Spektrum, d. h. aus einem Genie 2 


sur Bromabsorptionsgrenze hin, also von etwa 0,2 bis 0,92- 10-%cem; 
und zwar in einer Intensitätsverteilung, wie sie annähernd auch 2 
in dem kontinuierlichen Spektrum der sogleich näher zu er- ~~ 
örternden Fig. 1 der Taf. I auftritt. Von allen diesen Wellen 

wirkte also derjenige Teil, welcher zwischen der Anregungs- —__ 
grenze und der Silberabsorptionsgrenze (A = 0,485 - 10”®cm) 

liegt, — wegen der bekannten selektiven Absorption, welche = 
die Silberatome der photographischen Schicht auf diese Strahlen aoe EI 
ausüben — ganz erheblich viel kräftiger als der Teil mit längeren : 
Wellen, so daß man daher auch sagen kann, daß der Durch- _ 
schnittswert der bei den früheren Beugungsversuchen haupt- — 
sächlich wirksamen Wellenlängen etwa 0,35-10—* cm betrug. 
Ich erwähne dies, weil A. Sommerfeld!) bei seiner bekannten, 
auf die photometrische Ausmessung einer unserer älteren Bu- 
gungsaufnahmen durch P. P. Koch?) gegründeten ausführicen 
theoretischen Untersuchung für die wirksame Wellenlänge bei _ 
jener Aufnahme eine der obigen sehr naheliegende Größe, näm- er Rp 
lich 0,4.10® cm, d.h. tatsächlich die richtige Größenordnung 
gefunden hat. ae ist aber eine derartige Untersuchung 


die Herstellung v von Spaltaufnahmen mit Beugungsfransen in % er 

höchstem Grade wünschenswert. ee 
Zur Erzielung derartiger Aufnahmen benutzte ich nun ~ 

zunächst eine Röntgenröhre mit Molybddnantikathode*), da dieses 


1) A. Sommerfeld, Ann. d. Phys. 38. S. 473. 1912. 

2) P.P. Koch, ebenda $. 507. 

8) Die Röntgenröhrenfabrik von C. H. F. Müller, Hamburg, hatte 
die große Freundlichkeit, uns sowohl von dieser wie auch von der später 
noch zu beschreibenden Röhrenart mit Kupferantikathode sogar mehrere 
Exemplare gratis zu liefern, wofür ich derselben auch hier unsern besten 
Dank sagen möchte. 
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Metall, wie das in der Fig. 1, Taf. I, dargestellte Spektrum der 
Röhre zeigt, neben dem bereits beschriebenen kontinuierlichen 
Spektrum ein sehr starkes und auch für diese Verhältnisse ge. 
nügend monochromatisches Linienspektrum (2, = 0,71- 107° em 
i, = 0,63. 10”® cm) aussendet, dessen Strahlen ferner auch noch 
ein so großes Durchdringungsvermögen haben, daß man hier. 
bei eine von der Pumpe abgeschmolzene Röhre mit gewöhn. 
licher Glaswand benutzen kann. 


An den mit einer solchen Röhre hergestellten Spaltauf- 
nahmen konnte ich nun aber gegenüber unseren älteren Auf. 
nahmen unter dem Mikroskop keine Unterschiede erkennen; 
und so beschloß ich, zu größeren Wellenlängen überzugehen, 
Hier kam nun zunächst eine Röhre mit Kupferantikathode in 
Betracht, deren Linienspektrum hauptsächlich aus den Strahlen 
2, = 1,54-10-% cm und A, = 1,39 - 10”® cm besteht und also 
bereits mehr als die doppelte Wellenlänge von dem des Molyb- 
däns besitzt. Dieses Kupferspektrum stellt nun jedoch schon 
eine so weiche Strahlung dar, daß zunächst die Glaswand der 
Röhre hierbei mit einem Fenster aus einem durchlässigeren 
Material versehen werden muß, als welches Zindemannglas ge 
wählt wurde. Diese Substanz läßt sich nämlich mit dem Röhren- 
glas verschmelzen, so daß auch diese Röhren nach Herstellung 
des richtigen Vakuums von der Pumpe abgeschmolzen werden 
können. Sämtliche Röhren waren nach dem Coolidgetypus ge 
baut, und ihr aus Wolfram bestehender Kathodenglühdraht 
hatte die Form einer linearen Spirale, die von einer kasten- 
artigen „Sammelvorrichtung“ aus Nickelblech umgeben war 
(Modell „Media“ der genannten Röhrenfabrik). Die Antikathoden 
der Röhren waren in bekannter Weise zu einem Wassergefäß 
ausgebildet und es konnte deshalb u. a. auch aus der Menge 
des hierin verkochten Wassers ein ziemlich guter Rückschluß 
auf die geleistete Röntgenenergie gemacht werden. 


Die Röhre befand sich während der Aufnahme in einer 
größeren Kiste aus Bleigummi, die sich in Höhen- und Seiten- 
richtungen bequem verstellen ließ und nach der Seite der Spalte 
za eine beliebig weit zu schließende viereckige Öffnung hatte, 
die aus zwei gegeneinander verstellbaren Paaren von 4 mm 
dicken Bleiplatten bestand. Der Betrieb der Röhre erfolgte 
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mit Hochspannungsgleichrichter, die Heizung ihres Kathoden- 
glühdrahtes durch einen Heiztransformator. 

Die Anwendung einer Röhre mit Kupferantikathode hatte 
nun gegenüber der mit Molybdän außer der mehr als doppelt 
so großen Wellenlänge ihres Linienspektrums für unseren Zweck 
noch den weiteren Vorteil, daß zur genügend starken Anregung 
dieses Spektrums die Spannung so niedrig genommen werden 
konnte, daß dabei derjenige Teil des kontinuierlichen Spek- 
trams, dessen Wellenlängen kürzer als die der Silberabsorptions- — 
grenze sind, und der, wie wir gesehen haben, dann stets den 
stärksten Teil des kontinuierlichen Spektrums bildet, überhaupt 
nicht mehr angeregt wurde. Ein solches Spektrum meiner 
Cu-Röhre, das also in dem kurzwelligen Teile seines kontinuier- 
lichen Spektrums kurz vor der Silberbandkante abbricht, ist 
in der Fig. 2,I der Taf. I dargestellt; und man sieht nun tat- 
sichlich, daß hier das kontinuierliche Spektrum im Vergleich 
zu dem Linienspektrum schon erheblich mehr zurücktritt als 
in dem Molybdänspektrum der Fig. 1 der Taf. I. 

In noch etwas höherem Maße erreicht man dieses Ziel 
ferner dadurch, daß man irgendwo in den Gang der Strahlen — 
ein dünnes Nickelblech einschaltet, eine MaBregel, die übrigens 
noch den zweiten Vorteil mit sich bringt, daß dadurch, wie 
das Spektrum II der Fig. 2, Taf. I, zeigt, zugleich auch die 
ß-Linie des Cu-Linienspektrums ganz erheblich viel stärker ge- 
schwächt wird als die «-Linie desselben, so daß man also auf 
diese Weise ein fast streng monochromatisches Strahlenbündel 
erhält. Die Erklärung dieses übrigens bereits mehrfach an- 
gewandten Kunstgriffs ist bekanntlich die, daß die Absorptions- _ 
grenze des Nickels zwischen den beiden genannten Cu-Linien 
liegt und daher die kurzwelligere von ihnen in dem Blech die 
Eigenstrahlung desselben anregt, die langwelligere aber nicht, 
die erstere also erheblich stärker absorbiert wird als die letztere. 
Die Dicke des Nickelblechs betrug bei meinen Versuchen rund 
0,01 mm. 

Mit der dem Spektrum II der Fig. 2 der Taf. I entsprechen 
den Strahlung wurden nun meine sämtlichen neueren Beugungs- 
aufnahmen gemacht, d. h. es handelte sich also hierbei im wesent- 
lichen um die Wellenlänge 4 = 1,54 - 10~* cm der K -Linie des 
Cu, oder genauer um ein Gemisch der K_-Linie (i. = 1,537 10° cm 
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und der etwas schwächeren K '-Linie (A = 1,541- 107° 
Metalls. 

Eine derartige weiche Strahlung würde nun aber, wem 
man sie wie bei den früheren Beugungsversuchen ohne weiteres 
durch die Spalte hätte schicken wollen, in der dabei zu durch 
setzenden Luftstrecke von über 150 cm Länge eine sehr be 
trächtliche Absorption erfahren haben — nach den Absorption. 
formeln von Glocker?) beträgt z. B. für einen Weg von 140cm 
Länge in einem Gasgemisch von der Zusammensetzung der 
atmosphärischen Luft die von der Strahlung A = 1,54-10-%cm 
durchgelassene Intensität nur etwa 13 Proz. der auffallenden -, 
so daß es daher sehr wünschenswert erschien, hierbei die Luft 
nach Möglichkeit aus dem Strahlengang zu entfernen. Die 
hätte sich z. B. dadurch erreichen lassen, daß man die ganz 
Eisenschiene, auf der die beiden Spalten s, und s,, sowie auch 
die photographische Platte befestigt werden (vgl. die Fig. 1 der 
angeführten Abhandlung von Walter und Pohl) in ein zu 
evakuierendes Gefäß brachte; wesentlich einfacher war jedoch 
der von mir zu diesem Zweck eingeschlagene Weg, der daria 
bestand, sowohl zwischen den beiden Spalten s, und s,, wie auch 
zwischen s, und der photographischen Platte je eine evakuierte 
Röhre von etwa 70 cm Länge einzuschalten, deren jede beiderseits 
mit einem für Röntgenstrahlen leicht durchlässigem Fenster ver 
sehen war. Die Röhren hatten etwa 5cm inneren Durchmesser, 
bestanden aus Glas und waren an ihren Enden mit nach außen 
umgebogenen Wülsten versehen, welche abgeschliffen wurden, 
und auf welche dann mit Pizein je eine Messingplatte auf 
gekittet wurde, in die eine Öffnung von etwa 5 cm Höhe und 
0,5 cm Breite eingeschnitten war. Auf diese Öffnungen wurden 
dann ebenfalls mit Pizein 0,11 mm dicke Zelluloidfolien aut 
gekittet. Die beiden Vakuumröhren waren ferner mit je einem 
seitlichen Ansatzrohr versehen, das zur Verbindung mit der 
Luftpumpe diente und einen Hahn zum Absperren, sowie eine 
kleine Geißlersche Röhre zur Prüfung des Vakuums trug. 

Die Röhren wurden vor Beginn der Versuche auf Kathoden- 
strahlenvakuum evakuiert, ihre Luftleere ging aber stets schon 
innerhalb 24 Stunden auf Geißlersches Vakuum zurück und 
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hielt sich dann monatelang anscheinend unverändert in diesem 
Stadium, ein Beweis, daß die Gasentwicklung in den Röhren 
nicht etwa von einer Undichtigkeit der Kittungen herrührte. 
Sie dürfte vielmehr aus dem Zelluloid der Fenster stammen. 
Auch die Farbe des Lichtes der Prüfröhren beim Durchgang 
der Ströme eines kleinen Induktors sprach in diesem Sinne. 

Um aber ferner auf alle Fälle den Mindestgewinn an Strah- 
lungsenergie, der durch die Einschaltung der beiden Röhren in 
den Strahlungsgang erreicht wurde, festzustellen, wurde damit — 
etwa zwei Monate nach ihrem letzten Evakuieren — unter Be- 
nutzung der in Frage kommenden Cu-Strahlung zunächst eine 
Aufnahme des Beugungsspaltes und seiner Umgebung mit be- 
simmter Expositionszeit gemacht; dann wurden die beiden 
Röhren entfernt und nun auf der gleichen Platte dasselbe 
Objekt unter Anwendung verschiedener, passend abgestufter 
Expositionszeiten mehrmals nebeneinander aufgenommen. Auf 
der fertig entwickelten Platte wurde dann die Schwärzung einer 
bestimmten Stelle des aufgenommenen Objekts in den ver- 
schiedenen Aufnahmen mit dem Registrierphotometer verglichen. 
Dabei ergab sich, daß die Expositionszeiten, welche man zur 
Erreichung der gleichen Plattenschwärzung mit den beiden 
Vakuumröhren einerseits und ohne sie andererseits nötig hatte, 
sich wie 1:1,63 verhielten, so daß also durch die Anwendung 
der Röhren immerhin ein erheblicher Gewinn an Schwärzung 
erreicht wurde, wenn auch freilich lange nicht in dem Maße, 
wie ich erwartet hatte, so daß ich in Zukunft ein anderes 
Fenstermaterial als Zelluloid zu verwenden gedenke. Meine bis- 
herigen, nunmehr näher zu beschreibenden Beugungsaufnahmen 
mit der genannten Cu-Strahlung sind jedoch noch sämtlich mit 
Benutzung der beiden soeben beschriebenen Vakuumröhren her- _ 
gestellt, so daß also die Strahlung dabei außer dem Lindemann- 
glasfenster der Röntgenröhre und der Nickelfolie von 0,01 mm 
Dicke auch noch die vier Zelluloidfolien der beiden Vakuum- 
führen zu durchsetzen hatte, von denen die erste sich dicht 
hinter dem ersten Spalt, die zweite in einigen Zentimeter Ab- 
stand vor dem zweiten, die dritte dicht hinter diesem und die 
vierte endlich in einigen Zentimeter Abstand vor der photo- 
graphischen Platte befand. Diese war ferner in eine doppelte 
lage schwarzen Papiers eingehüllt. Die von den Strahlen 


: 
& 
arch. 
be 
10 
der 
om 
n R 
Luft 
Dies > 3 
j 
zanze 
auch 
| der 
n zu 
doch 
daria 
uierte 
‚rseils 
ver 
esser, 
außen 
ırden, 
auf 
> und 
K 
urden 
auf 
* 
t der 
eine 
ug. 
. 
oden 
schon 
und 
if 


B. Walter. 


668 


außerdem noch zu durchquerende Luftstrecke hatte eine Länge 
von im ganzen 19 cm, wovon 9 cm auf den Raum zwischen 
dem Fenster der Röntgenröhre und dem ersten Spalt, 5m 
auf den zwischen diesem und dem zweiten Spalt und 5 cm auf 
den zwischen diesem und der photographischen Platte fielen, | 

Schon die erste meiner derartigen Beugungsaufnahmen nun, 
die bei einer Belastung der Röhre mit etwa 25000 Volt Maximal- 
spannung (aus der Grenzwellenlänge des Spektrums der Fig.2,] 
Taf. I nach der Einsteinschen Formel berechnet) und 20 Milli 
ampere Strom 7 Stunden dauerte, und bei der aus dem Anti. 
kathodengefäß der Röhre ungefähr 41 Wasser verkocht wurden, 
zeigte unter dem Mikroskop an einer bestimmten Stelle des 
Bildes des Beugungsspaltes, da wo dieser eine Weite von 4—5 4 


hatte, auf der einen Seite schwache Andeutungen eines yon 


ganzen Verlauf sehr wohl als eine Beugungsfranse aufgefaßt 
werden konnte, zumal da sonst auf der ganzen Länge des Spalt- 
bildes nichts derartiges an demselben zu sehen war. Eine 
Wiederholung der Aufnahme mit 12stündiger Expositionszeit 
(verkochte Wassermenge 71) und etwas anderer Einstellung der 
Vakuumröhren, aber unter sonst gleichen Bedingungen, ergab 
die gleichen Erscheinungen mit wesentlich größerer Deutlichkeit, 
während bei einer 30stündigen Exposition (verkochte Wasser- 
menge 18 ]) das Hauptbild sich an der genannten Stelle nach 
den Seiten hin schon so stark verbreiterte, daß die Streifen 
hier wieder undeutlicher wurden. Dafür zeigten sich aber 
wieder in dieser letzten Aufnahme an einer Stelle des Spalt 
bildes, die einer etwas weiter nach unten, d. h. nach der spitzen 
Seite des Beugungsspaltes zu gelegenen Stelle dieses Spaltes 
entsprach, deutliche Anzeichen von Beugung, wie ja auch nach 
der Theorie zu erwarten war. 

Auf der Taf. II ist deswegen außer den in den Figg. 1-3 
dargestellten sich aneinanderreihenden drei Abschnitten des 
Spaltbildes der zweiten der drei genannten Aufnahmen in Fig.4 
auch noch derselbe Abschnitt wie in Fig. 3 nochmals dar- 
gestellt, der aber hier aus der zuletzt erwähnten, 30stün- 
digen Aufnahme stammt. Die vier Figuren sind durch Mikro- 
photographie der betreffenden Teile der genannten Aufnahmen 
gewonnen und stellen 76fache Vergrößerungen dieser Teile dar, 


dem Hauptbilde deutlich abgesetzten Streifens, der nach seinem | 
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| älteren angeführten Abhandlung in Frage kommenden Auf- 
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d.h. fast die gleiche Vergrößerung wie die in der Taf. VI der 
genannten Abhandlung von Pohl und mir wiedergegebenen 
Bilder, mit denen sie also unmittelbar verglichen werden 
können. Es kommen jedoch hierfür nur die beiden, in den 
Figg. 8 und 4 unserer damaligen Tafel dargestellten Ab- 
schnitte des Beugungsbildes in Frage, wie aus den über den 
Figuren angegebenen Werten der entsprechenden Weite s, des 
Beugungsspaltes hervorgeht. 

Aus diesem Vergleich ersieht man nun zunächst, dab bei 
meinen neueren Aufnahmen die Schwärzung in dem mittleren 
Hauptstreifen ganz außerordentlich viel größer ist als an den 
entsprechenden Stellen unserer älteren, was natürlich eines- 
teils auf die größere Expositionszeit bei den neueren Auf- 
nahmen und andernteils vor allem auf die wesentlich höhere 
Belastungsmöglichkeit der modernen Röhren zurückzuführen 
ist, wie am besten daraus hervorgeht, daß bei der in unserer 


nahme während der siebenstündigen Expositionszeit aus dem 
AntikathodengefiB der Röhre nur 470 ccm, bei der in den 
Figg. 1—3 meiner jetzigen Taf. II z. T. wiedergegebenen da- 
gegen in 12 Stunden ungefähr 7 1 Wasser verkocht wurden. 

Eine so starke Schwärzung des Mittelstreifens der Auf- 
nahme, wie sie die neueren Bilder zeigen, ist nun aber auch 
nötig, wenn man auf das Auftreten von Beugungsfransen neben 
jenem Streifen rechnen will; denn nach der Fresnelschen 
Theorie beträgt die Intensität im ersten Beugungsmaximum 
mur 4,5 Proz. und die im zweiten sogar nur 1,6 Proz. von der 
im Mittelstreifen. 

Was nun aber die gesuchten Beugungsfransen in den in 
der Taf. II dargestellten Abschnitten meiner in Rede stehenden 
Aufnahme anbetrifft, so sieht man schon beim einfachen Drauf- 
blicken auf die vor einem auf dem Tische liegenden Figuren, 
daß sich zunächst in der Fig. 1 zumal auf der rechten Seite 
des Mittelstreifens von diesem von oben nach unten zu all- 
mählich ein Streifen abzulösen beginnt, der dann in der sich 
daran anschließenden Fig. 2 ganz deutlich als gesonderter 
Streifen neben dem Hauptbilde daherläuft, um sich in der 
einem noch engeren Teil des Beugungsspaltes entsprechenden 
Fig. 3 noch weiter vom Hauptbilde zu entfernen und sich 
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hier allmählich in dem Plattenschleier zu verlieren: das letztere 
offenbar, weil hier die Intensität der durch den immer enger 
werdenden Beugungsspalt hindurchgedrangenen Strahlung schon 


B. Walter. 


zu gering wurde. In der mit 2%/,facher Expositionszeit her- 
gestellten, in Fig. 4 der Tafel wiedergegebenen Aufnahme des. 
selben Teiles des Spaltes, sieht man nämlich den Streifen an 
dieser Stelle wieder deutlicher als in Fig. 3, wie man besonders 
dann erkennt, wenn man die 
Figuren in der in der neben- 
stehenden Textfig. 1 dargestellten 
Weise betrachtet, d. h. also in der 


Aah 


Licht a Richtung des Mittelstreifens mög- 
bid schräg über das betreffende 
Bild hinwegsieht und dabei die 


walten Tafel so hält, daß hinter dem 
Bilde in der Blickrichtung ein 
dunkler Hintergrund liegt, da sonst 


Fig. 1. der vom Papier zurückgeworfene 
ays Glanz stért. 

Diese Betrachtungsweise der Bilder empfiehlt sich nun 
aber besonders auch fiir die beiden Figg. 1 u. 2 der Taf. Il; 
denn einesteils tritt dann der Verlauf des oben beschriebents 
Beugungsstreifens derselben noch mit viel gréBerer Deutlich- 
keit hervor; vor allem aber deswegen, weil man nur so die 
hier geradezu fundamentale Tatsache beobachten kann, daß 
rechts neben dem genannten Streifen in der ganzen Länge der 
Fig. 1 noch deutlich ein zweiter, schwächerer Streifen einher- 
läuft, der sowohl seinem Abstand wie auch seiner Intensität 
nach durchaus dem Streifen zweiter Ordnung entspricht, Tat- 
sachen, die m. E. die Richtigkeit der hier vorgetragenen Auf. 
fassung dieser Bilder über jeden Zweifel erheben.!) 

Hierzu kommt noch die weitere Tatsache, daß die Streifen 
auch — der Theorie entsprechend — entgegengesetzt der Keil- 
form des Beugungsspaltes verlaufen, also z. B. in der Fig.l 
der Taf. II unten einen größeren Abstand von der Mitte haben 


1) Auf meinem in der Einleitung erwähnten Vortrag in Braun 
schweig hatte ich diesen Beugungsstreifen zweiter Ordnung noch 
nicht entdeckt, weswegen ich damals meiner Sache noch nicht 80 
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als oben. Zur genaueren Feststellung dieser Tatsache kann 
hier etwa so verfahren werden, daß man sich, indem man die 
Figur in der beschriebenen Weise betrachtet, in der Verlängerung 
der seitlichen Streifen sowohl wie auch in der der Mitte des 
Mittelstreifens oben und unten je einen Bleistiftstrich macht 
und dann die betr. Abstände bestimmt. In dieser Weise er- 


gab sich nun der fragliche Abstand in Fig. 1 der Taf. II. wine 


oben unten 


für den ersten Beugungsstreifen 2,0 2,5 mm 
und „ ,, zweiten a 3,5 4,1 


” 


woraus man zunächst ersieht, daß beide Streifen unten weiter 
yon der Mitte entfernt sind als oben, wie es die Keilform des 
Bengungsspaltes verlangt, und ferner vor allem, daß auch die 
Abstände des ersten Streifens zu der des zweiten sich beider- 
seits annähernd wie 3 : 5 verhalten, wie es ja nach der Theorie 
der Fall sein muß. Nimmt man noch hinzu, daß auch die 
Intensität der beiden Streifen, soweit hier eine Schätzung 
möglich ist, durchaus mit der von der Theorie verlangten Größe, 
d.h. also mit bzw. 4,5 und 1,6 Proz. von der des Mittelstreifens 
übereinstimmt, so sind damit alle Bedingungen, die man hier- 
bei an die fragliche Erscheinung stellen kann, erfüllt. 

Der letzte und wichtigste Prüfstein derselben ist nun aller- 
dings noch die Größe der Wellenlänge A, welche sich aus dem 
Abstand dieser Streifen von der Mittellinie für die hier in Frage 
kommende Röntgenstrahlung ergibt. Die Theorie liefert hierfür 
bekanntlich die Gleichungen 
worin d, und d, bzw. die Abstände des ersten und zweiten 
Beugungsstreifens von der Mitte des Bildes, 5 den Abstand 
der photographischen Platte vom Beugungsspalt und s, die 
Weite dieses Spaltes an der betr. Stelle bedeuten. Zur Be- 
rechnung von A gehen wir nun hier am besten von dem Ab- 
stand d, des Beugungsstreifens zweiter Ordnung der Fig. 1 
der Taf. II aus, für den sich, da es sich hier um eine 76fache 


Vergrößerung des Originalbildes handelt, nach den oben an- 


gegebenen Zahlen im Mittel d, = Oe = = 0,50- 107? cm ergibt. 


Der Abstand 5 ferner ist schon ER zu 75 cm angegeben, 
und es fehlt also zur Berechnung von A nur noch die Spalt- 


tere = | 
on 
her- 
ies 
ig 
ies 
die 3 
den- 
ay 
lten 4 
der 
nog- 
ande 
die 
d 
em 
: 
onst 
| 
ene 
nun 
die 
daß 
der 
ıher- 
silät 
Tat 
Auf- 4 
: 
eifen “4 
Keil 
aben 4 
3 
ht m 
noe Rice - 


672 B. Walter. Beugungsfransen an Spaltaufnahmen usw, 


weite s, an der betr. Stelle. Diese ergab sich nun durch 
mikroskopische Ausmessung im Mittel einer größeren Reihe 
von Versuchen, die teils vor, teils nach der Aufnahme an. 
gestellt und bei denen auch Mikroskopobjektive verschiedener 
Vergrößerung benutzt wurden, zu 4,9 u, so daß sich also 


2 9850-10 9.10-+ = om 


5 75 
ergibt. 

Diese Größe ist nun allerdings um etwa 15 Proz. kleiner 
als der erwartete Wert von 1,54-10”® cm; indessen kann dieser 
Unterschied in Anbetracht der vielen bei diesen Versuchen 
auftretenden Schwierigkeiten wohl kein großes Bedenken er- 
regen; und man wird sich zunächst damit begnügen, daß 
wenigstens die Größenordnung der auf diesem Wege gefun- 
denen Wellenlänge die richtige ist. Im übrigen gedenke ich, 
die Versuche weiter fortzusetzen und hoffe auch, bald über 
bessere Ergebnisse berichten zu können. 


Zusammenfassung: Bei Aufnahme eines keilförmigen Spaltes, 
der sich von etwa 40 u Weite bis zum vollständigen Schlusse 
verjüngte, mit den Strahlen einer Röntgenröhre mit Kupfer. 
antikathode ergaben sich auf der in 75 cm Abstand hinter 
dem Spalt aufgestellten photographischen Platte in denjenigen 
Teilen des Bildes, welche einer Spaltweite von 5—2,5 u ent 
sprachen — aber auch nur in diesen — ‘neben dem Haupt 
bild deutliche Maxima und Minima, die ihrem ganzen Verlaufe 
nach als Beugungsfransen angesehen werden müssen. In der 
Gegend 5 w sieht man sogar zwei Streifen, die bzw. die Frans 
erster und zweiter Ordnung darstellen, da sie nicht bloß das 
richtige Abstandsverhältnis vom Hauptstreifen, sondern auch 
das richtige Intensitätsverhältnis zu diesem haben. Für die 
aus ihren Abständen vom Hauptstreifen berechnete absolute 
Größe der Wellenlänge der Kupferstrahlung ergibt sich der 
Wert 1,31.10”® cm, also immerhin die richtige Größenordnung, 
wenn auch noch um etwa 15 Proz. zu niedrig. Der Unter 
schied wird auf die Schwierigkeiten der Versuche geschoben 


Dur Physikal. Inst. der Universität, im März 1924 
(Eingegangen 18. März 1924.) 
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